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ABSTRACT 


GRAVITATIONAL RADIATION THEORY 

BY 

Thomas L, Wilson 

A SURVERY IS PRESENTED OF CURRENT RESEARCH IN THE THEORY OF GRAVITATIONAL 
RADIATION. IRONICALLY/ THEORETICAL WORK ON GRAVITATIONAL RADIATION FIRST 
APPEARS IN THE ELECTROMAGNETIC (VECTOR) THEORIES OF GRAVITATION - WHERE IT 
HAS SINCE BEEN FORGOTTEN . THIS IS IN PART DUE TO AN OVER ENCHANTMENT WITH 
THE MORE GENERAL TENSOR THEORY OF RADIATION/ ALTHOUGH SUCH ENCHANTMENT IS 
WARRANTO) IN THE SENSE THAT IT MAY BE THE PROPER SOLUTION TO THE PROBLEMS 
IN RADIATION THEORY. 

Special care is taken to stay away from the metaphysics of gravitation 

THEORY AND TO STICK WITH TIC MATHEMATICAL STRUCTURE OF GRAVITATIONAL RADIA- 
TION. FURTHERMORE/ THE RADIATION PROBLEM IS TREATED IN A FASHION ENTIRELY 
DIVORCED FROM OTHER PROBLEMS IN GRAVITATION. THE DEVELOPMENT PROCEEDS 
CANDIDLY THROUGH THREE POINTS OF VIEW - SCALAR/ VECTOR/ AND TENSOR RADIATION 
THEORY - AND THE CORRESPONDING RESULTS ARE STATED. 
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S.CEUTO 


Horatio: 


Brutus: 


Horatio: 


Cassius: 


Brutus: 


Horatio: 


Cassius: 


Horatio: 


Cassi us: 


(Enter Horatio, Brutus, & Cassius) 


Alas. My experiment must not work. I have been 
searching for gravitational radiation for months now 
and I have found nothing. 


If you see nothing, then there's nothing there. The 
stars must not emit gravitational radiation. 


Not necessarily. Perhaps we are looking for the wrong 
thing. 


The fault, dear Brutus, is not in the stars but in 
ourselves - that we are poor Masters of human thought. 
And not until we conquer ourselves shall we ever master 
the secrets of the Universe. 


Perhaps. Nevertheless, if there's nothing there then 
the linear tensor theory is wrong. 


That is incomprehensible. If the tensor theory is wrong 
we have no basis for an expanding Universe. 


Horatio, dear Horatio, here you are - standing upon a 
speck of dust trapped in a sunbeam - and you have the 
audacity to state that the Universe is expanding. 


Cassius, I know the Universe is expanding. I have 
measured the cosmological redshift. Indeed, that is 
why the heavens are black at night. 


Horatio, what you have discovered is a footprint in 
the sand. And behold - it is your own. 

viii 



CHAPTER I 


INTRODUCTION 


The notion of gravitational radiation is an old one. It has as its be- 
ginnings the scalar theory of Newtonian mechanics, it becomes manifest 
and explicit in the electromagnetic vector theories of gravitation, and it 
is finally extended to its most general form in the tensor theories of 
relativity. The purpose here is to consider the nature of gravitational 
radiation. Thence, it is our hope to investigate and to understand the 
linear and the nonlinear characteristics of scalar, vector, and tensor 
radiation theory. 

Early considerations of gravitational waves can be founu in the work of 
Hooke (1671) on the physical cause of gravitation and of Laplace (1802) on 
Newtonian cosmology, but not until the advent of classical electromagnetism 
does gravitational radiation take a very definite form. This is the in- 
vestigation of Maxwell (1865) into an electromagnetic theory of gravitation. 
Subsequently, theoreticians are found applying numerous features of the 
laws of classical electrodynamics to planetary orbits in an effort to ac- 
count for several anomalies in celestial mechanics. The most popular 
problem to be found is the anomalous precession of the perihelion of the 
planet Mercury which is explained by mechanisms (all of them due to Newton) 
such as solar oblateness and variations in the inverse-square law. However, 
the concepts of electrodynamics (not forseen by Newton) are found to intro- 
duce suggestions of a velocity-dependent force fashioned by Hollzmuller (1870) 
after the electron theories of Weber (1846) and Riemann (1861). Adopting 
the precessional behaviour of Mercury, for example, as an experir ental 
basis for the determination of the velocity of propagation of a retarded 
gravitational potential, one can arrive at the value of Gerber (1895): 
c - 305,500 km/sec, within 2% of the speed of light. 
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The lasting dividend of these 19** 1 Century investigations is the 
gravitational wave. The shortcomings of Newtonian mechanics thus result 
in the electromagnetic (vector) theory of gravitation, definitively stated 
by Heaviside (1893) - not to mention the work of Lorentz, Poincare', Lodge, 
Ishiwara (1914), and H.A. Wilson (1921). Vector gravitational waves are 
treated explicitly by Gans (1905, 1912) and Abraham (1911, 1912, 1913), 
who are the first to identify the annoying features of negative Poynting 
vecto-s, negative energy densities, and particle instabilities under grav- 
itational radiation reaction. But the vector theories of gravitation are 
found to be inadequate by some and consequently this is one of the justifi- 
cations for the more general theory of gravitation, the tensor theory due 
to Einstein (1915, 1916, 1918). Einstein's treatment of mass quadrupole 
radiation is the first for tensor gravitational waves. As closer study 
demonstrates, however, there is nothing intrinsic in the tensor theory which 
excludes the possibility of negative energy states discussed by Maxwell 
in 1865. 

The nonlinear nature of Einstein's tensor theory of gravitational radia- 
tion is its most difficult, provocative, and yet rewarding feature. But 
most of the success in its application has been due to the many analogies 
with electromagnetism and the numerous linearization schemes, simply be- 
cause the exact solutions of the nonlinear field equations in this theory 
are formidable, difficult to interpret, and beset with mathematical singu- 
larities which cannot be associated with radiative sources. 

There is even controversy in the tensor theory of gravitational radiation 
as to its very existence, due to its covariance and often intractable nonlin- 
earities. In the vector theory, however, there is really no such question; 
rather it is the physical interpretation of the difference between electro- 
magnetic and vector gravi Jtional waves and their resonant interactions 
(it is a unified field theory). Regardless, the concept of gravitational 
radiation now offers an exciting new vista in theoretical physics and may 
even become a definitive test of the wealth of existing theories of electro- 
magnetism and gravitation. 
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CHAPTER II 

ELECTRCMA6I€TIC THEORIES OF 6RAVI iATIONAL RADIATION 


The Scalar Theory Of Gravitational Radiation 


The concept of a gravitational wave follows from a wave ec t . .• Recall 
from the Newtonian theory of gravitation that a force £ ex- ed upor, point 
mass m at a distance r from another mass M is the equation of motion 


F 

*0 




where G is the Newtonian constant of gravitation. This relation in turn 
defines a field intensity g(x) which is a force per unit mass. 
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For a multi pv e distribution of matter p(x') it follows that the field inten- 
sity g is derivable from a scalar potential as 


g = -y<t> (?3) 


where 



Taking the divergence of Eq (2-3) one obtains the differential form of 
Gauss' law and a zero curl 


V .y- 

Vxq = O 



(2-Sa) 

(2-SW) 
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Hence, the field is irrotational. Eqs (2-5) and (2-3) together constitute 
Poisson's Equation 

V 2 $ = -4-irGP <?-«) 

for this static and stationary case of a scalar Newtonian potential <t>. 
However, the theory is as yet nonradiative. It provides only for the in- 
ductive transfer of energy between Newtonian masses, such as by the mech- 
anism of tidal friction (Appendix G & S). 

To obtain a radiative scalar theory one borrows from his understanding of 
the vector theory - which has a scalar component - and argues that a non- 
stationary scalar potential due to an oscillatory multipole distribution 
p = p 0 exp(iwt) extends Poisson's Equation to the wave equation 


□ 4 * -Hf> fc-f) 


where Qb^-c 1 ^ , where * = 4^G, and where c is the velocity of propagation 
of the scalar gravitational wave 4>. The solution of this wave equation is 
essentially (2-4), 



(l*8a> 


except that the brackets qualify the retarded solution of Lorenz (1867), the 
advanced solution of P.itz (1908), or standing waves (both). 


A multipole expansion of h in (2-8) states that 





The Newtonian (Coulomb) scalar potential, furtherr're, has a spherical har- 

1 L 

monic decomposition (equivalent to 2-8b) into the 2 L -th multipoles which 
derives from an expansion of the denominator in (2-8a) 



(Z-8c) 



(z*w) 


where 
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and where the time-dependence has been suppressed. (This spherical harmonic 
decomposition is given for the sake of comparison with the tensor decomposi- 
tion in Chapter 3). 

Arguing that the dipole p contribution cannot radiate and conserve momentum, 
the first radiative contribution 

p ~ ff <&fi V («> 

is due to the quadrupole (L=2) 


The Vector Theory Of Gravitational Radiation 

The vector theory of radiation derives from Maxwell's theory of electro- 
magnetism. One can argue against establishing an electromagnetic theory of 
gravitation on the basis that the Newtonian field equations (2-5) are irro- 
tational. They have no vector potential. 3ut from the rotation of the plan- 
ets in the solar system one can present a reasonable argument for the exis- 

* 

tence of a gravitational vector potential. Indeed, the tensor theory of 
gravitation does just this (Appendix E, Eq E-10 & E-18). 

Maxwell's Equations will be written as follows: 


V Sr * T f 

(Ml,) 

V*E. =->§At 

(MW 

VB * 0 

(»««) 


(MM) 


where u is the permeability, e is the permittivity, E is the elect, ic f-'eld 
intensity, B is the magnetic induction, T is toe current, and o is the ci srge 
density. 

Comparison of + he Newtonian field equations (2-5) with the Maxwell field 
equations (2-11) shows that the gravitational analogue of E is g. Similarly, 

* *r 

*This also derives from the notioi of the Coriolis force as well as the remark 
of Wilson (1923) that there is a relation between the rotation of mass and the 
origin of magnetic fields (see Appendix F). 
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an entire development, of a gravitational vector theory follows that of 
Maxwell. The resultant field equations (using lower case symbols for grav- 
itation) are 


* p 


vx* ^ - awit 

m 9 

(l-UAj 

O 

il 

• 

> 

CZ-Ut) 

yxb s *4i6^j 

* 

(Mil) 

where - and • are the gravitational permesbil ity and permittivity. 

and 

i * m q s 


b = VxA 

00 +* 


are derivable from scalar and vector potentials. The axial vector 
associated with the Coriolis force for example. 

The constitutive relations are 

b car. be 


(Z-Ma) 

b * jbS 

(t-MW 


These are important in the nonlinear theory of electromagnetic gravitation 
(in the fashion of nonlinear optics and Maxwell-Dirac Spinor electrodynamics). 
The Lorentz equations of motion are 


f>t * i f l-Vf - & p-»s) 

I A9 +j*fe 

There is likewise a gravitational Ohm's Law (See Appendix F) 

j = Z 9 s /» x ( il6) 

which defines a gravitational conductivity providing diffusion equations 
for gravitational vector waves as well as gravitational hydrodynamics. In 
what follows, ^ and its associated diffusion are neglected. (See Appendix F) 


The continuity condition is 
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7-j + ^| a O (Ml) 

while the gravitational Poynting vector is simply 


§« 3*b . (t - B> 


It is straightforward to determine the wave equations for g, b, A, and t. 

~ 

They are all of type (2-25). Their solutions are $ in Eq (2-8) as well as 



The quadrupole radiation is (e = ellipticity of the ellipsoid) 




In 4- vector notation the gravitational field equations are the Maxwell 
ones deriveable from the field tensor 


1 - A _ k 

yiv * y 

These field equations are (letting < = 4ttG/e) 

W + W * " 0 • 

The gravitational stress-energy tensor is 


(2- 20a) 
(2-20W) 

(2-20 
(2-22 a) 
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where the field energy density c*t®® is 
and where the Poynting vector is 

C*t i0 . (**») 

Explicitly in terms of the fields, the stress-energy tensor is 

(iuw 

Subject to the Lorentz gauge condition 
the 4-potential A u = (4>/c , A) has the wave equation 

0 fit = -*/ (MSI 

where j p = pv u obeys the continuity condition 



in Eq (2-17). Note that (2-25) contains the scalar case (2-7) from 
Newtonian mechanics. 

This, then, is the electromagnetic theory of gravitation. It is a vector 
theory whose equations of motion are (2-15), whose field equations are (2-12), 
and whose wave equation is the vector one (2-25) subject to the gauge condi- 
tion (2-24). The theory is readily quantized, consisting of Spin-1 gravitons. 
If it be desired the development can be extended to Proca's Equations, pro- 
viding for a massive graviton and a cut-off (Lapla:e-1846 ) in the gravita- 
tional field. Furthermore, if the development is assumed to satisfy the re- 
stricted theory of relativity, the entire theory is Lorentz covariant. 

Such a theory does not account for the apparent nature of gravitation that 
all matter attracts; indeed, it maintains that in the physical world there 
could well exist states of negative matter (negative mass) or antimatter which 
are repulsive. Such speculation, in fact, has been pursued by Foppl (1896, 
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1898) and Schuster (1898). However, the principle embarassment in such a 
Maxwellian vector theory of gravitational radiation is that Maxwell himself 
abandons it (In 1865). He does so because the energy density of the grav- 
itational field in Eq (2-23) is negative-definite, if the force between 
masses is to remain attractive. Maxwell accordingly invokes a universal 
background £ 0 in order to maintain a positive-definite energy £ ;, 

£ 55 

The error in Maxwell's judgement is that he has divorced the gravitational 
theory from the electromagnetic one. This is to say that the total stress- 
energy tensor is composed of the electromagnetic one T' MVJ and the gravita- 
tional one t iJV . The total stress-energy tensor must be required to remain 
positive- definite 

(T >V + V*) > O . («»*) 

which (along with the continuity equation) couples the gravitational vector 
fields with the electromagnetic ones - or else the force between masses must 
be allowed to become repulsive upon the appearance of negative- energy den- 
sities. More explicitly, Eq (2-27a) gives the coupling between electromag- 
netic and gravitational fields as 

(**»>) 

This coupling is also supplemented by the continuity condi t‘ . (2-17) for 
both fields, 

( 2 -*> 

or 

In the absence of charge (p 6 =0), for example, an electromagnetic current can 
manifest itself in (2-27d) due to the flow or rotation of mass - constituting 
a basis for the origin of magnetic fields. Further queries into the nature 
of the electromagnetic theory of gravitation are given in Appendix F. 



10 


The Tensor Theories Of Gravitational Radiation 

Recalling now the scalar and vector wave equation*, these can be readily 
generalized and one further step can be taken by introducing a tensor wave 
equation 


for gravitational waves. In this equation T MV represents the stress-energy 
tensor for all matter except those effects due to the gravitational fields 
themselves. Wave Eq (2-28) therefore represents a linear tensor wave theory. 
However, in order to make the energy of the gravitational field positive- 
definite, the auxiliary condition 

A^ >v - 0 (MS) 


must be imposed, 
from (2-28) that 


But if condition (2-29) is imposed then it must be true 



(i-») 


which is not true. The total stress-energy tensor must be divergence-free: 

(T>’ ♦ t"), v = o , (».) 

where t MV is the contribution from the gravitational fields themselves. 
Hence, the tensor wave equation must satisfy 


subject to the divergence condition (2-31). If we assume that t MV is at 
least quadratic in the gravitational field variables, then the tensor wave 
equation (2-32) becomes nonlinear. If ever quantized, it should contain 
Spin-2 gravitons. 

In sunnary, then, the scalar, vector, and tensor theories of gravitational 
radiation are respectively 
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□4> = -*/> 

(i-D 


ix-is) 

UhT= -x(T r +tO . 

tt-wj 


They all constitute gravitational radiation, subject only to appropriate 
boundary conditions of incoming and outgoing radiation at infinity. These 
are field equations, and not necessarily equations of motion. It is the" 
latter which address the meaning of radiation reaction and the energy con- 
tent of such waves. 

By far the most sophisticated theory is the tensor one, which we are now 
prepared to discuss. 



CHAPTER III 


GENERAL RELATIVISTIC THEORIES OF GRAVITATIONAL RADIATION 


Introduction 

In the development of General Relativity, Einstein's theory is confronted 
with the task of accounting for a number of accepted gravitational phenomena, 
known to most any student of natural philosophy. Amongst these are the 
deflection of light in a gravitational field as predicted by Newton in Book 
III of Opticks (1704) and calculated by Soldner (1801), the anomalous preces- 
sion of the perihelion of Mercury which preoccupies more than 100 years of 
physics, as well as the inadequacies stated by Abraham and others of the 
restricted theory of relativity for an explanation of gravitational behaviour. 

Gravitational radiation is another such example. By 1905 this is being 
explicitly addressed in the work of Gans (1905, 1912) and several years later 
by Abraham (1911, 1912, 1913), as we have mentioned. The theoretical problem 
at this time has developed to the point of demonstrating the existence of a 
gravitational instability due to gravitational radiation reaction (the radia- 
ting source gains energy - a problem which still plagues the theory). 

When Einstein does proceed with his own theory of gravitation, it is to 
borrow an idea from Harry Bateman and to "geometrize" gravitation, that is, 
to put the theory in the form of a relation between geometry and the sources 
rather than between the field and the sources. The resultant theory of grav- 
itation is described by Einstein's tensor field equations 



where R is the Ricci tensor, g is the metric tensor, R is the scalar 

UV U V 

curvature, and 




T' +t 

> »v jut 




The equation ds -g^ dX l ‘dX v might properly be called "Bateman's Equation." 
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is the total stress-energy tensor consisting of all nongravitational effects 
T’ v as well as the stress-energy pseudo-tensor components t ^ due to the 
gravitational fields themselves. The left-hand side of the field equations 
(3-1), referred to as the Einstein tensor G , is subject to the Bianchi 
identities S liv » v - 0 which implies that the right-hand side must also be 
divergence free 

* (T /u ' t V“).„ * o (»« 

The right-hand side of the field equations is made nonlinear by the presence 
of t in (3-1) and (3-2) since the latter is at least quadratic in the 
field quantities. Furthermore, the equivalence of mass and energy due to 
Poincare (1904) 

E. » Me* 

is assumed throughout, although this postulate has never been demonstrated to 
have an experimental basis for gravitational phenomena. In effect, this 
postulate depletes the mass monopole source of gravitational radiation, in 
contrast to the invariance of classical electromagnetic charge under electro- 
dynamic radiation where (3-4) has a firmer basis (Cockroft & Walton, 1932). 
The dilatation of mass 

M * = (>*■!? 

is also postulated, which presents complications in the theory of gravita- 
tional radiation reaction. It is the origin of the "slow motion" approxima- 
tion which will be discussed shortly when the linearized theory is addressed. 

Note that Einstein's cosmological term Ag /2 is not included in the field 

yv 

equations (3-1), simply because its presence revokes all of the existing 
results in the general relativistic theory of gravitational radiation. It 
will not be discussed further. 

Approximation Methods & Coordinate Conditions 

Although there do exist exact solutions of Einstein's nonlinear field equa- 
tions (3-1) (Chapter 4), the physical interpretation and meaning of such re- 
sults is difficult to determine. As a consequence, most of the theoretical 
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work on gravitational radiation has been performed with approximation 
methods which are valuable because they provide some understanding of the 
nature of the problem - although they are actually inadmissible because they 
destroy the nonlinearity of the equations. 

These approximations consist of linearization and perturbation techniques 
which make various assumptions about the nature cf the metric tensor g . For 
example, the "weak-field" approximation assumes 


3 a* * V v + y 




where h constitutes a weak (h<<g) perturbation of the background geometry 
such as the Minkowski or the Schwarzschild backgrounds. Evaluation of 


pv 


the Ricci tensor R to first order in h then transforms the field equations 

|iV yv 

(3-1) in vacuum (T = 0) into a system of linear differential equations. An- 


pv 


other method, known as the "k-approximation" or "fast motion" method. 


v£ k V‘V v ,k v ** 

assumes that the metric tensor g^ v can be represented as a power series ex- 
pansion in the parameter k, which is proportional to the gravitation constant 
G. Next, there is a "slow motion" approximation which expands g in inverse 


<•) 


<0 




uv 


(4) 


(a-8> 


powers of the speed of light c: #, 

3« s *' 1 

I***! 1 * 

3J, = fj, * f J. ♦ 

This technique is the basis of the so-called "EIH approximation" which has 
played a significant role in the development cf Ihe equations of motion in 
General Relativity and the determination of what are referred to as "post- 
Newtonian" results (the precession of the apse, the deflection of light, and 
the gravitational red-shift). The numbers in parentheses in (3-8) reflect 

-n _ o _ o 

the "order" of the approximation - namely the power of c , such as c , c . 
One last technique that must be mentioned is the "double series" approach 
which consists of an expansion of thenetric tensor g in powers cf a mass 


parameter m as well as powers of c* n like in (3-8): 


MV 




■v 
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where now the sum of the numbers or exponents in parentheses represents the 
"order" of the approximation. 

Coordinate conditions constitute the next significant aspect of linearized 
General Relativity. Just as the Lorentz gauge condition (2-24) 



is invoked in classical electrodynamics in order to obtain the vector wave 
equation (2-25) 

a (ms) 

so do similar gauge conditions exist in general relativity. Einstein (1918) 
defines the auxiliary function 

which is subject to the Einstein coordinate condition (Appendix B) 

V\ v * o <?-«> 

resulting in a tensor wave equation of the form (2-32) from the field 
equations (3-1) 

OY r =-HT'“\ 

Such gauge conditions and coordinate conditions have led to much controversy 
in General Relativity, primarily because they destroy the general covariance 
of the theory. Fictitious gravitational waves also manifest themselves 
which must be removed by coordinate transformations. Furthermore, gravita- 
tional waves can be created in the linearized theory or they can be annihi- 
lated, simply by a coordinate condition. Upon this basis, some authors 
question the very existence of gravitational radiation. 

Two more features of the linearized theory are significant: (a) Boundary 

conditions, and (b) Radiation reaction. Boundary conditions such as the 
Sonmerfeld radiation condition (outgoing radiation at infinity) are important 
because they license the theoretician to obtain whatsoever experiment implores 
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- radiation, or no radiation. Radiation reaction is important because many 
authors employ the geodesic postulate of Einstein whereby gravitational rad- 
iation does not result in radiation reaction, in the linear theory. To be 
sure, such solutions are inconsistent with the general nonlinear theory. 

It should be evident that there is a wealth of theoretical work, both con- 
troversial and thorough, which exists in the tensor theory of gravitation. 

It is now best to begin with Einstein's original derivation of gravitational 
radiation in the linearized approximation. 

Einstein's Original Derivation Of Gravitational Radiation 

The geometrical nature of the field equations (3-1) is made more transparent 
if the latter are contracted with g yv 

R * XT , 

where < = 8?rG/c^ and T = T y is the trace of the stress-energy tensor (3-2). 

u 

Substitution of this relation (3-13) for the scalar curvature R in the 
original equations (3-1) gives 

(v,4) 

This form of Einstein's field equations states more explicitly the corres- 
pondence between the structure of geometry in the Ricci tensor and the 
structure of the sources of matter implicit in the stress-energy tensor T . 

Einstein (1918) makes a "weak-field" approximation and thereby linearizes 
the Ricci tensor in (3-14), as shown in Appendix A. Employing the Einstein 
coordinate condition (Appendix B) in Eq (B-8), the first-order Ricci tensor 
(A-6) becomes (B-9) and the field equations (3-14) above simplify to 

R* 

In other words, the wave equatior and its coordinate condition (3-11) for 
the tensor theory are, to first order in h. 
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(vi« 

(>•») 


The wave equation (3-16) is still nonlinear due to the nonlinearity of 
T , as discussed in (3-2). Its solution is represented by 




IVJVQ jV 



where r' = |X - X' | and where the brackets represent retarded, advanced, 
or standing wave solutions. To these solutions can be added those of the 
homogeneous wave equation 


Oh„*o (we) 

which are plane waves. The stress-energy tensor is now linearized by making 
a "slow motion" approximation such that M»m in Eq (3-5) while T^ v - T 
is represented by t . For distances r' much greater than the extent of the 
sources, In the radiation zone, r'» r whereby (3-17) becomes 


Imposing the boundary condition that -i = 0 at r = * and using Gauss' theorem 
one can demonstrate that 


p 

Since t smc £ in a slow-motion approximation and X =ct, then (3-19) and (3-20) 

oo o 

give 



(*-«) 


It is advantageous to use the mass quadrupole tensor (the mass monopole and 
dipole make no radiative contribution) 

0„ * jm (ax'V- 
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and express h in (3-21) as follows (recalling that * = 8^G/c H ): 


yv 


V - & { y w 

This completes the determination of the tensor wave h in the linearized 


tensor theory. 


By means of a coordinate transformation one may go to the principal axes 
of the quadrupole and choose an axis, say x\ to be an axis of symmetry such 
that Q.j v = 0. All longitudinal-longitudinal (h^, h^Q, hpp) and longitudinal- 
transverse (h-|2> h-|3> h^Q, hgg) perturbations vanish under such a gauge, 
leaving only a transverse-transverse set of perturbations h^g, ^3) 

propagating at velocity c in the direction: 


L - JG 

n *3 " 5?i 




2G 


l»» s 


(s-mU 

One now wishes to determine the radiation spectrum. The gravitational 

stress-energy is given by the pseudo-tensor t yv . For the perturbations h 59 , 

1 10 ^ 
h 33’ ^23 P ro P a 9 at1 ' n 9 along X , the energy flux ct of the stress-energy 

pseudotensor (Appendix C, Eq (C-l 1 ) ) simplifies to 
which becomes for retarded solutions xl=ct 

c ' c ' # = ilr« { C + . ** 

Substitution of (3-24) into (3-26), noting that * = S^G/c 4 , gives one 

(v27) 
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Generalizing this result to an arbitrary direction of propagation n^, the 
radiated power per unit solid angle d^'is (this particular representation is 
due to Landau & Llfshitz, 1962, but their work offers nothing that is not 
contained in Einstein's paper) 


if 



(I'M) 


The total power radiated is now determined by avers. ing this energy flux 
over all solid angles. The resultant energy loss cn:e to linearized grav- 
itational radiation is thus 


d£ 

At 


• N 


45 c ^ 


CLQ** 


(v») 


The mass quadrupole and the moments of inertia are related as Q- j = ^ 6 1j"^I-j * 
Hence for a spinning rod (Q- •> ~ I& , whereby 

1 J 


At _ 32 £ T*r»* 

At " 5 C* A , 



twice the value obtair.d by Einstein (1918), as corrected by Eddington (1922). 
[Eq (3-30) is correct only for an ellipticity of unity; see Eq (5-1).] 

For the particular care of a rod spinning about the z-axis (“=$)> Einstein's 
results are demonstrated by Park (1955) to give the radiation pattern shown 
in Figure 1. Transforming Eddington's (1922) Eq (18) into spherical coor- 
dinates, the power radiated per unit solid angle in the direction n is 

““ 5 I* 20 + -as oj £-30 

Integrating this over the azimuth Park obtains the radiation pattern 
in Figure 1 : | 
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FIGURE I 



Averaged over c, the remaining part of the solid angle, (3-32) yields 
(3-30). Figure 1 also applies to the circular Kepler oroblem. as w= shall 
see in Chapter 5. 


. * 

The angular momentum of the radiation i- likewise determinable from 


the gravitational pseudo- tensor t" v in (3-25), using Eq (D-4) 
- although Einstein neglects this aspect of the problem: 

from Appendix D 

ijjx't'*- XV*l*j, JS, 

(>♦') 

Einstein's original treatment of the linearized theory also 
the plane wave solutions of 

considers 

°V*° • 

(a-te) 

Defining appropriately the Einstein polarization tensor e (Einstein calls 
it « ; Weyl and Eddington call it a ) the general plane wave solution of 

the homogeneous equation (3-18) is 


. ik* * -ikx 

y° s v* * v c 

(3-3$ 

which satisfies (3-18) provided 


* o 

(3-34) 

W 

and (3-11) if 


*/■- hs , . 

(31S) 

Under a gauge condition (B-l) where X ,lJ = X w + the polarization tensor 

t transforms as 

Vi V 


* - V + V* + ^ • 

(3-34) 


For an Einstein-Eddington-Wey I plane-wave (3-33) propagating in the X 1 
airection with k^=k0~k while k^-k^~G 


*See Hansen (1972). 
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Letting : Q = ‘ 0<J / 2 k. ^ = -^/k. 

1 1 / 2 k , 

, and Cj = ail compon 

ents of (3-37) vanish except and 

CL. 

• 23. Those that remain (^ and ^ 3 ) 


represent the physical components of the gravitational plane waves. Those 
that vanish are "fictitious" because they can be produced by mere coordinate 
transformations. The tensor plane wave (3-33) has helicities 0, + 1, +2. 

But only helicity + 2 cannot be transformed away with a coordinate trans- 
formation. Linearized gravitational waves are therefore considered to have 
Spin-2 (although they do possess "fictitious" Spin-1 components). 

Einstein's work on the theory of linearized gravitational radiation is 
characteristically definitive in a number of respects. He obtains a wave 
equation from the linear field equations and therefrom obtains retarded 
solutions in strict analogy with the radiative solutions of c^ssical electro- 
dynamics. There is no gravitational dipole radiation because he invokes 

★ 

Newton's principle of equivalence wherein gravitational and inertial mass 
are the same, and whereby a system of masses necessarily has ?k_* same grav- 
i taticnal-to-inertial mass ratio and the dipole radiation vanishes under 
conservation of linear momentum - just as it ooes in electrodynamics for a 
system of charges with identical charge-to-mass ratios. The velocity of 
propagation of the remaining quadrupole radiation is the same as that for 
electromagnetic phenomena, upon the basis of correspondence between the 
general theory and the restricted theory of relativity and Newtonian 
mechanics (However, see Appendix G) . 

Einstein's development is not definitive, however, in the sense that he 
fails to address the real nature of his theory, its nonlinearities. Although 
there does appear to exist radiation from a mechanically driven Jacobi el- 
lipsoid possessing a time-varying mass quadrupole, there is nothing yet in 
the theory that demonstrates radiation exists foi an isolated, gravitationally 
bound system such as the two-body Kepler problem. Neither has it been es- 
tablished that the decrease in energy of the radiating mass quadrupole (3-30) 
is equal to the energy carried away by the grav i Lational radiation, a problem 
intimately related to the question of radiation reaction and the equations 


In order to write the equation of motion m.r - -GMm r/r^, Newton had to 
assume rrij = m^. This form of the principle! then, Ts due to Newton. 
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of motion. These more subtle aspects of gravitational radiation theory 
remain untouched, and are left to plague the proponents of Einstein's theory. 

The EIH Approximation 

If one were considering possibilities for tensor field equation? , the.e 
would necessarily involve the total stress-energy tensor T on tne basis of 
a generalization of classical field theory. Furthermore, the conservation 
conditions T u> ’, = 0 from classical electrodynamics would necessarily gen- 
eralize to the covariant relation 

T /V .,„ * 0 . 

Hence, if one takes the field equations in the form 

G^ V = (5-39) 

then the covariant divergence of this form (3-39) guarantees by virtue of 
(3-38) that 

G r ;v = o . (s-*) 

Because (3-40) is simply the contracted form of the Bianchi identities, tht : 
the reasonable choice for G is the left-hand side of (3-1), namely 

G fV = . (*♦»> 

Th s argument, partially due to Hilbert (1915), constitutes the teleological 
basis of Einstein's nonlinear field equations (3-1). 

The surprising result of this argument is that the equations of motion 
follow from (3-38), although this seems not to have been apparent at the time. 
Einstein (1916) originally introduces the assumption of geodesic motion, but 
Einstein & Grommer (1927) later succeed in demonstrating that the motion of a 
singularity in an external field must he along the geodesics of the external 
fieil. Finally, Einstc :n, Infeld, & Hoffmann (EIH) formulate an approximation 
method by which both the gravitational fiel< J and the equations of motion for 
its sources can be calculated from the field equations. 

This EIH approximation method is essential ]y 3 CjUCiGl S*3 tionary weak-field 
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approximation which has resulted in the so-called "post-Newtonian" results 
of General Relativity (see Appendix E). It is best understood by drawing 
an analogy due to Infeld (1938) and Trautman (19S8) with the scalar wave 
theory. If we consider the "near zone" (in the region of the sources) of 
the scalar wave equation 

□ 4 S-U f> 0-9) 


and further consider that $ may be expanded in powers of c small parameter 

■ - c-\ 

4 ®Si e "4 B = 


then operation upon j with the d’AlemLertian Q results in a series of simul- 
taneous differential equations determined by setting coefficients of the same 
powers of e equal. Because of the quasi -stationary condition that 
we obtain upon equating coefficients: 


p-e*: V’4>„ = -Xf> 





These near zone equations (3-43) exhibit the salient feature of the EIH 
approximation method. The wave equation associated with any particular power 
of e n is determined in the near zone by the solutions of lower order approx- 
imations, in this case ». " . This behaviour is represented by the arrows to 
the left of the equations which demonstrate how the approximations couple 
together. There is likewise a set of such equations for the homogeneous 
case in the wave zone where * = 0. But the interesting characteristic of 
the radiation zone is that it couples back into the near zone, acting on the 
sources as radiation reaction. 


( 3 - 43 ) 


Recalling tnat in General Relativity the scalar wave functions $ are replaced 
by the metric tensor components g , the above analogy must be reconsidered 

” M v 
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from the point of view of Einstein's field equations (3-1). As in the 
scalar case (3-42), Einstein's field equations for empty space may also be 
expanded in a power series 



( 9 - 44 ) 


However, in the case of (3-44) the coefficients of e cannot in general be 
set equal to zero, as was done in (3-43). A word of caution about the lit- 
erature is likewise in order due to a certain amount of confusion in notation. 
Most of the original EIH development is in terms of > = c"^ , rather than e. 

The treatments are equivalent, except that a 1 ” = e n+ ^, for the following 
reason (see also Appendix E). From the geodesic equation of motion, radia- 
tion reaction on the particle in motion is described by 

V . -r£x*ff (.-« 

u> «“ 

The equations of motion are of order X, but the metric components g in the 
Christoffel symbol of (3-45) are of order e. The two time derivatives 
(a 0 = c‘\ = X3^) raise the order ofe by two. To make matters even a little 
more miserable, some authors also get an extra factor of c"^ because they ex- 
press energy as P°=c"^E=XE (they get radiation reaction in x^® instead of x®). 


The EIH method is summed up by expressing g^ y as in Figure 2. This is a 

representation of the expansion of g yv in (3-8) as powers of e = c""*, where 

the numbers in parenthesis reflect the order of e. It illustrates the "post- 

Newtonian" effects in General Relativity and their origin in the metric tensor 

as depicted by the EIH method. There are two types of contributions, radia- 

2 5 

tive and non-radiative. By examining, for example, the quadrupole or "post - 
Newtonian" terms connected by the line indicated, one can see that g^ and 
g Q0 are odd-powers of e while g . is even in e. For non-radiative contribu- 
tions, g Q0 and g^ are even in e while g^ is odd. By virtue of conservation 
of linear momentum and the equivalence of inertial and gravitational mass the 
dipole radiation (represented by the dotted line) of seventh-order in >■ can 
be made to disappear. It also can be annihilated by a coordinate transforma- 
tion, as demonstr.ted by Infeld and Scheidegger (1951, 1953, 1955, 1960). 
Boundary conditions likewise determine the behaviour of Figure 2. Standing 
waves, as an example, manifest themselves if all of the radiative lines cancel 
and vanish from the diagram. 
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Inf e I d (i9J») estaoiishes *he radiation terms in Figure 2, demonstrating 
that they do not contribute to the equations of motion (radiation reaction) 
up to the seventh-order in x. Hu (1947) carries the approximation on to 
9^-order, being the first to study explicitly the problem of radiation 
reaction in the two-body Droblem. Inf eld & Scheidegger (1951) investigate 
the elimination of radiation reaction lines by coordinate transformations , 
but Goldberg (1955) finds solutions which cannot be eliminated in this 
fashion. By 1958, the definitive paper of Trautman (1958) firmly estab- 
lishes the EIH formalism depicted in Figure 2 as a linear perturbation meth- 
od with which to approximate the nonlinear field equations of Einstein 
in the weak-field, slow-motion case. 

The final step in understanding the EIH method is one of re-establishing 

Figure 2 in a manner that reflects the coupling discussed in Eq (3-43). 

This is done in the schematic representation of Figure 3 (partially due to 

Thorne, 1969), which reflects the work of the Infeld school through 1960. 

As with Eq (~.-43), the arrows in Figure 3 represent how the various levels 

of the approximation couple together. In the near zone there is the same 

scalar coupling of the even orders in except that an additional scalar 

2 5 

coupling in the odd orders manifests itself at the post level. Likewise, 
the solutions of the homogeneous wave equations (far from the sources) couple 
in the radiation zone. And with proper matching - the method of EIH or 
asymptotic matching, to be discussed later in this chapter - at the boundary 
between the near and radiation zones, the solutions are consistent. 

However, the most significant aspect of the EIH method, and any other lin- 
earized approximation of the nonlinear field equations for that matter, is 
the radiation reaction which occurs in the 9 th -order (A®) or post^-Newtonian 
approximation of Figures 283. If the boundary condition at infinity is the 
Somroerfeld condition (outgoing radiation) then the energy of the radiative 
sources must be depleted by an amount equivalent to that carried off in the 
gravitational radiation, or soaked up by some nonlinear mechanism in the 
near zone. This behaviour constitutes radiation reaction. 

A demonstration of this result is first conducted by Trautman (1958) and 

2 r 

Peres (1959). Peres succeeds in carrying the EIH method to the "post - 
Newtonian" approximation (7^-order in *-) , doing so at a critical time when 
the integrity of Einstein's linearized quadrupole formula (3-29) and the very 
existence of gravitational radiation has been a controversial subject. This 
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same result has been obtained before as we have seen, by Hu (1947) who works 

L. Q 

the radiation reaction problem out to 9-order in *(* ), at the suggestion 
of Pauli at Princeton. Both Hu (1947) and Peres (1959) consider the applica- 
tion of their results to the two-body Kepler problem and ascertain that the 
★ 

system gains energy under gravitational radiation, corroborating the re- 
sults of Gans (1905) and Maxwell (1865) - although Hu does not point this out. 
That the problem is a boundary-value problem, however, is maintained by 
Peres (i960), whereupon he obtains a decrease in energy as depicted in (3-29) 
for the Sommerfeld radiation condition. For the Kepler problem, Einstein's 
quadrupole radiation formula (3-30) reduces to 


ii 

a 


32 G 4 (m.+nQ 
5 r 8 » 




which is the result obtained by Peres (1960) - and by Hu (1947), except for 
the sign. The source mass likewise decreases by this amount. 

Radiation reaction in the EIH mass quadrupole and two-body problem is 
subsequently pursued by several authors, such as Ryten (1963), Demianski & 
Infeld (1963), and Infeld & Michalska-Trautman (1966, 1969). The EIH meth- 
od, furthermore, is employed by Chandrasekhar, et aj[ (1965-1970) and Thorne, 
ejt al_ (1967-1970) in certain astrophysical applications where the earlier 
work is extended to the hydrodynamics of perfect fluids. The significance 
of all of these results is that the EIH approximation method, while account- 
ing for some of the nonlinearities of the theory, corroborates Einstein's 
linearized quadrupole radiation. 

The derivation of Ryten (1963) of the mass quadrupole radiation in the 
EIH approximation provides, for the interested reader, a detailed example of 
how the method is pursued to the 9^-order where quadrupole radiauion 
reaction appears. Her final radiation formula, 


It . G_ 

At "5/ 





(S-45) 


is identical with the linearized approximation (3-29). 


It is of interest to note that Eddington (1922, P.251) encounters this 
problem also. 
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The " Fast Motion" Approximation 

The EIH approximation is both a "weak-field" and a "slow-motion" approxima- 
tion, particularly suited to the study of such behaviour as planetary motion 
and post-Newtonian General Relativity. But such a technique is limited in 
its applications. Consequently, another method known as the "fast motion" 
or the "k-approximation" method has established itself. It does not make 
a "slow-motion" assumption. 

This procedure consists of a scheme of successive approximations, like EIH, 
but with a power series expansion of the metric tensor g in terms of a small 

y V 

parameter k 



where k is proportional to the gravitation constant G: k = GM/R, M is the 

mean particle mass, and R is the mean interparticle separation. It is not 
quasi -stationary; that is, a The energy-momentum tensor is likewise 

expanded in a power series in k. 

Einstein (1916, 1918) is the first to use it, assuming that the field is 
sufficiently weak that the nonlinear terms can be neglected. This is the 
so-called "l s ^-approximation" in k, which is again addressed by Bertotti 
(1956) and deWitt & Ging (1960). Bertotti & Plebanski (1960) then establish 
a generalized Green's function method for nonlinear field theory and develop 
the equations of motion up to the "2 nt *-approximation" in k. Their work, 
however, contains infinite self-action terms - although these are removed 
with a renormalization process by Kuhnel (1964). Bock (1957, 1959) and 
Bonnor (1958) also discuss solutions of the "2 n< *-approximation." 

More recently, the "fast motion" approximation is utilized by Havas (1957), 

Havas & Goldberg (1962), and Havas & Smith (1965) to address the question of 

the existence of gravitational radiation from a freely gravitating system, 

such as a Kepler problem. Havas, Goldberg, & Smith confine themselves to the 
s t 

'1 -approximation," investigating its particular contribution to radiation 

damping. They obtain antidamping and a gain of energy for the radiating 

two-body problem, as does Hu with the EIH method. However, their work is 

criticized by Peters (1970) as improperly neglecting the stresses in the 
★ 

system. 


*See also Lind (1972) and Peters (1972). 
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2 5 

The radiation reaction associated with the "post ' -Newtonian" EIH effects 

r rl 

appears in the "3 -approximation" which none of these authors has addressed. 
Furthermore, Infeld (1961) has criticized the "k-approximation" method as a 
"step backwards," partly because it does not give the EIH results, which are 
well established, in the limit of low velocity. Nevertheless, the "fast 
motion" approximation will always have theoretical appeal because it has 
many realistic astrophysical applications. 


The " Double Series " Approximation 


In connection with the "fast motion" approximation just discussed, there is 
another method which is very similar - due to Bonnor (1959) and employed by 
Carmeli (1964, 1965). Using the "double series" expansion of the metric 
tensor introduced by Bonnor (1959) 



where M and m are the inertial masses of a two particle system, Carmeli also 
imposes the "slow-motion" assumption in order to expand the "k-approximation" 
terms of the metric ( §J) as a power series in c~ 

(pO <pV 

(•» 0> tf> 

just as in the EIH approximation. The result is effectively a synthesis of 
both methods with a proper correspondence between the two, thus addressing 
the criticism of Infeld (1961). Carmeli obtains the equations of motion in 
the "3 r< ^-approximation" which are to the 9^-order in c~^ (x®), where the 
radiation reaction manifests itself. Although he does not solve the re- 
sultant equations, he does present an argument maintaining that the radia- 
ting two-body problem is damped, losing energy, under an outgoing radiation 
condition. 


The Method Of Mat ched Asymptotic Expansions 

A variation of the EIH approximation has been presented by Burke & Thorne 
(1969, 1970, 1971) based upon the methods of singular perturbation theory for 
nonlinear differential equations. Their technique establishes no new results 
in gravitational radiation theory, but it does succeed in reproducing the 
Jesuits of the EIH method for gravitational radiation from the two-body 
problem found by Peres (1960) and Ryteri (1963). The work of Burke & Thorne 
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borrows extensively from the EIH met 3d, and the paper of Infeld (1938). 

The formalism is the method of matched asymptotic expansions, an under- 
standing of which follows from the continued analogy with radiation in 
electromagnetism originally due to Infeld (1938) and Trautman (1958). It 
begins with exactly the same EIH representation of the metric tensor g ^ 

as in (3-8) or Figure 2, 

a = -I t e(-H+tH‘,)+ • 

-e*V + ••• (HO 

* * *» + 

except that the various components have been expressed explicitly as certain 
functions: the are scalars, the g„, are vectors, and the g . are ten- 
sors. (These choices follow by analogy with the components of the energy- 
momentum tensor where T 00 is the scalar energy densiry, T oa is the vector 
momentum flux, and T ab is the Maxwell stress tensor.) Plugging (3-46) into 
the Einstein tensor 6 in (3-41) gives the follc«ing EIH result that 


(?-«) 


G.. * -±«‘VN> + - 

Goa • t«*KV (H^U ♦••• 

Gab = e 

where H is a 3-dyadic. J J 

Burke (1971) then imposes the auxiliary condi tons 

W ♦ V* * ° &•«) 

H + *° 

which succeed/simplifying the Einstein tensor considerably. The Einstein 
field equations (3-1) or (3-39) 






y*V * zX J (Ms> 

« lX l § **§1 , 

T 


now reduce to 
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where 0 ■» T** 

J a * -T„ ^-50) 

Sab s Tab 

defining ^ as tne energy density, d g as the momentum flux, $ ab as the matter 
stress, and G S ab as the gravitational stress. 

Figure 2 for the EIH approximation becomes in this method of matched 
asymptotic expansions Figure 4. Due to the auxiliary conditions (3-4A'< 
and the EIH field equations (3-49), the sources must satisfy 

+- d.p 

V- (§♦*§)*>. I 

as conservation laws. 

The gravitational force is determined by 

-T g S * 



where ip= 44) , with * representing the Newtonian potential. 


In effect, gravitationally bound systems such as the two-body problem are 
found to raciate quadrupole radiation while creating additional resistive 
fields which couple back from the radiation zone into the near or induction 
zone as pictured in Figure 3, causing radiation damping and a loss of energy. 

The method of matched asymptotic expansions and singular perturbation theory 
in effect criticizes the EIH method as being strictly valid only In the near 
or inner zone, due to the fact that "slow-motion" expansions are not nec- 
essarily valid at distances far from the cources in nonlinear theories. A 
separate outer zone expansion is required there where the stresses are the 
nonlinear ones due to the gravitational wave? themselves. The inner and 

outer expansions must then be properly matched. Although this technique has 

EIH 

only succeeded in duplicating the results of thg' method at the present time, 
its authors have argued that it simplifies calculations while providing a 
consistent and systematic frame work for the EIH theory. It could prove to 
be a significant step in the direction of a solution of some of the more 
formidable problems that remain unscathed ty existing techniques. 



33 



Forett Forets Forets Reaction 

, IRE 4 



34 


The Regge - Whe n ler-Zeri Hi Formalism 

Another interesting perturbation technique for solving Einstein's nonlinear 
field equations appears in the very important pape 1 by Regge & Wheeler (1957). 
They address the question of the stability of the Schwarzschild (1916) metric 
by considering a nonspherical perturbation of the Schwarzschild back- 
ground n. , which is defined by 

. JXW = V r 2 (*S«S»tf) 

t = , ms GM/f . 

Based upon a v;eak-field approximation as before in (3-6) 


= V * Vv . ( 3 *‘> 

the perturbations h are decomposed into a generalized nonsphe*' cal form 
involving tensor spherical harmonics and comprising a separabi solution 
of the type 

_m -ikt 

h M>.t) * Hfr)X (**>» • <~ 3S *) 

JkV JW 

The Einstein field equations are solved to first order in h, and it is 
determined that the time-dependence exp(-ikt) cannot diverge because imag- 
inary frequencies k are disallowed, thereby demonstrating the stability 
- or so maintain Regge & Wheeler - of the Schwarzschild solution. Their 
technique is contested by Peres & Rosen (1959) as unsatisfactory because 
it does not take int^ account properly the nonlinear effects of gravita- 
tional radiation. Indeed, Peres & Rosen maintain that small oscillations 
cannot be stable for any gravitational field assumed to be asymptotically 
flat at infinity - but this criticism has been either ignored or forgotten. 
The work of Regge & Wheeler is pursued further by Manasse (1963), Brill & 
Hartle (1964), Doroshkevich, et aj[ (1965), Peters (1966), Vishveshwara (1968, 
1970), Vishveshwara & Edelstein (1970), and Zerilli (1969, 1970). 


Regge & Wheeler note that Einstein's field equations (3-1) for the exter- 
ior Schwarzschild background • are simply 



(i-ssa) 
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and they argue that a perturbation of the metric as in (3-6) likewise re- 
sults in a perturbation of the field equations which still equal zero 

= Fy<|) ♦SR A> 0»hO. (3-5SW 


By virtue of (3-54), then (3-55) reduces to 


R/v = $ R,» • o , 


(J-S u) 


which is given in Appendix A, Eq (A-6). This perturbation can likewise 
be represented by a relation due to Eisenhart (1926) 






noting that is a tenser although r p is net. Eq (3-56) is the covariant 


uv 


generalization in curved space of the Schrodinger equation for a massless 
Spin-2 particle in flat space (Regge, 1957). 

The most general form of h consists of an odd-parity (-1)*" + ^ case 


h - 

J" 


o o 
o o 

0 o 

* * & 


* 




as well as an even-Daritv (-1)“" case 


(3-S7a) 




V 

Jt 


H.M \<t fib 


n$» 


where the asterisk ★ means h..=h... By performing a gauge transformation 

I J J < 

- known as the Regge-wheeler gauge (Appendix I) - such as is discussed in 


(3-57V4 
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Appendix (B-5), the perturbations h of (3-57) can be reduced to the fol- 

'j 

lowing canonical form: 

Qdd -Pari ty Pe rturb a tions 

0 o O K,(r) 

o o o 

o o o o 

« * ° ° J (3-SM 

Even - Pari ty Per turbations 

V M. 0 0 1 

* l/o O 
0 O »*k O 

O O 0 » a KsJeJ (**&,) 

where aximuthal symmetry has been assumed (M=0). This particular form of 
h is what is meant by the "Regge- Wheeler gauge." It consists of two odd- 
parity (lower case) radial functions h Q and h-| as well as four even-parity 
(upper case) radial functions H Q , H 1 , H 2> and K. 

The Regge- Wheeler canonical perturbations of the Schwarzschild metric, 
then, substituting (3-58) into (3-6) and (3-53), amount to 

Odd-Parity Schwarzschild Metri c 

<U l = UV ** ir 1 ♦ r*(«* * 

+ 2 (h^2t+h,WJir)^ a **** (*ui® ^)Pfc»sft) 

E ven - Parity Schwarzschild Metric 

Js l = (ltH 0 P L a' k Vit l + (» + H t P L A 4i(t )jtJir 2 ^ 

+r 1 (l+KP u A lkt X^»«»W) + 2Hf|* 4tir, 






37 


The problem remaining is to deletin' n 2 tne odd-parity and even-parity 
radial functions. These are obtained by solving the differential equations 
(Appendix H) which result from tne substitution of h in (3-58) back into 
the perturbed Ricci tensor and field equations [2-ZC). For odH-narity 
case Regge & Wheeler obtain the following "Schrodinger" type equation 

9 . t(k'-yjQ = 0 <H> 

where Q = eh^/r and r is given by 
The effective potential V e ^ is 



For the odd-parity quadrupole oscillations (L=2) the effective potential 
V eff in (H-8) can be shown to be that in Figure 5. 



Vishveshwara (1970) treats the scattering of gravitational radiation from 
such a Schrodinger- type effective potential. 


The even-parity radial equations (Appendix, H-10) prove to be more dif- 
ficult and are not reduced to a single second-order differential equation 
until the work rf Zerilli (1970b), which is presented in detail in Appendix J. 
He obtains a homogeneous "Schrodinger" equation (J-7) for the radial function 
) from which h^ in (3-58b) can be determined: 

Km Til \i . 1 A _ a (j”7) 



[r-'Hv :0 
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where r is given by (ri-5) and where 




2X*(X-H)r |3 + IsXmVfjBv? * 






V^(r) in (J-8) is the even-parity dual of (H-8). Eqs 'J-7) and (J-8) 
are sometimes referred to as Zerilli's Equations. 


Ter .or Harmoni c Decomposi tion Wi th Sources 

The original treatment due to Regge & Wheeler (1957) is the first decom- 
position into tensor harmonics of perturbations on a spherically symmetric 
background. Tensor harmonic decompositions are useful because they provide 
for a separation of variables in the differential equations, and for this 
reason are further' analyzed by Mathews (1962), Peters £ Mathews (1963), 

Thorne, et al_ (1967-1970), and Zerilli (1970a, c). Such decompositions are 
also important in generalizing the linear perturbation technique used on 
the Schwarzschild background to include sources. 

Recall that Regge & Wheeler do not account for sources. Theirs is a 

canonical nonspherical perturbation h of the Ricci tensor such that R^ - 0, 

pv pv 

as in (3-56). They thereby find Green's functions from which solutions of 
inhomogeneous equations with sources can be constructed. On the other hand, 
in order to establish realistic perturbations with astrophys - cal applications 
one can account for sources with a non-zero stress-energy tensor T . For 

IJV 

example, this treatment of the problem is reasonable in the study of the per- 
turbation due to a small particle falling into the Schwarzschild geometry 
from infinity. 


Such an analysis is formulated by Peters (1966) who considers first-order 

linear perturbations of the Einstein tensor rather then 6R , and 

pv pv 

accounts for sources of small mass m (m<<M s ) by means of a variation in the 
stress-energy tensor T . It is a weak-field approximation (3-6) and must 
be classified as a slow-motion one also, to the extent that the dilatation 
of mass (3-5) retains the condi tion m(particle)<<M s (Schwarzschild) . Peters' 
linear perturbation of the Einstein equations (3-1) results in the following 
field equations for a background n of constant curvature (Ricci flat with 
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In (3-60), f represents an arbitrary gauge vector introduced in order to 
reduce these field equations to a canonical form, once the unperturbed 
background metric n is chosen. Furthermore, in order that (3-60) be con- 
sistent, 6T must obey a conservation law. Taking the variation of (3-38) 
shows that the representation of the perturbing source i must satisfy 



Peters represents the part of ■iT due to the perturbing particle as 

sT r H - m 

rm*. . («*• 

z y (s) is the space-time position of the particle in the metric. Peters 
finishes his treatment with a Green's function method for curved space and 
specializes his results to the Schwarzschild background. 

Peters' field equations (3-60) and conservation law (3-61) assume the 
stability of the background metric n ; this work is not a stability anal- 
ysis as is that of Regge & Wheeler. Furthermore, the perturbation is of 
first-order and thereby neglects radiation reaction. The perturbing par- 
ticle is assumed to follow a geodesic in the unperturbed background for 
the computation of 5T in (3-62) - an assumption which is not correct. The 
perturbing particle follows a geodesic in the perturbed metric g and not 
a geodesic in n , according to Einstein's theory, (-ee Figure 20) 

UV 

Nevertheless, Peters' representation of Schwarzschild perturbations is a 
very useful approximation, as has been demonstrated by Zerilli (1970a, c), 
who decomposes Peters' field equations (3-60) into tensor harmonics using 
the Regge-Wheeler gauge. 

Zerilli (1970c) develops an orthonormal set of tensor harmonics (Appendix K) 
for application to gravitational radiation theory, based upon the earlier 
work of Regge & Wheeler (1957) and Mathews (1962). Because Peters' field 
equations (3-60) are of the form 

QlVl*-” 5 T» . <*» 

* Peters and Zerilli assume h,„ is small if T, is small. 

U v 
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where Q is a rotationally invariant operator, Zerilli expands both sides of 
(3-63) in tensor harmonics in order to separate the angular variables 0,$ 
from the radial equations. Choosing the Regge-Wheeler gauge (Appendix I) 
the h decompose as 

u V 




where tensor subscripts "uv" have been replaced by the double stroke in 
order to simplify tensor notation. The superscript (e) corresponds to the 
electric or even-parity case while (m) stands for magnetic or odd-parity. 
Referring to the tensor harmonics in Appendix K, one can verify that (3-64) 
is essentially the same as the Regge-Wheeler gauge in (3-58). 


(3-M) 


Likewise the stress-energy tensor for the perturbing particle ST can be 
expanded into the tensor harmonics of Appendix K: 



The respective coefficients A^,. . . ,F. M are defined in Appendix L and repre- 
sent the amplitudes of the various tensor harmonic components of the stress- 
energy tensor which "drive" the perturbation of the metric. 


Substitution of the tensor decompositions (3-64) and (3-65) back into 
Peters' field equations (3-63) or (3-60) results in one equation for the 
magnetic (odd or L+l ) parity and one for the electric (even or L) parity 
case with tensor harmonics on both sides. The coefficients of the tensor 
harmonics on each side of these equations muct then be set equal, which 
gives the radial equations for the Peters-Zerilli analysis in the Regge- 
Wheeler gauge with sources (Appendix 0). The sources, in turn, must satisfy 
the conservation conditions specified in (3-61) and Appendix M. 

Zerilli next takes the Fourier transform of the radial equations (Appen- 
dix P) and then reduces them (P-1 and P-2) to second-order "Schrodinger" 
equations with sources, in the same fashion as is performed with the source- 
free Regge-Wheeler problem discussed earlier. Introducing two auxiliary 
radial functions and R^ he obtains Eqs (3-66) & (3-67): 
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Magnetic-Parity ' Sc hrodinger 11 Eq uation Wi th Sources 


.2n (m ) »- On) -1 Am) (jn) 

(s-tta) 

where 





(j-UO 


+ Om} . 


Electric - Parity " Schrodinger" Equation With Sources 



where 




1 (j-UW) 

V 0 »*(Xr+3»V i 





X * i(L-»XL«) . 


The original radial functions h Q , h-j , Hq, , H^, and K can now be derived 
from in (3-66) and (3-67). These relations as well as those for the 
auxiliary source coefficients C-j , C 2 , and B are all given in Appendix Q. 
Equations (3-66) and (3-67) are also referred to as "Zeriili's Equations." 
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Upon evaluation of the coefficients of Appendix L for a Schwarzschild 
geodesic, Eq (L-8), one finds that for a particle falling radially (say 
along the z-axis) into the Schwarzschild geometry, all magnetic (odd or 


L+l) parity contributions and all M f 0 electric perturbations vanish. For 
this case, the source function S^q in (3-67) reduces to 


riW _ 

“ X**3W 


w!-?)[ic • =gy“ 1> ’ m 


where 




(3-C8W) 



These equations are of current astrophysical interest, in that several 
groups are evaluating them for higher L by numerical integration techniques, 
as we shall see in Chapter 5. 


Zerilli solves his equations in h for L=0 and L=1 in the above case of 

pv 

a particle in orbit about a Schwarzschild geodesic. The magnetic monopole 
(L=0) is identically zero. The electric monopole (L=0) case correspond to 
a mass perturbation which simply adds the particle mass m to the Schwarzschild 
mass M (y Q is a constant), depending upon the position r of the observer: 


vfr.tt 

Jl = 


J I- 23 , r<R(t) 

| , r>R(t) . 


(V&a) 


The magnetic dipole (L=l) contribution is the angular momentum of the 
equatorial orbit assumed, which is 

l_, = mo. 1 ,, 

a«R*sin©$ J , 

while the electric dipole (L=l) represents a coordinate transformation of 


g by h^^, as we know it should - but not contributing dipole radiation due 
to conservation of momentum. For L>1, analytic solutions have not been obtained. 
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The Bondi News Function & Multipole Expansions 

An entirely different approximation technique is the method of mult’oole 
expansions. Multipole approximations have historically provided a great 
deal of insight into the study of radiation theory, certainly in the clas- 
sical vector theory of electromagnetic and gravitational radiation, as dis- 
cussed earlier. Likewise, multi pole approximations have found their place 
in the tensor theory of gravitational radiation. 

In order to illustrate the multipole approximation in General Relativity, 
it is best to consider the work of Bondi (1960, 1962, 1965) who has been 
expecially concerned about the meaning and the physical existence of grav- 
itational radiation from gravitationally bound and otherwise isolated sys- 
tems. On the basis of causality arguments, Bondi (1962) considers only re- 
tarded solutions and he places great emphasis upon the requirement that 
gravitational radiation must result in a loss of energy and hence mass of 
source in order to be consistent with the spirit of General Relativity. 
This means that for gravitational radiation from a mass quadrupole, there 
must be a secular change in the mass monopole. Thus the monopole must be 
coupled to the radiative multipoles, such as the quadrupole, of the multi- 
pole approximation. 

Bondi's development (Appendix R) treats the stationary-radiative-stationary 
transition of an outgoing gravitational "sandwich" wave (Figure 6) as it 
passes through some retarded hypersurface "u=t-r/c=constant" in an asymp- 
totically flat region of space: 

STATIONARY - — RADIATIVE — STATIONARY 



EI5UBEE 
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In Figure 6, regions 1 and 3 are stationary and Minkowski flat with no 
curvature. Region 2, however, is radiative and empty, behaving in much the 
same fashion as the Regge-Wheeler treatment presented earlier. Such a trans- 
ition is assumed in order for the Bondi formalism to be consistent with 
Huygens' principle. It is tantamount to the statement that r : -avitational 
waves in a stationary-radiative-stationary transition have no "tails": 

Figure 7 

-fWlr A/ VWlr 

® ! ® !© ® ! ® ! © 

! i ! 1 

(A) Bondi s(Huygens)As8umprion (B) Gravitational "Tails' 1 


The existence of "tails" (Figure 7b) means that the background geometry 
continues to oscillate as a consequence of the "shock" perturbations induced 
by the passage of the gravitational wavefront. It has since been demonstra- 
ted, however, that Bondi's assumption of no "tails" (Figure 7a) is indeed 
invalid. Papapetrou (1969) and Hallidy & Jam's (1970) have investigated 
the problem of the existence of a final stationary state (Region 3 of Fig- 
ures 6 & 7a) and they have concluded that stationary-radiative-stationary 
transitions cannot occur for axially symmetric (Bondi) radiation of a finite 
multipole expansion. Nevertheless, it is informative to consider Bondi's 
representation because Sachs (1962) has demonstrated that the condition of 
axial symmetry can be relaxed, and the "tails" eliminated. 


The critical and controversial question of outgoing radiation at infinity 
is addressed by Bondi by choosing the axially symmetric (independent of <t>) 
metri c 



Axially Symmetric Metric 


AJU* + 2Bl«4r 


(3-TO) 


which is extremely well suited to the problem of the Sormerfeld radiation con- 
dition. The Bondi coordinates are X M =(u=t-r,r,u,0* Asymptotically the Bondi 


metric reduces to its Minkowski form 
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For sufficiently large values of r, Bondi (Appendix R) assumes for the 
arbitrary functions in his metric (3-70) or (R-l) the following: 

X = f + ©(*) 

F = f •*■ 0(W 
B = [1+ £] + ©(*») 

A = [,-“]+©(« 

He shows that f=b=0. Then he derives from the field equations R =0 

MV 

in Region 2 of Figure 6 a differential equation for the time rate of change 
of the function M(u,o) in (3-71) 



Based upon correspondence arguments with the static case, the "mass aspect" 
M(u,Q) is shown to be related to the mass of the source m(u) by 

r* 

m(u}= -H M<n,*l *>"•<*« . (3-u) 

Differentiation of m(u,Q) in (3-73) with respect to u and substitution of 
M q from (3-72) gives Bondi's result that 

te s *in^*'' new • 

0 

,di calls the function c=c(u,Q) a "news function" because from it alone 
can be determined the entire behaviour of his metric. In particular, from 
relation (3-74) one can see that knowledge of c(u,°) on some hypersurfaci. 
u=constant would determine the secular change of the mass monopole source. 
Furthermore, if the news function is not zero, then from (3-74) there must 
be a monotonic decrease of the mass m('i) of the source. Bondi's method is 
effectively a method of asymptotic linearization, and is illustrated in 
Figure 8. 

In a paper due to Penrose A Newman (1965) a set of ten geometrical quan- 
tities, defined for asymptotically flat space-times, are shown to be conser- 
ved during the emission of gravitational radiation. One of these quantities 
is particularly relevant to the Bondi representation arid states that a 
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Figure 8 



4 ? 


transition between two stationary states (Figures 6 & 7a) is allowed if 

D* - mQ = Conserved } (VIS) 

where D is the dipole, m is the mass, and Q is the quadrupole moment of 
tne source. (3-75) is net sufficient tc insure that the fir, a’ state is 
stationary, however. 

Bondi's work is generalized by Sachs (1962) who relaxes the condition 
of spdvial symmetry on the fields, allowing for tv/o polarizations instead 
of Bondi's one. Sach's generalization is Used upon a complex r.ews function 
c, from which he determines the rate of loss of mass due to gravitational 
radiation (Figure 8) as 

m - -<CC*> > (3r74) 

which is identical to Bondi's (3-74), except that carets arc used to rep- 
resent the integral or average over the sphere at infinity. Morge- & Peres 
(1963) identify the Sachs news function with the mo-.s quadrupoTe tensor 
and then from (3-76) derive the relation 

m =-h V* > (377) 

in close agreement with the results of Einstein's linearized ouadrupole 
radiation (3-29) and of the EIH approximation (3-41). Sach's paper (1962) 
has gained more significance with the work of Pa ( jdp^trou (1969) and Hallidy 
& Janis (1970) mentio: ad earlier because the problem of gravitationrl "tails" 
in the Bondi representation does rot exist if the condition of axial sym- 
metry is removed. 

Early criticism of Bondi's approximation method is also presented by 
Bonnor (1963) and Bonnor & Rotenberg (1966) who make use ' i nly of Bondi's 
metric (3-70) & (R- 1 ) . Bonnor (1963) argues that only part of the news func- 
tion ctu/J), the linear part, is known to the observer i n Bondi's *pproxima 
tion method. The nonlinear and critical part is indeterminable. If you 
plug only the linear part of c(u, ) for a quadrupole oscillator, as an 
example, into Bondi's formalism (R-7) the oscillator gees ,nst. . e. Bonnor 
& Rotenberg (19b6) eliminate this shortcoming of the Bondi method by using 
instead Bonnor's (1959) "double series" approximation (3-9). Up„n passing 
to nonlinear approximations, they demonstrate that under forced oscillations 
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the source loses maos due to quadrupole-quadrupole interaction and exper- 
iences radiation recoil under quodrupole-octupole interaction. They also 
point out that the axial symmetry of the Bondi method produces "tails. ‘ 

Bonnor & Rotenberg conclude, furthermore, with the statement that the ex- 
istence of radiation from a freely gravitating system is still an outstand- 
ing problem (1965). 

The error in Bondi's original approximation and the basis of the origin 
of the "tails" is a naive one. Using the multipole approximation methc* in 
the linear vector theory of electromagnetism, one can determine the solu- 
tions for a stationary-radiative-stationary transition (Figure 6) K y satis- 
fying a finite number of linear conditions. Bondi errs, however, in extend- 
ing this result to a nonlinear tensor theory of radiation. In the general 
relativistic tensor problem, Bondi needs an infinite number of nonlinear 
conditions to make his approximation work for the axially symmetric case. 

Nevertheless, the Bondi formalism and its association with the notion of 
gravitational news or data on a radiative hvpersurface has inspired a great 
deal of research into methods of asymptotic linearization - in particular, 
the tetrad Tormalisms of the next chapter. 

The High Frequency Approximation 

Relevant to any discussion of approximations in General Relativity is the 
high-frequency approximation of Isaacson (1968) who arrives at an effective 
stress tensor by analogy with the geometrical optics limit of electromagne- 
tism. In the high-frequency limit the gravitational fields uncouple from 
their sources and attain an existence all their own. The Isaacson stress- 
tensor is gauge-invariant, second-rank, symmetric, and (most important) 
not a pseudotensor (like the Landau-Lifshi tz pseudotensor of Appendix C). 
Treatments using the I .aacson approximation include Price & Thorne (1969, II), 
Ipser (1971), and Thorne (1968). 

The Isaacson effective stress tensor linearizes the radiation zone, placing 
the burden of the nonlinearities upon the region of tie sources - as does 
the Bondi formalism. 



CHAPTER IV 


EXACT SOLUTIONS & ASYMPTOTIC APPROXIMATION 


Looking for exact solution's of Einstein's nonlinear field equations is 
very much like taking a peak inside Pandora’s box. In fact, one m i g ht sim- 
ply make this statement and leave the whole thing at that. Nevertheless, 
any treatise un gravitational radiation theory must assess the physical 
meaning of radiative solutions, in particular those which have been found. 


Indeed, it is precisely the physical interpretation of the exact solutions 
which makes progress in the nonlinear theory of radiation extremely awkward. 
All of the radiative solutions either contain naked singularities, or for 
those singularities which can be identified with a source, there always 
appear singular 2- surfaces which transport energy and/or momentum from in- 
finity into the source, where it can then be radiated outward again. Con- 
sequently, the most profound theoretical question about nonlinear gravita- 
tion? 1 radiation - that of energy transfer - remains unanswered. And until 
it is answered, the very existence of nonlinear tensor radiation will 
remain in doubt. 


Exact Plane Wave Solutions 

Plane wave solutions of Einstein's empty space field equations G = 0 
or R =0 are investigated by Brinkmann (1925) with the metric 

pv 

is* = iu{2iV JixSiY*] (4-0 

although he does not characterize the solutions as radiation. They are more 
formally addressed by Rosen (1937) whose metric is derivable from (4-1) 

is* - 2Liniv - (f , Jx 1 +G4y) (4-*A 


or as stated by Bondi (1957) 
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is* = , C^-*W) 


f % 9 

where u = 1 - 5 , e=fc(u), a=a(u), and a u = 5 u a= uSy. Rosen concludes, however 
that the solutions of ( 4 - 2 ) cannot exist because they possess physical 
singularities. That this conclusion is too severe is pointed out by Bondi 
(1957) and Bondi, Pirani, & Robinson (1959) who maintain that Rosen does not 
distinguish adequately between coordinate singularities and physical singu- 
larities. Rosen's metric is empty, and it is flat ( R aBY( 5 = 0) if 


+ . («-» 


As pointed out by Bondi (1957), the coordinate transformation 


**T-S«ct-x 






is nonsinguiar for u>0 and (4-2b) becomes 


is* * C*a* - (AM,W) - /C (5-s) 

which represents a non-flat region between two flat ones provided e. t 0 in 

* u 

the last term. This is a "plane-wave" zone of finite extent whose amplitude 
is determined by e=e(u). A more general plane-wave, with a variable plane 
of polarization is given by Bondi (1957) as 

is* = iVtf- i J*) - (ct-J)*[cwh ifi (^\ it) + 

Smh2jU«20(A| t -JS^-2»Mi2psm2eJh|lt| | ^ 

where Q=s(u). The empty-space condition for (4-5) is 

(et-1) [|£+ 8* s»h* ij £\ . M 


A discussion of singularities is presented by Bonnor (1957), while Kundt (1961) 
investigates the general nature of "plane- fronted" gravitational waves. 

Recent work by Szekeres (1970, 1972) and Penrose & Kahn (1971) has made 
significant progress in the study cf exact solutions for colliding plane 


★ 

They are not really "plane" because the departure from flatness R^ /Jt 
depends upon y and 2 , fcr propagation along x. 
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waves (essentially a sort of "graviton-graviton" interaction). The nonlin- 
earities are taken fully into account, which is of the utmost importance, 
because the very notion of '.inear superposition (Maxwell's theory) is at 
stake. Indeed, as we have seen in Figure 8, the Bondl-Sachs asymptotic 
representation essentially linearizes the field equations in the radiation 
zone by virtue of the fact that there is sjperposition of therews functions 
(c = Cj + Cg) - which places the burden of the nonlinearity upon the near 
zone In the region of the sojrces. Szekeres (1970, 1965) maintains that 
superposition Is simply invelld for colliding "sandwich waves," in conflict 
with the Bondi -Sachs formalism. Penrose & Kahn employ impulsive fi-f unction 
plane-waves, arriving at similar results as Szekeres. In either case singu- 
larities exist, which make the physical interpretation difficult to assess; 
perhaps the singularities disappear for more realistic, curved wave-fronts. 

Exact Cylindrical Wave Solutions 

Exact cylindrical waves for R = 0 are presented by Einstein & Rosen (1937) 
and Rosen (1937, 1954, 1956, 1958), based upon the static axi ally-symmetric 
metric of Weyl (1918) and Levl-Civita (1919), 

is = - ***Uf*U) (+«* 

which reduces to the Laplace equation (in culindrical polar coordinates) 

?**\+y*f * o (4-8« 

Vvy*’* . (4-ad 

Rosen's cylindrical metric is obtained by merely Interchanging the roles of 
z and t In (4-8> , whereby 

is* - *f)- -i*V . (4-9* 

Substituting in R^ v =0, the Rosen metric (4-9a) gives 


VfV **« s0 

(4-9W 


(4-3e) 

v*pV t • 

(4-9* 
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For the Sommerfeld radiation condition, the solutions of (4-9) for 
k = mp/c are (Rosen (1954)) 


i * itfkftoii .« ♦ WM4 ♦ ktfw l U-|g) 

♦UWtyW -MMfl «««* [ 

The last term in (4-10) is a secular variation or the metric which can be 
interpreted as a loss of energy due to radiation. Rosen (1956) objects to 
this inte» tation, claiming that there is no transport of energy and 
momentum because the pseudotensor vanishes. Weber & Wheeler (1957) demon- 
strate that it is nonzero. Harder (1958, 1961, 1969, 1972) studies exact 
cylindrical waves also. 


Although the singularity along the axis of symmetry can be identified with 
a source such as a thing rod, which makes it more tenable, this rod is 
infinite in length. Thus the nice cylindrical solutions which radiate en- 
ergy outwards require an infinite source. Furthermore, there is reason to 
believe that energy is transported back around and down the axis of symmetry 
along the source singularity with zero net transfer of energy. 


Exact Spherical Wave Solutions 

Birkhoff (1927) proves that any general spherically symmetric solution of 
the empty space field equations can always be transformed into the static 
Schwarzschild metric. By virtue of Birkhoff's theorem, no spherically sym- 
metric gravitational radiation can exist, an important result for astrophys- 
ical theory. Spherical gravitational waves in the linearized theory have been 
investigated by Boardman & Bergmann (1959) and Bonnor (1959). As stated by 
Weber & Wheeler (1957), however, spherical gravitational waves can never be 
truly spherically symmetric, just as is the case with electromagnetic ones 
(the fixed point theorem of topology). 

Robinson & Trautman (1960, 1961), nevertheless, have discovered a class of 
exact solutions which correspond to a form of expanding radiation which they 
classify as spherical gravitational waves. 

Other Sol ut ions 

Additional solutions can be found in the work of Peres (1959) and Takeno 
(1956, 1957, 1958). Takeno treats the problem of plane, cylindrical, and 
spherical waves in the nonsymmetrical unified field theory. 
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Invariant Formulatio n Of Gravitational Radiation Theory 

As we have seen in some detail, coordinate or gauge transformations can be 
used to create and destroy apparent gravitational radiation. But that radia- 
tion which cannot be transformed away is cons’ iered physical. The Riemann 
curvature tensor R „ . likewise cannot be transformed away by means of a 
change of coordinates. As a consequence, theoreticians - enchanted by the 
prospect of relating the Riemann curvature tensor to the existence of 
physical radiation - have begun to classify its radiative characteristics. 

Recalling that the transversal ity of E-H = 0 and ! E j = | H j i n the radiation 
zone of electromagnetism can be expressed in terms of the null vectors k u 
(defined by k a k a = 0) as 

k“F = 0 k a F* = 0 , 

aS aB 

it is argued that the Riemann tensor possesses similar syirmetry properties. 
Based upon a theorem due to Debever (1958) and Sachs (1961), a similar set of 
null vectors k a satisfy the following equations (see Bcnnor (1963) for a 
discussion) : 


Petrov Metric Type Equations For Nul 1 Vectors No . Distinct Rays 


I 

II 
D 

III 

N 


k ra R 3]ij[y R 5 1 k kJ = 0 
\*[« k e] k V = ° 
R aBy[6 k E l kV “ ° 

R o6y [s m E ] n 

'Ws'^e] = ° 

R ‘ k J = 0 

aByo 


r m im a m Y = 0 


4 

3 

2 

2 

1 


where the Petrov (1954) types (Appendix U) have been applied by Pi rani (1958) 
as a means of classifying the radiative properties of the Riemann curvature 
tensor. This can be stated in terms of a Penrose (1960) diagram: 


* X Penrose Diagram 

D— H 

0— 



54 


which represents the substructure of the Petrov metrics. There are three 
spaces (three vertical columns) which may or may not be composed of a sub- 
space. The left-hand column (reading down) is composed of "Petrov-type I 
nondegenerate," “Petrov-type I degenerate" (D), and flat-space (0). The 
middle column is "Petrov- type II," and "Petrov- type II null" (N). All the 
metrics are "algebraically special" except type I (the most general type). 


The Riemann tensor has the form 


or dropping the subscripts 


(Ml») 


R 



(Mfc) 


Hence, Petrov-type N is radiative (it goes as r' 1 ). The same formalism also 
applies to the Weyl tensor. In terms of the Bondi news function, the Petrov 
types in Eq ( 4-1 lb) are (Trautman (1962)) 

Appendix U contains details on all of the exact solutions discussed earlier 
and their associated Petrov-Pirani classification. 

The Petrov classification of fields using the Riemann tensor, then, is 
introduced by Pirani (1957) into gravitational radiation theory with the 
hope of constructing an invariant basis for establishing the existence of 
gravitational radiation. Referring to Appendix U , one can see that this 
endeavor is no* entirely successful. 


Invariant Formulation & Asymptoti c Approximations 

The embodiment of the Pirani invariant formulation of gravitational radia- 
tion using the Petrov classification can also be designated as a tetrad cal- 
culus. Witten( 1959) and Penrose (1960) likewise develop a two-component 
spinor formalism. The tetrad and spinor calculus are then combined together 
by Newman & Penrose (1962) in a formalism with emphasis upon an expansion of 
the Weyl tensor rather than an expansion of the metric tensor in inverse 
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powers of r, due to Bondi at a! (1962). Hopefully, the Newman- Penrose 
method will offer a reconciliation of the invariant formulation with the 
asymptotic approximation of Bondi. 


Exact Solution For A Gravitationally Bound Sys 

An exact radiative solution from a system representing a bound self- 
gravitating state is now available, and is of great importance in that it has 
physically meaningful singularities associated only with the material sources 
and constitutes the first exact solution of any kind for a system in motion. 
It is the work of Bonnor 4 Swaminarayan (1964), Swaminarayan (1966), and 
BiSak (1968, 1971). It is an analogue of the classical Born (1909) problem 
of the uniform acceleration of two charges (see also Rohrlich (1965), 118). 

The solution is "Petrov-type I, nondegenerate." It comes from the 
metric (Bicak, 1971) 

“ here 

i * ^ jfV-t 1 ) t (Jiy-Oft-f-lii) ' 

and where a-j, h ( >0, I^O, and k>0 are arbitrary constants. The solution 
describes two pairs of mass points which are represented by the world lines 
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and are thus uniformaly accelerated by ± » / (2h i ) in a Minkowski background 
There are four possible cases for particle mass: 




,a 


2 X >0 

a ‘ 

a*«o 


} 


v * i ( kw ^ 

v * i (VMO 

( 4 - MA ) 


a, arbitrary 
a,>o 

a, *0 a»>o 

Eq (4-14a) corresponds to two mass dipoles (positive and negative mass); 

Eq (4- 14b) corresponds to four masses of positive mass with stresses; 

Eq (4-14c) corresponds to a pair with negative mass moving under stress, while 
the other two particles move freely; and Eq ( 4-1 4d) represents only two 
particles with the same positive mass, connected by stress. 


Transforming to spherical coordinates, the line element (4-12) becomes very 
involved, but Bicak gets for the Bondi news function 




(iHSs ft 

where u = t - r and 

*=-1=1 - i*L « 

U. tt, v “ m > 



U; » J u* ♦ 2H; 

t»isa 


ftriaaii 




[ tt.TIx ’J 

nrwM iMt. ww. 


The Bondi news function (4-15) begins with an asymptotic expansion contain- 
ing r* 1 , irrespective of the choice of a.. , h^ , k (that is, for all cases 4-14) 
This is to say, the asymptotic expansion of the Riemann tensor tetrad is rad- 
iative according to Bondi's method. Bi<*ak (1971) demonstrates that the radia- 
tion pattern, furthermore, for the Born solution of electrodynamics is 

while that for the Bonnor-Swami narayan solution of gravitational radiation is 

( Jp / to) 6w = t»*e . d-iuj 

BicSk, in addition, points out that the particle masses in the Bonnor- 
Swami narayan solution are the same as those occurring in Bondi’s (1957) paper 
on negative mass in general relativity. There is strong evidence that the 
particles in the Bonnor-Swami narayan solution are gravitational monopoles. 



CHAPTER V 


GRAVITATIONAL RADIATION IN ASTROPHYSICS 


Let us turn now to astrophysics. A wealth of astrophysical applications 
exist? for the linearized solutions of the tensor theory of gravitational 
radiation - even though the exact, radiative solutions which are known hove 
been difficult to interpret. Such applications must be astrophysical in 
nature simply because the order of magnitude of the radiation is so small 
(G/c 1 ' = 2.7xl0' 60 ). For quadrupole radiation alone, 

« tl *»“ <5^ 

Until exact solutions which are physically meaningful are established, it is 
nevertheless informative to adopt the linearized theory, to forget any contro 
versy which may exist in regard to radiation from the Kepler problem, and to 
consider the observable phenomena which might result. 

Gravitational Radiation From A Spinning Jacobi Ellipsoid 

Consider a spinning Jacobi ellipsoid (triaxial inertias 1^, ^ 2 ’ *33 anc * 

semi-axes a, b, c) whose body axes are chosen to coincide with the principal 

2 

axes and whose angular velocity Q is along the X-axis. The Einstein- 
Eddirgton quadrupole radiation formula (3-30) becomes in this case 

iVa* 

-fe * I ►•» 

where I = I ^+1 and e = (*11**22^ * (a-b)//(ab) is the equatorial ellip- 
ticlty. For an ellipticity of unity (e=l, I22 = Q) formula reduces to the 
original quadrupole radiation (3-30) for a spinning rod. On the other hand, 
if I-j -|=I 2 2 th® Jacobi ellipsoid degenerates into the Maclaurin spheroid and 
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ceases to generate gravitational radiation because it has zero ellipticity 
(e=0 in 5-1) . 


The radiation reaction torque acting an the ellipsoid due to gravitational 
radiation would appear to contribute a loss in angular momentum L by an 
amount 


4k- 

At ~ £ At 


3 Z 

s 




(S-2^ 


For the triaxial ellipsoid, this is (Chandrasekhar, 1970) 


<L 

At 






Naively treated, Eq (5-2) should result in a spinning down of the ellipsoid. 
Consequently, the exponent n of n n in -d£/dt is known as the "slowing expon- 
ent" and has been used to indicate the presence of multi pole radiation 
mechanisms, 



The Crab pulsar has been observed to have (e.g. Ruffini & Wheeler, 1969) 


n = 5.76 + 0.65 (Optical) 

n = 3.6 + 0.6 (Radio) 

which would imply the existence of a quadrupole radiation mechanism (n=6) 
such as (5-1). Assuming a mass of 0.786 Mq, a mean radius of 9.75km, and 
eccentricity of 8xl0 -4 , and a period of 33 msec, then a "neutron" star spin- 
ning as a triaxial ellipsoid in (5-1) can account for the observed rate of 
change of the Crab pulsar's period of 4x10 sec/sec. The power radiated 
would be 

-d£/dt = 2 x 10 38 erg/sec. 

However, energy dissipation and angular velocity cannot be related in this 

fashion. Chandrasekhar (1969, 1970) has pursued the evolution of the Jacobi 

2 5 

ellipsoid under gravitational radiation in the context of the "post ' - 
Newtonian" (E1H) approximation. He has arrived at the surprising result that 
a triaxial ellipsoid increases in angular velocity as it loses angular momen- 
tum and radiates energy (5-1). It asymptotically approaches the Maclaurin 
spheroid, whereupon it ceases to radiate, at the point of bifurcation (see 
Figure 9). The Maclaurin spheroid is then dynamically unstable under 
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gravitational radiation reaction, exhibiting the possibility of fragmen 
tation. 





(1 - b/ A) = 


FIGURE 9 


The dissipation of energy, then, under gravitational radiation can be 
derived from potential and internal energy, and reed not be at the expense of 
rotational kinetic energy. Chandrasekhar's results have significance for 
gravitational collapse and the formation of "black holes" in astrophysics, 
because asymmetries during collapse should be radiated away as gravitational 
radiation if his approximations are applicable. The collapsing object 
should evolve into a spheroidal (Maclaurin), nonradiating "black hole" - 
which is secularly and dynamically unstable. It could fngment and even 
bifurcate under this EIH analysis, in contrast to the results of Penrose 
(1972) and Hawking (1971) . 

Figure 9 also explains the "glitches" or sharp drops in the rotational 
periods of pulsars. Starquakes briefly create a Jacobi ellipsoid which 
quickly becomes Maclaurin. 

Gravitational Radiation From The Kepler Problem 

Another system of astrophyslcal interest is tie iwo-body Kepler problem which 
has already been discussed at some length in the EIH approximation. It is 
particularly interesting because it is a case of a self-gravitating system 
not undergoing forced oscillations but yet experiencing a time-varying 
quadrupole moment. 


*"Black hole" is the acronym attributed to the cpometry in (3-53) when the 
radius of the mass M is less than the Schwarzschi la radius 2GM/c 2 . 



2 ^ 

Recalling for a circular orbit that =G (m-j+iT^) /r' ana that the elllpticity 
is unity (e=l), then (Landau & Lifshitz, 1962) Eq (5-1) becomes 


it _ 31 G* 
~ it = S ? 


This formula should represent the grav ional mass quadrupole radiation 
from a stellar binary, provided such radiation exists and boundary conditions 
allow for it. Noting that the total energy ofaKeplerian conic is 



then substitution of (5-3) into r of (5-4) gives 


(A G* m l m 1 (m,tmQ 

S c s f 5 * 





which represents the rate of decay of the binary - but only approximately 
because the circular orbit assumption must be treated adiabatically as the 
conic starts to collapse. Eq (5-5) is often used by astrophysicists when 
considering the lifetimes or relaxation '.imes of b’ ry systems undergoing 
gravitational radiation. Such aiguments, for exampte, are used to dispel 
binaries as possible pulsar mechanisms because they collapse too rapidly. 

This interpretation of gravitational radiation reaction (binary collapse) 
effectively ignores or considers invalid the results of Hu (1947), Peres (1959), 
Havas & Goldberg (1962), and Havas & Smith (1965) who all find a dilatation of 
the binary system under gravitational radiation - which is consistent with 
the behaviour of observed pulsar periods. Thinking of the binary as an 
ellipsoid, however, it should increase in angular velocity (collapse) under 
Chandrasekhar's interpretation vf the EIH method. 

The first adequate treatment of the general Kepler problem subjected to 
gravitational radiation is due to Peters & Mathews (1963) and Peters (1964), 
who study the case of eccentric orbits using Mathews' (1962) tensor harmonics. 


★ 

Hu also argues that such a behaviour of gravitational radiation reaction 
should account for the cosmological expansion of the Universe. 
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Instead of Eq (5-3), they obtain 



<l£ _ G* (ni l m») l (m l 'Hna) j» , 

“it - * 7 — & — * (e) 

(s-6a) 

where 

r \*mt*w* 

(Kb) 

while V’ is 

the eccentricity and "a" represent' the semi-major 

axis or apse. 

If v/e let 

r, 1 G* 

*© ir c * 

(s-t*) 

the,'. (5-6a) 

is simply 




(sr-ca) 

But recalling our circular orbit assumption, where the enhancement factor 

f (e = 0)=l , P c 

( in (5-6c) is actually 



n L iL T 2, 

*o TT c* 1 Xl 

(S-7») 

which is equivalent to (5-1) (and 3-30 for e=l): 



_ ii » 2LSL P * 

It s r « e • 

(5'7W) 


Ip terms of (5-/a), the Einstein-Eddington-Park power formulae (3-30) and 
(3-31) that we obtained earlier can now be abbreviated as 


jL p( 9> *) = (^[4 oo» 8 * s«?e .m2*] (S-8 Ji 

i P(fl) = 4p 0 [co a J e + i^e]= W 
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The results of Peters & Mathews (1963) reduce identically to these angular 
power distributions in the case of gravitational radiation from a circular 
orbit. The circular Kepler orbit, then, like the spinning rod of Einstein, 
is represented by Figure 1. 

For the more realistic case of an eccentric orbit, the enhancement factor 
f(e) becomes important, as shown in Figure 10a. The power radiated in the 
n^ harmonic. Figure 10b, is 



(S-9aJ 


where g(n,e) is derived from the Sessel functions O n (ne) as ^ ^ 


The enhancement factor f(e) is merely the sum of all of the harmonic con- 
tributions in (5-9b): 

t(«) = £ 

n«i J 

Figure 10 illustrates the important result that most of the power radiated 
is emitted at the higher harmonics for greater eccentricities. 

Before leaving the question of gravitatio al radiation from the freely 
gravitating Kepler problem for two or more bodies, the role of tidal friction 
must be assessed. The treatments presented so far have been adequate only 
for point masses (actually, they are even marginal for point masses - if we 
recall Appendix G). Consequently, the energy dissipation represented by the 
radiation - dt/dt in (5-3) necessarily exhibits itself as a loss in orbital 
angular momentum. However, for macroscopic bodies tidal phenomena manifest 
themselves, providing an energy dissipation mechanism deriving from spin-orbit 
interaction. The resultant tidal friction has the opposite effect on the 
binary system as does gravitational radiation reaction (Appendix S). The 
mechanism of tidal friction, furthermore, cannot be differentiated from 
variations in the velocity of propagation of gravitational radiation, as we 
have already stated (Appendices G & S). 

Gravitational radiation and tidal friction remain in equilibrium until the 
spin angular momentum of the primary is extinguished. Acting as a primary 
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in a binary system, a "black hole" can play host to another companion which 
will radiate away its angular momentum as gravitational radiation, and the 
system will not collapse. Local nonlinear effects can vary the velocity of 
propagation and cause the same result. 


Gravitational B reins s trah lung & Synchrotron Radiation 

A natural extension of the bound (e<l,£<0) two-body Kepler problem just 
treated is the hyperbolic scattering of a small mass m as it passes in the 
vicinity of a large mass M such as a "black hole, at an impact parameter b 
(Figure 11a). On the basis of an electromagnetic analogy, one expects from 
the linearized tensor theory of gravitation that tie small mass m should emit 
a "burst" or "splash" of multipole radiation as it experiences a transverse 
acceleration due to M. 



Such an analysis is presented by Peters (1970) as a generalization of the 
non-relativistic bound orbits treated by Peters & Mathews (1963). Peters 
removes the "slow-motion" assumption employed in the derivation of the quad- 
ruple radiation (3-29), (5-1), and (5-6), and he considers the re’ativistic 
case of high velocity - assuming that the impact parameter is large enough 
(b>>M), which keeps the scattering looking like the classical one of Figure 
11a rather than the general relativistic behaviour at small impact parameters 
( JM<_h£6M ) such as depicted by Figure lib. (The orbital theory of Figure lib 

is due to Forsyth (1920), Morton (1921), and Darwin (1959) who qualify the 

★ 

conditions for capture, stability, and scatter. ) 


See also de Felice (1968) tor a gravitational Stormer-type representation of 
general relativistic orbit theory. Also read Ruffini & Wheeler (1969, 1971). 
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Peters develops his own representation of relativistic radiation at high 
velocities, different from the "fast-motion" approximation which has already 
been discussed. In fact, hi? paper (1970) presents an excellent critique of 
the "fast-motion' method used by Smith & Havas (1965), who treat tne Kepler 
problem by improperly neglecting the stresses ,n cne system. 

In the "slow-motion weak-field" limit, Peters shows that the energy 
radiated in Figure 11a is (-<<1 or vi<c) 

AE. * Tf § tlV^s , *«c (?-»> 

which can be derived from the Einstein quadrupole radiation formula (3-29) 
by calculating the quadrupole tensor Q.. of the mass m from its Schwarzschi Id 
geodesic in Figure 11a. The same result is obtained by P.uffini & Wheeler 
(1969), being discussed further by Ruffini (1973). However, for the relativ- 
istic case Peters' approximation gives 

g*mV 

wVG-ffp 
. _ 6* W*»* 

for the radiated power P and energy At, where the precise coefficients of 
Eqs (5-11) must be determined by numerical integration. The results are 
given in Figure 12, evaluated in the eouatorial plane (o = ~/2 on the angular 
distribution of the bound conic in Figure 1) of the hyperbolic orbit. The 
forward beaming or bremsstrahlung is manifest in Peters' analysis. 

In light of Peters' relativistic gravitational bremsstrahlung, the next 
reasonable question to consider is that of the highly relativistic (high vel- 
ocity) bound conics in close circular orbits about a massive object M, such 
as a Schwarzschild "black hole.” The tighter the circular orbit, the higher 
its angular velocity (neglecting the Roche 1’ and consequently the greater 
its chances of behaving like a rotating "searc., qht," radiating out Peters' 
bremsstrahlung (Figure 11c). Such a mechanism constitutes gravitational 
synchrotron radiation. 

Circular Darwin (1959) orbits are allowed for radii r>3M, although they 
are unstable in the region 3M<r<4M. Consequently, Misner et a_[ (1972) have 
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made the suggestion that gravitational synchrotron radiation should be 

emitted by such orbits, although they are unstable and rapidly plunge into 

the "black hole" under radiative perturbations. Nevertheless, Misner and 

★ 

associates present an analysis of scalar gravitational synchrotron radiation 
based upon the Regge-Wheeler-Zerilli formalism, using standard JWKB methods. 
They obtain the result that particles in high relativistic orbits which are 
coupled to a scalar gravitational field radiate strongly into narrow 
synchrotron beams at high harmonics of the orbital frequency. 

Because of the importance of the synchrotron mechanism ir. the experimental 

verification (the energy of the source can be much less) of gravitational 

radiation, Misner's work on the scalar theory has prompted Davis, Ruffini, 

Tiomno, & Zerilli (1972) and Breuer, et al^ (1973) to investigate the general 

problem of scalar (Spin-0), vector (Spin-1), and tensor (Spin-2) synchrotron 

radiation from r = (3+6 )M circular, unstable orbits. If u=mu is the freq- 

0 3 

uency of the radiation and <o o represents the frequency of the orbit u) 0 -/(M/r o ), 
then the power emitted is given in all three cases (scalar, vector, tensor)by 

in the Regge-Wheeler formalism. The functions RjjJJj and rJ™ are of magnetic 
(odd) and electric (even) parity respectively, being defined explicitly in 
Appendix R. They are computed In the asymptotic region (r -*■ +<») from a 
Green's function technique using the solutions u(r # ) and v(r*) of the 
"Schrodinger" equations (3-66) and (3-67) without sources: 

fe * (.*• -«» * 0 ) 

^ » (.‘-v*). -0 J . 

u must be outgoing at r* = + » and v r,.ust be ingoing at r* = - «° (into the 
"black h^le" at r=2M). The structure of V eff depends upon the field 
(scalar, vector, tensor). 

The results of Davis, £t (1972) are given in Figure 13, which represents 
the power radiated at a circular orbit r=3.Q5M (6=.05). The scalar spectrum 
corroborates that of Misner et aj^ (1972), determined by the JWKB method, but 
the vector and tensor power spectra do not at all behave as do conventional 

^Scalar radiation occurs in the Brans-Dicke (1961) theory. See also 
Morganstern & Chiu (1967). 
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Figure 13 
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synchrotron mechanisms. Synchrotron radiation is enhanced at high multipoles, 
which accounts for the beaming effect. From Figure 13 the vector and tensor 
contributions are significant at lower multipoles with an associated radiation 
at wider beam angles. As a consequence, vector (electromagnetic or gravita- 
tional) and gravitational tenser radiation do n.t concentrate energy very 
effectively into the orbital plane of tie source. 

The critical frequency is defined as 

(J . * <J W\ i x x ^ i 

cr,t • ir(v^9 t? . 

The power spectrum P(^) in (5-12) is found to vary as 

P(«) ** «•>' 5 

for all three spins: s=0 (scalar), s=l (vector), and s=2 (tensor). 

Particle Falling Into A Schwarzschild Black Hole 

Also of wide astrophysical interest has been Zeri lli's (1965, 1970) study 
of the radiative behaviour of a particle falling radially into a Schwarzschild 
"black hole," deriving from his spectral decomposition of the Regge-Wheeler 
and Peters formalisms into tensor harmonics (Chapter 3). This particular prob- 
lem is very much like Peters' relativistic gravitational bremsstrahlung of 
Figure 11a, except that Zeri 1 1 i is treating a geodesic which is captured 
(impact parameter b<3M) - more like trajectory 2 of Figure lib. In particular, 
Zeri 11 i treats a radial, unperturbed Schwarzschild geodesic (Eq 3-68, along 
the z-axis) to characterize his particle (which neglects and avoids the crit- 
ical issue of radiation reaction). This is the case of zero impact parameter 
(b=0 in Figure 11a) . 

Zerilli's equations (the electric, even parity ones in 3-67) have been 
numerically integrated by several authors. Davis & Ruffini (1971) and Davis, 
Ruffini, Press, & Price (1971) investigate the asymptotic behaviour of the 
outgoing burst of gravitational radiation due to a particle falling from rest 
at Infinity - as suggested by Zeri 1 1 i . The results are given in Figure 14 for 
the effective potential, the energy flux, and the tidal stresses of the 
Riemann curvature tensor (these are of interest to experimentalists). Forward 


&-I4) 
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bremsstrahlung is not manifest (Figure '• 4d) as it is in Peters' analysis of 
larger (b>>M) impact parameters (Figure 12). Davis, et a^ find an outgoing 
radiated energy of 

~ 0.0104^) *nc* (s-lfca) 

2 

where me is the particle's rest energy and M is the mass of the “black hole." 
The energy radiated into the "black hole" proves to have contributions at 
all higher multipoles L (hence divergent), but assuming a cutoff L =tt/2 (M/m) 
the total Ingoing radiation is finite: 


which is independent of the mass of the "black hole" (?). (E^+E^) in 

(5-16) is still a small percentage of the particle's rest energy mc^, however 
which is an important result due to the Schwarzschild geodesic assumption. 
This is to say, neglecting radiation reaction (by assuming the radiated 
energy is small enough for a sufficiently small mass m) results in a small 
loss of energy - if you can rationalize away the ingoing divergence. The 
spectrum of gravitational radiation which Davis, et al_ obtain is given in 
Figure 15a for the quadrupole (L=2) contribution. 


Ruffini (1973), however, has addressed the more practical case of Zerilli's 
problem wherein the particle is assumed to possess a nonzero kinetic energy 
at infinity. Solving Zerilli's Equation (3-67) for a Fourier- transformed 
source given by 



where 


* « (i- to) it/dr 
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Ruffini demonstrates that the spectrum does not vanish at low frequencies 


(as it does in Figure 15a for y=1). His results (Figure 15b-d) illustrate, 


This result appears contradictory if one allows Peters' scatt^ing problem of 
Figure 11a to evolve Into trajectory 2 of Figure lib (capture, * ^3M) and fin- 
ally lets b*0. But for b=Q there is no curvature of '< h > trajectory. 
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in fact, that large amounts of energy are emitted there. In Figure 15b, 
his comparison of quadrupole (L=2) contributions for varying > demonstrates 
the dramatic effect of assuming that the particle has a velocity at infinity 
(y> 1). Furthermore, as shown in Figure 15c-d, the gravitational tensor rad- 
iation emitted outward is flat at the lower harmonics for quadrupole (L=2) 
and higher multipoles. This further substantiates the earlier discussion 
(Figure 13) of the inefficiency of tensor radiation as a synchrotron mech- 
anism. There is :.o bremsstrahlung in the sense of Figure 12. 

Umklapp & Gravitational Pi pole Radiation 

The multipole expansions considered earlier contain dipole contributions 
which are usually argued away, based upon the Newtonian principle of equiv- 
alence that inertial and gravitational mass are the same. This means that 

if the gravitational-to-inertial mass ratio is identical for all components 

* 

of a massive system, gravitational dipole radiation vanishes in the center 
of momentum system - due to conservation of linear momentum. Although neither 
the Newtonian nor the Einstein theory of gravitation assumes that mass must 
be positive-definite, introducing the concept of negative mass does not nec- 
essarily change this argument against gravitational dipole radiation (the 
sign of the gravitational constant G is also important). Nor is the question 
made simple it we recall Bondi's multipole representation wherein the mass 

monopole, dipole, and quadrupole are all non-zero and coupled together, as 

o 

seen in the Newman-Penrose conserved quantity D -mQ in Eq (3-75), or in 
Eq (R-15) . 

To certain observers momentum simply need not be conserved, and in such 
frames gravitational dipole radiation may serve a meaningful purpose. The 
quickest example of such a process is found in solid state physics where 
thermal resistance cannot be accounted for by the normal convention of mom- 
entum conservation. The notion of "umklapp," due to Peierls (1929, 1955) 
must be invoked. 

Conventional conservation of energy , ■ amentum in the linear quantum 

theory are represented by (assume o • ^-vectors) 

Hu, ♦ * Hoj 

« IS, ♦ %\ g • * lij ((MOW 


*This does not rule out mass dipole pulsations, however (Campolattaro 
& Thorne, 1970). 
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If you now consider a crystal, and you feed momentum into one side of it 
in the laboratory, then according to the normal process of (5-18) all of 
that momentum has to come out the other iide. There is nothing n the theory 
which accounts for the thermal i nation of phonons; the crystal lattice, hence, 
has infinite thermal conductivity or zero thermal resistance. Eq (5-1 8b) 
cannot explain the known behaviour of crystals for an observer in the 
laboratory center- of -momentum frame. 

Admittedly, this equation for momentum conservation is a linear one, 

implying the need for a nonlinear mechani c m to account for the thermal iza- 

★ 

tion of phcnon populations. But Peie is presents a linear mechanism in 
stead of (5-18b), which works, and it is known as "umklapp": 


k +k. *)!»♦§ , 

where 6 is a reciprocal lattice vector. One way of visualizing umklapp is 
to recall moments in old Western movies when the wagon wheals appear to spin 
backwards. This is an oversimplification, *-Twever, because thermal izati on 
and entropy are real. 

For a periodic lattice, the wave-vectors k^ and (k^+G) are indistinguish- 
able. Normal momentum conservation ( 5- 18b) and umklapp (5-18c) are repre- 
sented in figure 16 for the l st -Brillouin zone i. two-dimensional k-space. 



•Sc 


(a) Normal 


(b) Umklapp 


Whensoever the normal process of momentum conservation (5-18b) exceeds the 


★ 

The M axwell-Dirac theory (second quantization) is nonlinear, for example. 
There, Eqs ( 5- 18a ,b) are inadequate due to photon-photon interaction. 
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boundary of the l st -Bri 1 loui n zone, umklapp ensues and brings it back into 
the zone, as ■‘n Figure 16b. Umklapp destroys momei. urn and actually changes 
the direction of „ne "Poynting vector" or eiergy flow of £q (E- 3a). The 
onset of umklapp also constitutes the onset of an increase in system entropy. 


The notion of umklapp can also serve as a linearization of some of the 
nonlinear effects of nonlinear radiation theory - and hence, its apf' 1 '. 
to gravitational radiation. The nonlinear wave-wave interaction of cullid - ig 
plane waves can be treated in such a fashion (Figure 16b). The nonlinear 
scattering of gravitational radiation can be analyzed this way. And, in 
fact, scattering of gravitational radiation from the background geometry it- 
self (such as the Schwarzschild background) can occur wit', an umklapp mech- 
anism. Lastly, it can be used to analyze gravitational "tails." 


Umklapp does not necessarily have to be gravitational dipole radiation; it 
simply provides for a nonlinear source c momt turn in a linea’ fashion. 

Hence, it has a meaningful place in the astrophysical theory of thermaliza- 
tion due to gravitational interaction. If a deity feeds nonlinear gravita- 
tional plane waves into one side of a crystal (periodic structure) Universe, 
no one expects them all to come out the other si^e. There is also an obvious 
Implication of a lattice structure in geometrization itself, which can be 
characterized by gravitational umklapp. The perturbations h due to one 
gravitational wave give the background geometry a periodic structure, which 
behaves like a grating or moving mirror and scatters othe^ gravitational 
waves ( Br 1 1 loui n scattering, 192?). Thus, there is a Bragg condition for 
gravitational waves. 


Gravitational umklapp, however, does readily treat dipole radiation mech- 
anisms. Two particular problems come to mind: Motion of the enter of ma.s 

in two-body systems (which is known to exist in certain theories), and 
theories which dispel the Newtonian principle of equivalence (mj^). 

Acceleration of the center of mass of the two-body problem is considereG 
to exist in general relativity by Levi-Civita (1937), but this is demonstrated 
to be In error by Rooertson (1938) and Rob.rtson & Noonan (1968). 1 ’w ever, 

many other theories of uravitation do predict such an effect. irg K as 


K 


ir 

c 1 LmaO-c*)} 1 * 


( 5 - 19 ) 


thffi from Whi trow & Morduch (1965) the parameterized v .:..’ r ’-as compare as 


*See also Schild (1362). 
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Cravitation Theory Accel erati op Of Center Of Mass 

Ger^al Relativity 0 

he ustrom (1912) 

Scalar Potential ^ Variable c) 3K or 5K 

Vector Potential - Parameterized (p) pK 

Kustaanheimo (1957) 0 

Birkhoff (1943) K 

Whitehead (1922) k 

Whitehead - Parameteri zed (b) (1-4 j)K 

Poincare (19C4, 1906) 0 


All of the theories which exhibit an acceleration of the center-of-ma'-t 
in the two-body problem also emit gravitational dipole radiation as a con- 
sequence of radiation reaction and conservation of momentum. The background 
metric can be thought of a recoiliny under gravitational umklapp. 

A breakdown in the equivalence principle (mass is a tensor, not a scalar: 
m ,/mc .) also results in gravitational dipoie radiation, in addition to 

Jt 3 i 3 

producing mass anisotropies. Recent work by Nordtvedt (1968, 1969), Dicke 
(1969), Will (1971), Thorne 4 Will (1971), and Will 4 Nordtvedt (1972) studies 
mass anisotropies in the motions of massive, self-gravitating systems. This 
'.iordtveut effect" results in gravitational dipole radiation and gravitational 
yir-’ •■, if the effect exists. In particular, Brans-Dicke-Nordtvedt 
"t .«<les" should emit dipole radiation during collisions. 

Backscatter (Parker, 1972), like gravitational "tails," and gravitational 
shock waves are a iso admirably treated by gravitational unklapp. Umklapp can 
also occur in the periodic lattice of a neutron star. This is pertinent not 
only for gravitational radiation but alsc for the propagation of optical and 
acoustical (sound) phonon branches there. 

Gravitational dipole radiation and umklapo are discussed further in Appen- 
oix G, based upon angular momentum conservation laws in general relativity. 

Gravitational Zitterbewegung 

Standing w*ves are important in physics, for example in characterizing the 
electrical properties of insulators in solid state physics (the electron vel- 
ocities form standing waves). Standing waves are also important in the theory 
of waveguides. Furthermore, they manifest themselves in the theory of gravitation. 

* 

See also Hawking (1972). 
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The standing wave solutions of the scalar, vector, and tensor wave 
equations in radiation theory give rise to zitterbewegung. The assumption 
of both incoming and outgoing radiation at infinity for boundary conditions 
results in advanced (Ritz, 1908) and retarded (Lorenz, 1867) solutions which 
together comprise standing waves with no outgoing radiation at all. The 
classical hydrogen at mi, for example, is stable for standing wave solutions; 
quantum mechanics is not necessary on this basis alone. 

Standing waves likewise occur in gravitational radiation theory. A small 
mass In the presence of such standing waves will go Into a jittering oscilla- 
tion between the nodes of the standing waves and resonate at twice (quadrj- 
pole) the frequency of the radiative source. 

Gravitational standing waves appear in the analysis of cylindrical waves 
by Harder (1958, II) as well as the work of Thorne (1969, III) on nonradial 
pulsations of neutron stars using the Regge-Hheel er formalism. Thorne (1968) 
also presents a rigorous treatment of .iravitational radiation damping in 
terms of a standing wave analysis which concentrates attention on the poles 
of the S-matrix (the complex frequencies **J+i/T at which the incoming rad- 
iation vanishes). For all resonances in the quadrupole standing-wave normal 
modes of various neutron star models, the poles all lie in the upper half of 
the complex frequency plane - which corresponds to damping (and not anti- 
damping) of the exponential exp(i“t). If a particle is ir the presence of 
such an astrophysical object, zitterbewegung ensues and the smaller particle 
oscillates in resonance until it radiates itself away. 


Scattering £ Absorption Of Gravitational Radiation 


A generalization of the vector (electromagnetic) scattering problem to ten- 
sor plane-waves results in an optical theorem and scattering cross-section 
for gravitational radiation. If we consider the scattering of h in (3-33) 
from a tensor scattering center f (an antenna, for example) 





♦ 




> 


(s-to) 


a Gegenbauer expansion of the incident plane-wave part gives in the usual 
fashion (using the Einstein polarization tensor t ) 

* 7*3 — , , V 

« («*«*,- kl*M) ’ 
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which is the optical theorem for tensor scattering. For an antenna with an 
axis of circular symmetry, the antenna cross-section reduces to (see Weinberg, 
1972, or Ruffini & Wheeler, 1969) 


tat 


»siV f 1 

4*1* L(^5 r ^n^i 


sm0 




where '- 0 is the natural frequency of free oscillation of the antenna, ; is 

the total decay rate of free oscillation, r„ is the decay rate of free 

grav 

oscillation due to the re-emission of gravitational radiation, and n = r_ v .,„/r. 

grav 

This cross-section is maximum when the antenna is perpendicular to the 
gravitational wave (e=*/2) . 


The significance of an optical theorem (also known as the Ewald-Oseen 
Extinction Theorem or the Bohr-Peierls-Placzek relation) for linearized tenser 
radiation is that it represents a unitarity of the scattering-matrix. It 
also constitutes Huyghen's principle, which we know is not valid if there 
exist gravitational tails - from our discussion of the Bondi news ^unction. 

An idea of the behaviour of gravitational tensor scattering is represented 
by Figure 17 from Price & Thorne (1969), for a neutron star. 


It is interesting to note that the scattered spectrum of Price & Thorne is 
identical to (5-8) for the Kepler problem and the spinning rod of Figure 1, 
provided L*2 and M=2 for even-parity, quasi-normal pulsations. Vishveshwara 
(1970) also considers such scattering from the effective potential of Figure 
5, maintaining that the scattered radiation contains the signature of the 
Schwa - zschild "black hole." The scattering of scalar gravitons from rotating 
Kerr "black holes" is discussed by Misner (1972) and Press & Teukolsky (1972) 
with the interesting conclusion that the scalar wave is amplified as it 
scatters off the hole. 


In the vector theory of gravitation, c> course, there are the direct ana- 
logues of electromagnetic scattering phenomena such as Mie scattering of 
gravitational (Spin-1) waves. 


Detection & Experimental Verification Of^ gravitational Radiation 

If you recall from Newtonian mechanics that for any nonuniform gravitational 
there exist gravitational gradients, you realize that you can readily disting- 
uish between gravitational and nongravitational acceleration effects. If you 
a^e in free-fall in an inverse-square field (in an elevator, if you wish), a 
measurement with inertial sensors of the gravitational gradient components ir 
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your elevator frame will not only distinguish such accleration from non- 
gravitational ones, but it will give you the value of the field as well. 
Einstein's principle of equivalence, important for correspondence with the 
Minkowski frames of restricted relativity, is valid only in a uniform 
gravitational field ("fictitious" forces arising from the choice of coor- 
dinates and hence g are not questioned here). But su<*h uniform fields 
are not known to exist in astrophysics; hence his principle is a tautology. 


The Riemar.n curvature tensor R , is essentially a glorified version of 

■Aw. & v 

the same thing, the .Newtonian gravitational gradient or tidal component in 
classical mechanics. Instead of taking the gradient of the scalar field 
intensity with multipoles (quadrupoles, etc.) in Eq (2-3), tidal stresses 
can be determined from the R 1 ^ components of the Riemann tensor (Eq 1-2 
of Appendix A) in the tensor theory of gravitation. These tidal forces or 
Riemann stresses lie at the basis of the construction of gravitational 
radiation detectors. To first order in the linear perturbations h of g , 
Appendix A determines that 


R 


0) 

VP 






(5-M) 


From Einstein's quadrupole radiation in (3-26) and (3-27), which lead to 
the linearized quadrupole power formulae (3-29) and (5-1), the non-zero 
Newtonian ti Je producing components due to a gravitational plane wave prop- 
agating in the X^-direction are caused by hgj and (l^-hj^) in (3-26) or 
(3-37). From Eq (5-23) above, using '3-24), 


'-*£.!•« - £(<VcQA 

_ CS-M) 

R*0»0 * • 

Note that has zero trace (R^q+R 2 q 2 o + ^ 03 o = ^ *hich means that the 

quadrupole detector Q. . responding to the Einstein-Fddington-Weyl plane-wave 
(in Chapter 3) maintains a constant volume under the Newtonian tide compon- 
ents of stress. These components of transverse stress simply deform the 
object by exciting its quadrupole moment Q... A cylinder (Figure 18a), for 
example, has two transverse modes of polarization, separated by 45 due to 
quadrupole symmetry. A disk (Figure 18b) also exhibits a radial stress 
which is indicative of a scalar component of radiation. No detectors have 
been designed for measuring stresses due to electromagnetic (vector) theories 


*Einstein (1911) did not understand this when he first established the principle. 
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of gravitation, nor is there consideration of dipole mechanisms. 



Sam* as (a) 

(B) SCALAR (C) VECTOR 

FIGURE 18 



Experimental work on gravitational radiation has been primarily fostered 
by Weber, et al (1960-1972), who measure responses of a gravitational 
quadrupole antenna to supposed radiation sources. Although they have 
found coincidences between their two cylindrical antennas, these measure- 
men*i' have not been duplicated. The unusually large flux (10 ergs/cm -sec) 
given by Weber's measurements (Figure 19) Indicates unphysically short life- 
times of astrophyslcal objects such as the Galaxy, and consequently a great 
deal of theoretical work on tensor synchrotron mechanisms Is a direct re- 
sult - as we have seen. One should also note in Figure 19 that the coin- 
cidences occur not only when the galactic center Is In view but also 12 
hours latar when it is occulted by the earth. Others (Beron & Hofstadter, 
1969) present argjments that high-energy cosmic rjys can excite phonon 
oscillations of quadrupole symnetry in gravitational antennas (cr telescopes) 
and be mistaken for gravitationally induced oscillations (which are of the 
order of a nuclear radius!). Some work on the generation of gravitational 
radiation In the laboratory is available (Weber, I960; and Forward, 1966), 
although the quadrupole power factor of 10"®^ erg/sec in Eq (3-29) is more 
than enough to discourage the most ardent experimentalist. A good survey of 
experimental work Is presented by Logan (1973). 

Attempts other than synchrotron mechanisms to explain the high fluxes 
measured by Weber Include the amplification of gravitons upon scattering 
from rotating Kerr "black holes" (Misner, 1972; and Press & Teukolsky, 1972! 

*This suggestion is rebuffed by Ezrow, et al (1970), however. 








as well as the collision of black holes which cannot bifurcate (Hawking, 
1972). Hawking's mechanism is reminiscent of using a power saw tc cut 
saplings. 
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But rather than build detectors, some observers have been using natural 
gravitational antennas - for hundreds of years. Hooke (1671) considers the 
attraction of two spheres, that is, gravitational attraction itself, as a 
"gravitational antenna" indicating the existen e of gravitational waves in 
order to explain the mechanism of gravitation. Laplace uses the planetary 
orbits of the solar system as a means of determining perturbations caused 
by finite velocities of propagation of fields, and hence of waves. And now 
experiments are being performed using the quadrupole modes of the earth, 
the moon, the earth-moon system, and the Mariner 6 spacecraft as gravitational 
antennas. Dyson (1969), in the process of considering the seismic response 
of the earth to gravitational radiation, makes the aiggestion that coinci- 
dences might be autocorrelated between seismic data and astrophysical pul- 
sars. This suggestion gives negative results from some authors (Wiggins 
& Press, 1969), but an interesting flury of experimental work is being 
carried out by Dror Sadeh (1972), who claims to have correlated terrestrial 
and lunar seismic activity with pulsar CP 1133 - much to the consternation 
of seismologists such as Mast et al^ (1972). 

In order to make this discussion of the tidal stresses complete, it is 
worthwhile to conside* the Riemann stress components in the Regge-Wheeler- 
Zerilli formalism For a wave propagating along the X^-axis as before, the 
asymptotic tide-producing, transverse components are 

R*,*w jywAr 

SM >• 

where 

VJjjMl* (& - « t8 » - 5?5 (*»» 

which appears in the tensor harmonics of Appendix K. The physical inter- 
pretation of these stresses is found by referring to Figure 14h. Tidal 
acce 1 -itions are also discussed by Finley (1971). Tidal deformation is 
trt j by Manasse (1963). 
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Astrophysical Miscellania 

There are numerous sources of gravitational radiation in astrophysics 
which have been examined with the theoretical results of Einstein's nonlinear 
tensor theory of gravitation. The nonradial, asynmetric pulsations of stars, 
in particular neutron stars, have been studied quite extensi by Thorne 
et aX (1967-1970), using the familiar Regge-Wheeler formal’! bu* applied 
to general relativistic hydrodynamics. The odd-parity perturbations des- 
cribe differential rotation, while the even-parity perturbations desc’ibe 
the pulsations. Thorne's work is a thorough treatment of the damping of 
stellar oscillations under gravitational radiaticn and a study of the as- 
sociated relaxation times. However, it is limited in that Ruderman (1968, 
1969) has presented moaels of neutron stars with crystalline mantels which 
support shear stresses. Ipser (1971) extends these investigations to 
slowly rotating neutron stars. 

Pulsating stars prove to have another significance in the astrophysical 
theory of gravitation, as pointed out by Morganstern & Chiu (1967). In the 
Brans-Dicke (1961) theory a transverse scalar (Spin-0) component o' the 
radiative Riemann tensor is predicted, whose source is a monopole oscilla- 
tion (see Robinson & Winicour, 1969). For a radially symmetric pulsation, 
scalar waves can be radiated away but not tensor waves. Such scalar radia- 
tion, furthermore, damps out the radial pulsation of the neutron star in a 
matter of seconds. The existence of the scalar field can hence be ruled out 
if neutron stars are found to pulsate. 

Gravitational collapse with asyrmetries is another aspect of gravitational 
radiation theory which has astrophysical importance. Symmetric collapse, of 
course, does not radiate tensor radiation due to Birkhoff's theorem. Asym- 
metric collapse is treated by de la Cruz, Chase, & Israel (1970) who deter- 
mine that ingoing radiation (down the Schwarzschild "black hole") causes 
external asymmetries to decay quickly with a brief relaxation time on the 

3 

order of 2GM/c . They conclude that "black holes" cannot serve as tensor 
synchrotron mechanisms, a result consistent with Chandrasekhar's tri axial 
Jacobi ellipsoid of Figure 9 wherein asymmetries are quickly radiated away 
c ing collapse and the "black hole" becomes a Maclaurin spheroid. The 
radiation of such asymmetries is also discussed by Israel (1971), Penrose & 
Floyd (1971), and Penrose (1972). 
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Excited "black holes" are discussed by Fackerell (1971), Peres (1971), 
and Goebel (1972). They are also relevant in the context of "black hole" 
collisions, wherein the gravitational radiation is entirely tensor radia- 
tion according to the analysis of Hawking (1972). This is simply because 
Brans-Dicke "black holes' nave no scalar field Hence, no Brans-Dicke scalar 
radiation can be emitted during the collision of two "black holes." 

In the many and varied efforts to study the behaviour of nonlinear tensor 
gravitational radiation, low frequency and high frequency approximations 
have been made. Low frequency work, of course, is what leads Ruffini (1973) 
to the important conclusions of Figure 15. Likewise, low frequency analysis 
is the basis of Bergmann's (1971) investigation and of the Mariner spacecraft 
tracking studies of Anderson (1971) (see also Gibbons, 1971) in an effort to 
observe scintillations due to gravitational waves. High frequency approx- 
imations, on the other hand, are discussed by Isaacson (IS68) as well as 
Halpern (1971), Isaacson & Winicour (1972), apH Kafka (1970). The signifi- 
cance of Isaacson’s work has already bee" stated in the discussion of his 
resultant stress-tensor. Kafka assumes the possibility that gravitational 
radiation might be of primordial origin. Primordial gravitational radiation 
is further pursued at very long wavelength by Rees (1971) and critiqued by 
Jackson (1972). This work all emphasizes the cosmological significance of 
gravitational radiation. 

Misner's (1968) mixmaster univer-;,; is another interesting embodiment of 
gravitational radiation in relativistic astrophysics. It essentially treats 
the Universe as a gravitational waveguide, much like the work of Campbell & 
Morgan (1971) which employs a gravitational Debye potential formalism and 
the work of Marder (1958, II) on cylindrical standing waves. Lastly, Lapedes 
& Jacobs (1972) maintain that tachyons should produce gravitational Cerenkov 
radiation, although tachyonic systems a^e not very popular among astro- 
physicists. 

Meditations Upon A Few Thoughts From Classical Physics 

There are a number of topics in the classical astrophysics of gravitational 
radiation which must not go unnoticed. They are controversial. And yet we 
know that controversial thought : . traditionally wherein lie many of the 
great advancements in human s'ience. Consequently, it is a greater heresy 


* They can be excited by a particle falling into the geometry, as in the 
Zerilli problem. 
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to ignore such though+s than to weigh their merit - although many of them 
should not be dignified by a great deal of detail. 

The concept of mass is fundamental to any theory of gravitation and hence 

to any theory of gravitational radiation. But is mass positive-definite, 

or does Nature provide for states of negative mass (Foppl , 1898) and negative 

matter or anti -matter (Schuster, 1898)? We know, for example, that there 

exist negative energy instabilities in nonlinear plasma theory (Sagdeev & 

Galeev, 1969). Can the Dirac theory of negative energy states be re-inter- 

★ 

preted, perhaps, in terms of negative mass (using the negative mass shell 
of p ),and if so are there such things as negative energy photons and 
gravitons (e.g. Misner, 1972) associated with jumping the energy gap? Such 
conjectures have been made in explaining the behaviour of quasi-stellar 
sources. Indeed, in view of the negative mass solutions manifesting them- 
selves in the only known c.. -t radiative solution of Einstein's nonlinear 
theory, we are forced to address this problem. It is no longer a metaphys- 
ical whim of Arthur Schuster. 

It is likewise believed that mass and energy are equivalent in relativity 
theory, as stated by Eq (3-4). This assumption results in an epistemological 
contradiction with the related postulate of classical electrodynamics that 
charge is invariant. Gravitational radiation is produced by mass, and it 
transports mass (according to Einstein's theory), while electromagnetic 
radiation is produced by charge - but does not transport charge. This 
irreconcilable difference is just cause for rejecting any grounds for an 
electromagnetic analogy in the linearized tensor theory of Einstein. It is 
an epistemological flaw of the first magnitude and deteriorates any founda- 
tion for a unified field theory based upon geometrization. The linear 
Maxwell theory and the nonlinear Einstein theory of energy transport cannot 
be reconciled. 

Also, we must note that negative energy densities characterize all gravita- 
tion theory. There is nothing intrinsic in the scalar, vector, or tensor 
theory that denies this. Time-symmetric arguments (Brill. 1960) for positive- 
definite energy are ta itological in the sense that they assume a reversal of 
the Poynting vector (g. Q * -g^) under time-reversal (reversal of incoming and 

outgoing radiation conditions). They assume what they attempt to demonstrate. 

_ 

Velocity is anti -para 11 el to momentum f or negative mass. The helicity j(>/p 
of anti-particles is directly accounted for If anti-particles possess negative 
mass. However, see Morrison & Gold (1957) and Schiff (1958) on this point. 
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Furthermore, mass is assumed to dilatate as stated by Eq (3-5), in the 
Einstein theory of relativity. The structure of gravitational radiation 
theory must not collapse just because this single postu'ate p^ov-’s to be 
untenable if there are future developments in electrodynamics. Likewise, 
there Is a "second postulate of general relativ'ty." In other words, con- 
vention has it that the velocity of a photon is independent of the velocity 
★ 

of its source. If the nonlinear tensor theory of gravitational radiation 
Is allowed to be linearised and characterized witn all of the other attri- 
butes of electromagnetic radiation, then must not gravitons likewise be 
independent of the velocity of thel * source? In contrast, there is a 
ballistic theory of gravitational radiation. The Ewald-Oseen Extinction 
Theorem (the Optical Theorem) renders the second postulate indeterminable 
by experimental method for both photons and linearized gravitons. 

Is a gravitational wave stable? The contribution of the Regge-Wheeler 
(1S57) study of the stability of the ihwarzsch Id metric is the implication 
that no physical theory can be considered seriously until it has been demon- 
strated to be stable. Whence, if the integrity of n^ v is questioned then 
why not question the stability of the gravitational waves h ? The- regen- 
erative nature of self-interaction in the nonlinear theory of general 
relativity is important in such considerations, 

I* the gravitc.i stable? This aspect of the question of stability is taken 
from the quantized linear point of view rather than the general nonlinear wave 
and resonance Instability discussed above. It becomes Important If experi- 
mentalists fail to detect such radiation or if there exists a "cosmological 
redshift" of the gravitational radiation power sp^-trum (which there should 
be in order to avoid an Olbers-Ches aux-Halley paradoxl. Decay mec'-anisms 
for gravitons, with and without rest-mass, can be readily iarovised. The 
dispersion relation for vector (Spin-1) gravitons Includes an effective mass 
term from the Proca theory or from the propagation of gravitons through matter 
using the Weyl tensor (in direct analogy with "plasmons," which are photons 
of effective mass propagatl"" through an electron gas). The vector (Spin-1) 
massive graviton is readily extended to the tensor (Spin-2) theory by addin' 
Proca terms to the tensor wave equation. 

Lastly, we mention that there is w ,rk in the literature on compo^ . j . ield 
theories, such as the neutrino theory of light. A candid discussion of 
these Is presented .;1thin the context of gravitational radiation theory in 
Appendix V. 

*This postulate has marginal experimental basis (Wilson, 1972). 
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Nonlinear Effects In Astrophysics 

Several aspects of the nonlinear nature of the tensor theory cf gravita- 
tional radiation are significant and must je kept in mind as one interprets 
the astrophysical properties of the Universe. The first characteristic of 
Einstein's geometrized theory of gravitation which asserts itself is the 
remarkable feature that it has no radiation reaction problem, as does the 
linear, classical theory of i.axwell and Lorentz. This is originally due 
to Einstein's geodesic postulate, but actually it is a consequence of the 
fact tha f the EIH equations of motion follow from the nonlinear field equa- 
tions because of (3-38). Referring to Figure 20, an empty Schwarzschild 
geometry is traversed by a geodesic on the background n • 



P— 

/ 

Figure 20 

However, upon inserting a particle of small mast m into tl.e Schwarzschild 
geometry, the particle does not follow the geodesic f> • n^ v but rather that 
for n + h The perturbations h are those imp".tu upon the st'ucture 
of the background by m and are determined by the very existence of m in the 
geometry. It you can aver solve Einstein's nonlinear field equations you 
get as a d^re.. dividend the new geodesic. Radiation reaction is automat- 
ically accounted for by the nonlinear nature of the theory. Linearizations 
of Einstein's theory, however, such as the Regge-Wheeler-Peters-Zerilli 
formalism are confronted with the problem of radiation reaction, because they 
neglect it. Zeri Hi's treatment of the particle falling in, if you recall, 
assumes the particle follows n •• which is not tie trajectory it actually 
takes. Only for vanishingly smaV mass ,..»0) is the li earized treatment 
tenable. 
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The regenerative feature of Einstein's nonlinear field e^ur ions is 
another crucial aspect of thair properties. The fields act as their own 
sources (Gupta 1957; Thiriing, 1959, 1961). 

Having emphasized the nonlinear nature of the tensor radiation theory, 
we must finally ask ourselves what it means to treat the tneory in a linear 
fashion, as with energy and momentum relations like Eqs ( 5-1 8a ,b) f~om 
quantum theory. Surely we can write for en' rgy 

E - h (fcv) > 

but what does that mean? From the work of Pauli & Fierz (1939) and Gupta 
(1954, 1962) tensor radiation is Spin-2, in e linearized theory of grav- 
itation. However, this does not mean tlr - Sp* i-2 "g-. avitons" exist in 
Einstein’s nonlinear theory, for it nas ne.°r beer quantized. A casual 
look at the behaviojr of colliding plane wav?*, (Szek.°res, 197c, 1972; ana 
Penrose & Kahn, 1971) and radiation scattering flor'ence & Jam’s, 1967) 
illustrates the naivete one must have in c der to consider energy and momen- 
cun in the normal linear fashion (superposition) 

Qj V la)^ — Oj (SHSfl) 

Vl, + = fcj . (s-IBb) 

The self-interaction (See e.g. Terrence & Couch, 1970) o* non’inear radia- 
tion and the behaviour of gravitational "tails" vDeWitt A Brehme, 1760; 

Couch el al_, 1968; and Hallidy & Janis, 1970) also are formidable problems 

- although significant progress nas beer made since the advent of Bondi's 
asymptotic approximation (Rondi et aj^, 1962) and cf the Newman-Penrose 
tetrad formalism (Newman & Penrose, 1962). 

We can invoke gray 1 'all onal umklapp (5 18c) but we .uusc a c '' ourselves 
questions about gravitational dispersion of the frequency v in Eq (5-26). 
Gravitational radiation must inverse interstellar distances through a 
"medium" of nonlinear wave-wave interactions. It would be foolish to sup- 
pose that the frequency detected by a terrestrial observer •/ would be the 
theoretical one for quadrupole radiation (3-29), namely: v'~2’. For this rea- 
son optical and radio (see Charman e_t aj_, 1970) correlations with invita- 
tional radiation detec.- coincidences arc of major significance bec.u-. 1 ^ney 
can help ascertain the dispersive nature cf i’ltergalactic space t; ora ■ i ta- 
tional radiation. 
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The failure to find any electromagnetic and gravitational coincidence* 
so far is a basis (assuming the gravitational data are real) for maintaining 
that there is a high degree of shift in frequency - meaning that the freq - 
uency of the source v (and hence its angular momentum) can be much less than 
that measured here on earth. Synchrotron mechanisms are then unnecessary 
in order to rationalize Weber's data. One over-estimates the energy of the 
source because he misinterprets the energy density of tie gravitational wave 
- due to the linearization of the theory. 

Another way of presenting the problem of ronlinear radiation is to turn again 

** 

to Eq (5-26). What is Planck's constant "h" for gravitational phenomena 

(which have never been quantized in the tensor theory)? Is it constant? 

S’-^ply be increasing the ratio h_„ /h„ one can easily account for apparent- 

y r a v Gfii 

ly high quadrupole fluxes. 

And what about the uncertainty principle? Is there an uncertainty principle 
for gravitational phenomena? 

AX AP > t> ( Sl7> 

Again, what is h? How do we know that we can even measure gravitational 
radiation? 

Indeed, one must not abuse relationships for the Maxwell theory and the 
quantum theory which have not yet been established for gravitation theory. 

The closest thing yet to gravitational quantization is the tensor harmonic 
decomposition of the linearized field equations. 

With this we end our query on the nature of gravitational radiation theory. 


Electromagnetic and gravitational coincidences are further complicated by 
the assumption that their velocities of propagation are equal - which may 
not be true. 




To demonstrate the epistemological inconsistency In gravitational and quan- 
tum physics, one can use the identical Bohr quantization condition employed to 
quantize the orbital angular momentum of the Coulomb Kepler problem, and derive 
h for the gravitational Kepler problem (the solar system). On the basis of the 
orbit of Mercury (this is really the second "Bohr" orbit If Bode's law is rel- 
evant), ^mva *1039 joule-sec, rather than lO" 3 ^. In making this analogy, one 
sees another instinctive feature of the gravitational quantization problem: 
charge is quantized and so is its field theory - while mass is not, neither 
is its gravitational field, at the present time. f\ 
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CHAPTER VI 


CONCLUSIONS 


Gravitational radiation theory has been discussed from the point of view 
of its scalar, vector, and tensor counterparts. This discussion has not 
ventured to ascertain the existence of such radiation in gravitation theory 
but rather to take a candid look and to explore its meaning in the context 
of contemporary physics. 

We have found that there Is a theoretical basis for such radiation in the 
vector and tensor theories of gravitation, but in doing so we have been forced 
to address all of the fundamental nuestions about radiation theory. These 
Include the notions of radiation reaction, the behaviour of advanced and 
retarded solutions, the choice of boundary conditions, the transport of 
energy, the meaning of energy, the nonlinear nature of radiation - and even 
the meaning of a wave. We have made analogies with classical radiation theory 
which are sound In the linear electromagnetic (vector) theory of gravitation, 
but they become unsuitable when extended to the nonlinear tensor theory of 
general relativity. There they have as yet no experimental basis. Nonlinear 
electromagnetism, on the other hand, does have some foundation in experimen- 
tal fact, such as Maxwell -Dirac spinor electrodynamics and nonlinear optics 
(both employ nonlinear constitutive relations 2-14). This in turn provides 
reason for analogy with a nonlinear vector theory of gravitational radiation. 
Nevertheless, we know and understand very little about nonlinear gravitational 
waves. 

To be sure, superposition must be sacrificed in the nonlinear theory. The 
linearization, furthermore, of the nonlinear field equations inherently des- 
troys any possibility of accounting for the interaction of the gravitational 
field with . tself. 

The exact radiative solutions of the tensor theory of gravitation likewise 
do not yet render an adequate understanding of the problem. The oily exact 
solution for moving sources in general relativity allows for negative mass, 
which like the concept of negative energy densities, may be very metaphysical 
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in nature. Indeed, enchantment with nonlinear tensor theory teaches us 
more about the tyranny of mathematics than it reveals about the physical 
nature of gravitational radiation. 

However, by means of linearized tensor radiation theory, we have studied 
the behaviour of mass quadrupole oscillators such as the Jacobi ellipsoid 
and the two-body problem. As paradigms in astrophysics, they have taken a 
useful place in the life and death of neutron stars, binaries, and "black 
holes." 

By means of linearization, we have also gained insight into characteristic 
features of the tensor theory which can eventually lead to a decision as to 
which of the scalar, vector, and tensor theories is most consistent with 
experimental data. 

Perhaps there has been undue interest here in the vector theory of grav- 
itational radiation, but it is intimately consistent with classical physics, 
quantum mechanics, and the quantum theory of radiation - which the nonlinear 
tender theory is not. There isJMIk a prodiguous amount of comparison 
of physical phenomena with electromagnetism - and there always will be. 
Accordingly, an understanding of the real nature of the vector theory is 
of fundamental significance. 

Nevertheless, this query on the nature of gravitational radiation has led 
us to understand better the characteristics of linear and nonlinear radiation. 
With a little bit of luck, we may even be on our way towards a unified 
picture of the nature of the physical world. 
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APP0CIX A: LINEARIZATION OF RICWNIAN GEOMETRY 


In an arbitrary geometry one is concerned with its affine connections. If 
the geometry is Riemannian, these affine connections are represented by the 
Christoff el symbols 

r; - W * ifl w i'f/ wl • M 

The Riemann tensor for this geometry is defined as 

RV»r;.-r^r;r;-r;r;. »-» 

The geometry is "flat" if R a =0. When contracted (a=$), R a gives 

ayBv yav 3 


the Ricci tensor R 


r sR - =p , ‘ -r" +rr c -r r 


For a "weak-field : approximation of the metric tensor 




(fc-4) 


then to fi^st order in h=h x =n' h a , the two right-hand terms in (A-2) and 
(A-3) are negligible, giving for (A-2) 

* l [ + 


and instead of (A-3), 


qW « d* - p* -P* (A-ft 

Evaluating (A-6) by taking the appropriate partial derivatives of (A-l) gives 

1$ * i[ Q v VI M 

to first order in h. This may be re-written as 

R® * i [ 0 (h^- (** . (A-8) 

“ y iC-* kh lW-W-VM , % 


Similarly 


to second order in h. 
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APPENDIX B: COORDINATE CONDITIONS & GAUGE INVARIANCE 


Gauge Invariance 

If one considers an infinitesimal coordinate transformation 

« x* «■ /w (B-o 


the metric tensor is changed as follows 




Because, furthermore, g yv »n uv -h u (raising indices) then 

A* 




Hence, if h is a solution of the wave equation (A-6) 


uv 


R® = 0 

”/»y 


(BrD 

(B-3) 

(•'4 


then so is 



This demonstrates the gauge invariance of the field equation (B-4). 
the perturbation of h in (B-5) is null 

1JV 




o 


If 

(B-4) 


one has what are referred to as "Killing's Equations." A gauge transformation, 
then, is comprised of inhomogeneous Killing vectors. 


Einstein Coordinate Condition 


Einstein (1916, 1918), in his original derivation of the gravitational 
radiation associated with linearized General Relativity, chooses to define 
a tensor wave function 


r* iirifki (M 
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subject to the following coordinate condition 

y* V = 0 (8- 8) 

B .V 

It appears to ue a common misconception that r< ation (B-8) is due to 
de Donder (1921), and the latter is often associated with his name. It 
should be more properly credited to Einstein. 

The virtue of the Einstein coordinate condition (B-8) is realized if one 
recalls the first-order linearized Ricci tensor in Appendix A. With (B-7) 
and (B-8), Eq (A-6) simplifies directly to 

- toy • (M> 

If y uv does not satisfy (B-8), then one can always find a y' uv and hence 
h' wv which do, by performing a gauge transformation (B-l) where Dc p =y uv 

V 

Harmoni c Coordinates 

Another coordinate condition much like (B-8) due to Einstein is the 
"harmonic coordinate" condition employed by de Donder (1921) and Lanczos (1922), 
and later exploited by Fock (1939, 1959). It is defined as 

rS/r^-o , 

which is not generally covariant and which thereby destroys the general 
covariance of the theory. It is "harmonic" because the coordinates themselves 
satisfy 0X^=0 due to (B- 10) . Defining a metric tensor desity § uv =/(-g)g uv , 
there Is the condition 

r„ * ° . <*''> 

equivalent to (B-1Q). The difference between (B-ll) and (B-8) is a very 
academic one. 
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APPENDIX C: ENERGY MO^IfrUM PSEUDOTENSOR 


The energy-momentum pseudotensor is demonstrated by Landau & Lifshitz 
(1962) to be 


M‘V(^r,;-rr,r^)}. <m 

Defining a metric density § UN W-gg vv , this is also 



In the weak-field linear perturbation approximation, however, the pseuao- 
tensor t can ! c found more directly than by Eqs (C-l) or (C-2). The pro- 
cedure is mentioned, for example, by Weber (1962). Multiplying the wave 
Eq (3-16) 


a V 21tT, v lM6) 


by h‘ e , one arrives at 



The covariant divergence in (3-3) of the stress energy tensor is 
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T p „ * (T> . t'V a 


Furthermore, the stress energy tensor density T 1JV and the stress energy 
pseudotensor density t uv are related by 




fc-5) 


From (C-5) and (C-4), then 


►V ' ' »v 

A 


t 1 ”.. ■ -wr's-ik/T* 

«VV- , V. 


Hence, (C-3) can be re-stated, with (C-6), as 


H^O^= *4H t/. 


(c-« 

(c-i) 

(c-») 


Expanding the left-hand side of (C-7), one obtains 

Direct comparison of (C-7) and (C-8) - recalling Eq (B-7) for y yv - gives 
the first-order stress energy pseudotensor density t pv =/-g t pv : CC"®) 

V* iK't'V' *■ V' 1 ?M> 1 - 


For the particular case of a perturbation h propagating in the direction 
(discussed in Chapter 3, Eq 3-25), and noting that /-g-1 for a flat back- 
ground the (1,0) component of t * 1 in (C-9) is 


K'v'V'V' >l . *• 


t°= J 

*• A*Fj 

Whence, the gravitational energy flux (Poynting vector) is determined by 




The first-order pseudotensor t y v in (C-9) provides the basis for studying 
radiation spectra and energy radiated in this approximation. 
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APPENDIX D: CONSERVATION OF MOMENTUM & THE EQUIVALENCE OF 
GRAVITATIONAL & INERTIAL MASS 


Angular Momentum & The Pseudotensor 


Having addressed the Landau-Lifshitz pseudotensor in Appendix C, Eq (C-l) 
and (C-2), we need to assess its relation to conservation laws. The conser- 
vation law 


(b-o 


is satisfied by the 4-momentum P p 


which can be written in the form of a 3-space integral 


(B-2) 

(D-3) 


Because the quantity in the integrand of (D ?) is symmetric in t;ie indices 
u,v we obtain a conserve tion law of angular momentum 

f tw> 

M fV =j(xV-xV) = 


for a closed gravitating system in General Relativity. Because of this 
conservation condition, one can define a center of inertia system moving under 
uniform motion. The 'enter of inertia is defined by M 01 = const in (D-4), 
whereby its coordinates are determined by 

jCf%trx-j)dv . 

By performing gauge transformations such as (D-5) one can create or annihi- 
late the pseudotensor t pv , which characterizes the energy and angular momen- 
tum of gravitational waves - the argument of Scheidegger & Infeld. However, 
t wv does not vanish globally under such transformations. 
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Equivalence Of Gravitational & Inertial Mass 

Landau & Lifshitz consider the purely gravitational case (T' yv =0) and 
define . 

■ fcMrt'-fYl fr-o 

whereby P y in (D-2) is an ordinary 3-space integral (i“l,2,3) 

• ( ®' 7) 


Demanding that the metric g^ v must be asymptotically flat (Mlnkowskian) at 
infinity, they evaluate (D-6) for the metric (where M is the total mass of 
the system; n is a unit vector along jr) 



which gives 
and (D-7) becomes 


P 



K° 



i 00i_ Me’ 

h 7 

(D-*$ 

0 

P°= Me . 

(P-9) 


In the center-of-inertia asymptotically-flat frame (D-8a), 4-momentum Is 
determinable and it is conserved. This is a statement of the Newtonian equi- 
valence of gravitational and inertial mass discussed in the text. There is 
no outgoing or Ingoing gravitational dipole radiation at infinity because 

~ P* = 0 6n») 


f. 




from (D-9). 

Failure Of Conservation Laws 

One cannot always formulate a principle of conservation of angular momentum, 
whereupon the notion of a center of inertia loses its mathematical basis. 

Such is the case for certain interactions, or if the integrand in (D-2) is not 
symmetric in the Indices u,v - which is precisely what happens when one 
Includes the Electromagnetic stress-energy tensor (T' yN 7o). Landau & Lifshitz 
determine (p. 311) that (D-2) becomes 
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p ,= # 

This result is relevant to an inequivalence of gravitational and Inertial 
mass. Eq (D-ll) is an expression for gravitational umklapp. 



101 


APPENDIX E: POST-NEWTONIAN GENERAL RELATIVITY 


We are concerned with the equations of motion (3-45) 

- - r£ tftf (3-45) 

® («> f 

In the EIH approximation of Figures ?. & 3 to some order e n * c" n . The equa- 
tion of motion (3-45) represents a geodesic defined by the affine connections 
r£ g which can be expanded as follows: 

tf * -r:-2r;x'-r.‘xY +[rl+2r;^r;xV]x i . 

(f.n >") n (O m m («-*) 

The powers of E on the right-hand side of (E-l) reflect the order to which 
the Chrlstoffel symbol must be expanded In g in Eq (A-l). Note that each 


time-differentiation (the dots) lowers the order of c by c . 

2 4 

Newtonian Approximation 0(e , * ) 

This is the left-hand line of Figure 2 in the text. Analyzing (E-l), to 
order 0( E ^,*^) there are the following contributions: 

e l : r£ 



• 4 


Hence 




Correspondence with Newtonian mechanics gives, then, 

h* ** % 

whereby from (A-4) g *•(!♦' !♦/<?) 


(s*l) 
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Post 0 ,5 - Newtcnian Approximation 0(f£,A£) 

3 5 

As depicted in Figure 2, there is no contribution at order 0(e ,x 3 ). 

This is due to the conservation of energy and momentum. 

Post 1 ’ Newtonian Approximation 

The post^ '^-Newtonian forces manifest themselves as the second line of 
Figure 2. They are comprised of the metric contributions g 00 , g ok> and g^. 
Frcm the equation of motion (E-l), we require 43 2 


e 4 : r£ 

A rr,.r; 


e*-. 

e‘: r; p 


For the Hicci tensor (a- 3) to this order of e , one has 


Calculating the Christoffel symbols r u g to these orders of c’ n using (A-l): 

j . {*14 ■ tfVVv-1 

( Pi • 

; | ;i 

6 : {P 4 *° 

For the Ricci tensor (A-3) to this order of e , one has 

(E S) 

Plugging (E-4) into (E-5) we get the following Ricci equations: 

$„= it f„ M 


(E-s) 
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Transforming (E-6) to harmonic coordinates (B-10), the Ricci equations 
simplify to 

it* 

«... rfj. 

«■?. < E, > 

«4 a tVJt, 

Einstein's field equations now follow from (3-14), using (E-7). Taking 
caution with the c"^ factors in the stress-energy tensor T yw , substitution 
of (E-7) into (3-14) gives to order 0(e 4 ,x 6 ): 



From Eq (E-8c) we arrive at the Newtonian result in (E-3) 

The Newtonian potential 4> is that defined by Poisson's Eq (2-6) 

V*4> = 4v6 f - 

g^ v r„jst be asymptotically flat. Hence, in order that^* Q vanish at infin- 
ity, the scalar potential * is identical to (2-8) 
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The solution of (E-8b) is the vector potential A.. * : 



(e-») 


which is identical to that in the electromagnetic theory of gravitation 
discussed in Chapter 2. 


In order that g-- vanish at infinity the solution of (E-8d) is 

" J 


% » 


(E-n) 


The remaining field equation (E-8a) is solved using (E-3* ) and (2-6') 

If we assume (Weinberg, 1972, following Miller, 1962) * 


then (E-ll) becomes 


* -*(£)-*+ , ^ 
if - AMP “I (**> 


which has the scalar solution 



(e-«J 


The harmonic condition used earlier imposes the following relation be- 
tween the scalar potential <t and the vector potential A: 

£^+V-fc= 0 . CEO 
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Using (E-3‘), (E-ll), and (E-13) in conjunction with the solutions 
♦ , A, and the Christoffel connections (E-4) become 

h - 

r; • ilwl-MH 

P; * 

Pi - (E-17) 

Pi * . 

Substituting (E-17) back into E-l) we obtain the "post^ Newtonian" 
equations of motion: fp-ttf 



POST-tEKTONIW EQMTICN OF rtUIQN 


Note carefully that these are simply scalar and vector potentials. Further- 
more, note that Eq (E-18) contains the Lorentz equation of motion (2-15) 
from the electromagnetic theory of gravitation in Chapter 2. 
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APPENDIX F: 

SOT't FEATURES OF TIE ELECTROMAGTCTIC THEORY OF GRAVITATION 

Having derived the equations of motion for the "post^ ‘^-Newtonian" approx- 
imation of the tensor (Einstein) theory (Appendix, Eq E-18), we observe the 
remarkable similarities which exist between it and the vector (electromag- 
netic) theory of gravitation - primarily the use of a vector potential. 

This sort of correspondence also manifests itself in the matched asymp- 
totic expansion method which decomposes the tensor field equations into 
scalar, vector, and Maxwell stress-tensor components in Eq (3-49). 

The Classical Basis For The Vector Potential In Gravitation Theory 

A simple hydrodynamic argument constitutes the teleological basis fcr 
invoking a vector potential into the Newtonian field equations (2-5) of 
Chapter 2. These equations for $ are irrotational. Hydrodynamically, a 
particle orbiting in an irrotational fluid acquires no intrinsic sp^n. If 
it does, the fluid is rotational, possessing a curl. Now extending this 
argument to the solar system, one is struck by the experimental fact that 
all of the planets possess an intrinsic angular momentum (S) in addition to 
their orbital angular momentum (L). 

Of course, we know that a ponderable, elastic body (with Love numbers) 
which is some form of ellipsoid (formed cosmogonically) will acquire an in- 
trinsic spin S which is commensurable with its orbital angular momentum L in 
the scalar potential $ ^ r of Newtonian theory. This is due to the grav- 
itational gradient, inductive energy dissipation, and the spin-orbit inter- 
action mechanism of, say, tidal friction (Appendices G & S). In other 
words, just because planets spin one does not necessarily need a vector 
potential to supplement Eq (2-5). 

Nevertheless, the gravitational hydrodynamics of the solar system does 
establish an experimental basis for invoking a vei r potential, and thereby 
achieve an electromagnetic theory of gravitation, as in Eq (2-12). The 
vector theory of gravitation in Eq (?-12) becomes even more interesting 
when one considers the origin of magnetic fields and the nature of the 
Coriolis force. 
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Rotation Of Mass & The Origin Of Magnetic Fields 

H.A. Wilson (1921, 1923) points out that there may exist a correlation 
between the rotation of mass and the origin of magnetic fields - basing 
his query upon the experimental fact that the ratio of the magnetic mom- 
ent to the angular momentun for the sun is the same as it is for the 
earth. Wilson also investigates a similar line of reasoning in the elec- 
tromagnetic theory of gravitation by supposing the "relativity" of grav- 
itational fields. In the same fashion that a uniformly moving electric 
charge exhibits an apparent magnetic field, then perhaps a moving neutral 
mass possesses an apparent magnetic field. Wilson then performs an ex- 
periment, but measures no such effect. 

Blackett (1947) pursues Wilson's ideas and advocates thefcrmula 

P=-p£u (F-0 

as the linear relation between the angular momentum U and the magnetic 
moment P (where B is a dimensionless constatn of the order of unity). 
Fuchs (1948) generalizes (F-l) as 

P * i/JiG ^ 

where v Is the gravitational permeability, e is the permittivity, and ug 
is the Bohr magneton. 


The Lorentz Force & The Coriolis Term 

Wilson's conjecture is too naive - as his experimental results determine. 
That is, he considers the analog of the Lorentz force stated in Eq (2-15) 

and surmises that the axial vector b for the "gravitational induction" field 
is a magnetic field B - just as in Maxwell's relations (2-11). 

However, from classical mechanics we already have a force term exactly like 
the Larmor- Lorentz term in Eq (2-15), namely the Coriolis force 


L*. = * z f« x X = -f *»(*«). 


Whence, we can make the identity 


b =-2u 

ftf A* 


> 


(F-3) 

(F-4) 
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whereby the Lorentz force (2-15) becomes 

Result (F-5) consists of the Newtonian (Coulomb) force and the Coriolis 
force. 

A re-interpretation of Wilson's conjecture, then, is to state that there 
is a relation between mass in uniform motion with spin and the Coriolis 
force. Likewise, in the same fashion that one has motionally induced 
electric fields 

E*-vxB (F-fc^ 


for equilibrium states (F=0 or high conductivity), one also achieves 
motionally induced gravitational fields 

<J =-Vx(2«) # 

Eq (F-6b) is the electrodynamic way of stating that "fictitious forces" 
behave like gravitational ones and are both characterized by g in 
(Eq A-l) of the equation of motion 

r i r*yjc-- o . hs> 

It should also be transparent that the electromagnetic theory of gravita- 
tion is Machian in nature, attributing inertia to the inductive effect of 
distant cosmological matter. 

The gravitational vector potential A in Eq (2-19) arises in the presence 
of a hydrodynamic mass current 


j = pv fc-n) 


★ 

Epistemologically speaking, however, the Machian school of thought finds its 
principal experimental basis in electrodynamics. Hence electrodynamics and 
the electromagnetic theory of gravitation stand on their own. Mach's work, 
in fact, is a philosophical derivative of Maxwell's electromagnetism - 
particularly in light of Maxwell's work on gravitation, available to Mach. 
The vector theory of gravitation does not have to answer to Mach. 
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in Eq (2-16). If v is rotational, that is, if 


V = Oxr ^ 

I 

2k) = VxV ( 
** * * J 


(F-7) 


we have directly from (2-19), if v can be moved through the integral sign, 

= -v<b/fc £•§ 

b = VxA = iofl*- M> 

<k (k C * C <■ • 


In order to establish a basis for the origin of magnetic fields in rota- 
ting mass, one should turn to the continuity equation (2-27d). For neutral 
matter (p g =0) , 

V.j = .*.j . ^ 

One can then argue that electromagnetic currents J arise in the presence 
of mass currents j. This sort of reasoning can be used to justify Blackett's 
relation (F-l). 

The Coriolis Vector Potential & Gravitational Radiation 

If we now temporarily adopt Eq (F-9) as a definition of the gravitational 
vector potential £, from which derives the Coriolis force, we can quickly 
surmise the behaviour of the field equations. 

First of all, the energy density in (2-23a) is 
The "gravitational induction" pressure (magnetic pressure) is 

(p e * -Is ) 

while (from plasma theory) the pressure of matter p is related to it by 



no 


? * (f-Kb) 

(p + ■55 * c *"* fe ) 

The cyclotron frequency (eE/m) becomes 
while the Larmor frequency is 


The spin u of the vector gravitational field is precisely the Larmor 
frequency of electromagnetism. 

The ExB drift velocity perpendicular to ^ or b is 



while there is an adiabatic invariant (Coriolis mirror) given by 


sin« _ sinaf 

b 2u 



Const. 


(NS) 


Likewise, a "plasma oscillation" follows in the usual fashion from the 
Maxwell equations (2-12) and the continuity equation (2-17), 




(v sg ) 

The magnetic moment is precisely the mechanical angular momentum L, 

(F-n) 

The torque acting upon a spinning mass of angular momentum L in the Coriolis 
field 2a; is 

T* Lxb (F-i® 



Ill 


The Poynting vector is S = g x h in Eq (2-18), or 

m A# 0m 



where the polarization of S is determined by c. 

Finally, it should be pointed out that the field intensity g is derived 
from (2-1 3a) 

3 = -V§- (2-Oai 

$ is really a centrifugal potential, containing both the Newtonian $ 
as well as centrifugal effects due to rotation of the source of the field, 
then includes the time-dependence of the Coriolis term. 


Arbitrary Vector Potential 


There is nothing implicit in the field equations (2-12) that demands a 
strict identification of the gravitational induction £> with the Coriolis 
effect 2u>, although it should be consistent with it. One can consider a 
ring current j of mass, using the vector potentials from electromagnetism, 
and determine the gravitational induction field b in spherical coordinates 
as a solution of (2-19) (SmythfjP. 294-295), 



where the ring 
For r>a. 


is at some colatitude a, and a is the radius of the mass. 



♦22fi 

<r a 




i-4-4— (MX**') 

Then integrating the current loop over the volume of the spherical mass, 
which is simply the moment of inertia of the mass, one gets approximately 



ii 2 L _ GI usw8 

~ 2 c* r 4 

UCOsd 

v 5- 

This is the tfector gravitational induction field about a mass spinning 
at angular velocity w. 

The orbital angular momentum L of another particle orbiting about this 
spinning mass interacts with b, creating a spin-orbit and a spin-spin inter- 
action. First, there is a precession of the Kepler conic by an amount 





fl L * <b r > 


for a polar orbit. Averaging over the orbit (where I is the Keplerian 

o 

moment of inertia I=Mr /2) 





(F-Mfc 


There is likewise a spin-spin interaction, similar to the gyroscopic effect 
predicted by Schiff (Phys Rev Lett 4, 215, 1960) in post- Newtonian General 
Relativity (Appendix E). If the second particle (satellite) possesses an 
intrinsic spin angular momentum $=1'^, it experiences a torque due to b of 


T * i s*le 

which results in a precession of the spin axis of 

n - T- j. 

“s* * * * t - 


(F-24) 


In an equatorial orbit 

Q s *i<S> - 


Compiling (F-23) and (F-25) 






(F-H) 



which is Schiff's result, except that the first term must be multiplied by -3/2. 
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Generalized Ohm's Law 

The classical relation for Ohm's Law (2-16) generalizes as follows: 

^ =/ > [3 + 5f^l -f** , * ll) 

where represents a pressure of matter. For the Coriolis interpretation, 
this is 

£ ' f [? +y*^ . *+ 


MacCullagh 's Equations & Gravitational Radiation 

It Is pertinent to remark that MacCullagh's (1838) equations can be used 
to derive Maxwell's Equations for electromagnetism. Likewise, they can be 
used in the vector theory of gravitation. MacCullagh's theory, in effect, 
is a hydrodynamic theory of electromagnetism, and in this context, of grav- 
itation. He assumes that velocity v is given by our (F-7) 


^fxv r 2u> 

(F-i> 

where u> ■ His equation of motion is then assumed to be 

• ^ 


II 

1 

X 

• 

(F-29) 

(F-7) and (F-29) are subjected to the zero divergence conditions 
sibillty or transversal ity) 

( incc^pres- 

g.v =0 
H*o . 

(F-30) 

From (F-7), (F-29), and (F-30) one can derive the source-free Maxwell equa- 
tions ( with displacement current in 1838!). If you recall the hydrodynam- 

Ics of current 



(1-IC) 
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then we have the implicit use of .’lacCullagh' s ideas if we simply state from 
(2-1 F) that 



For arbitrary constants of proportionality a, 6, 

Case to V* irff 

t'ift 


(F-30 



we obtain Maxwell's relations (2-11) for p=0 and j=0. We can Identify 

* 


Case (a) €.* i ZfA( 

Cas«Cb) /•* • 



Sources can be readily included. The purpose in going into this dis- 
cussion of MacCullagh's theory is to demnstrate that all of the "magneto- 
hydrodynamic" analogies in the vector theory of gravitational radiation are 
not magnetohydrodynamic at all. Magnetohydrodynamics was established by 
MacCullagh over a century ago. Again, the electromagnetic theory of grav- 
itation goes back to fundamentals. 


A Generalization Of Bateman 's Equations 

Starting with Maxwell's Equations (2-11), Bateman (1915) makes a simplifi- 
cation which in effect Is reminiscent of the original MacCullagh theory. 
Bateman assumes a complex electromagnetic field 

M = E ♦ i H (r-i4) 


whereby Maxwell's Equations reduce to and derive from 



(F-3S) 


These are simply a complex formulation of MacCullagh's Equations (F-29) 
and (F-30), with sources. (E.g. let M = v t i$ ). 
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The neat thing about Bateman's Equations is that they provide a simple 
way of synthesizing vector electromagnetism and vector gravitation, as a 
complex unified field theory. Ltt 


M ■ E * ft 

A» » % 

t = E+iJ 

& * H. 


(F-3C) 


Substitution of the complex fields (F-36) into the Bateman Eqs (F-35) 
gives the vector Maxwell Equations (2-11) for electromagnetism and (2-12) 
for gravitation. The Lorentz a 7£1 S ) has a complex component. 


Or 


F*F+F - i, , f -&F I (r-3B) 


One is Interested only in the real part of F in (F-38). 


The Interesting feature of this approach is the implication that charge 
and mass are the real components of the same thing - a single complex 
quantity. That is, the Newtonian and Coulomb forces are unified into a 
single complex operation. The same is true for the Larmor-Lorentz term. 


The field intensities are determined by 


A* _t 1 

* “ E M 
ft’* HMi 1 


There is no negative energy density. 


t *1 


e- r*A . 


(F-») 
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APPBCIX G: TIDAL FRICTION & THE VELDCI1Y OF PR7RWSATICN 

Eddington (1920) discusses an Important physical application of the 
retarded Lienard-Wiechert potentials when he considers the question of the 
existence of a couple or torque acting on the .wo-body Kepler problem. Re- 
ferring to Figure G-l, Eddington points out that the acceleration a acting 
on the present position of either component must be in the direction of 
the other binary component 



or a resultant torque T acts upon the system. Assuming a phase lag 5 
between the retarded and present positions, this osculating torque is 


where 


T * 

(6-0 

r 2 s * 


•i 
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Note that this is an inductive energy dissipation mechanism (it goes 
as r‘^), characteristic of scalar (Newtonian) gravitation theory. Accord- 
ing to the Lfenard-Wiechert scalar potential 



the force per unit mass or field intensity g = -v<|> in Eq (2-3) (assuming 

*0 

zero vector potential A) acting on either component in Figure G-l is 

3* * ^ (S-$ 

where M=ym for a circular orbit and where is a unit vector directed to- 
ward the present position and not the retarded one r. Hence, 

$=o tn 

and there is no torque. 

However, this equilibrium condition (no couple or torque) derives from 
the fact that the velocity of propagation of gravitational information c is 
just right. If we perturb the velocity of propagation 

c' • c*€ (®*s) 

the gravitational action is no longer in the “present' 1 direction ? 0 , 6^0, 

and a torque manifests itself. The conic then dilatates or contracts until 
once again it is in equilibrium, whereby the phase angle disappears and the 
torque vanishes. These results apply to point masses. 

The first matter that needs to be dispelled is the misconception of Laplace 
(1829) and others that variations in the velocity of propagation e necessar- 
ily manifest themselves as anomalous behaviour of the Kepler conic. This Is 
not true, because the retardation effect above cannot be distinguished from 
the mechanism of tidal friction (Appendix S) in the case of real, oonderable 
bodies. The second matter of importance is the assumption of a constant vel- 
ocity of propagation c in the derivation of the Lienard-Wiechert potential 
(G-2). In the presence of an index of refraction n (n * 1 + 2m/r for light 
in a Schwarzscnild background) whereby c' = c/n, the two-body problem 
expands or contracts until the torque disappears. 
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Hence, the effective inductive friction in the two-body problem due to 
retardation behaves a great deal like tidal friction, except that it 
occurs for point masses as well as elastic bodies with Love numbers. In 
turn, both of these regulate the behaviour of gravitational radiation 
damping (or anti -damping). c urthermore, a non. inear field theory can 
affect the velocity of propagation c', producing retardation, and in turn 
neutra’izing radiation reaction mechanisms. By means of this orbit-orbit 
interaction a nonlinear mechanism can vary the velocity of propagation in 
the vicinity of a "black hole" binary; there can be gravitational radia- 
tion (assuming it exists for the Kepler problem); and there does not have 
to be collapse. 
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APPENDIX H: DIFFERENTIAL EQUATIONS FOR PERTURBATIONS 
OF THE SCHWARZSCHILD METRIC 


The differential equations representing the nonspherical perturbations 
of the Schwarzschild metric 

SR^-srJ^+Sr^ (3»» 

are first derived by Regge & Wheeler (1957) with a number of errors w.'.ioh 
are later corrected by Vishveshwara & Edelstein (1970). They are derr.' t 
using spherical harmonics and are comprised of odd-parity and even-pa. i tv 
solutions. 


Odd-Parity Equations 

In the case of odd (magnetic) parity (-l)*" + \ Eqs (3-56b) take the 
following form in the canonical Regge-Wheeler gauge (3-58): 



where the subscripts (such as v^) represent derivatives (v r =dv/dr). 

Eqs (H-l) are then set equal to zero in order to determine (3-56a). For 
L=l, all of the angular factors in (H-l) vanish. Also for L=1 , fiRg^O. 
For higher L, Eq (3-56a) implies the terms in brackets vanish, providing 
three radial equations. The second-order radial equation resulting from 
(H-lc) is actually derivable from the other two, provided 


(h-2) 


which indeed is the case for the Schwarzschild metric where 



From these other two radial equacions { H- 1 a ) and (H-lb), h Q can be 
eliminated, resulting in a single second-order radial equation 


p + ftfc-vni§ ^ k‘/*^>*o. 


If one defines 


Q = V * (»-*?)* 

r*= r+2mtn|£i-i| O^-sW) 


then from (H-4) one obtains the odd-parity "Schrodinger" equation first 
given by Regge & Wheeler 


£q 

«»p** 


♦W-Xh»)Q-o 



where V e ^ is an effective potential given by 





In the case of the Schwarzschild background 



(ii-i) 

(M-a) 


Even-Parity Equations 

Continuing in the Regge-Wheeler gauge, and using upper-case letters to 
represent the even (electric) parity (-1)* - radial functions in order to dis- 
tinguish them from the lower-case odd-parity functions, the first-order 
perturbations of the Ricci tensor in (3-56b) are 


<H-a4 
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«*BMr-f *i 

s vM* *^*v® -£ *«Hafc 

♦[t^M-i x *i'-( w >lK]^ J r w (**« 


Derivatives of the Legendre polynomials vanish for L=0. For L=l, the two 
angular factors in of (H-9f) are not independent. For L>1 , substitution 
of (H-9) into (3-56a) gives the following even-parity radial equations: 

H«H.*H» (H-toa 

$ (mow 

♦*k(H*K) « O 

<■*** 

•itt4WfO 

-UMK + wV(i-ay l K - ItkrH, * O (H-IOj) 

2»w < z.k* ^-k*H v «ho9 

♦ (v-^)l*s ♦wW$jsO 
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Eqs (10e,f,g), of second order, can all be derived from (10b, c.d), of 
first order, provided 


(H-U) 


+ [iikP-*tfc!&JSL] H|a o. 


This relation (H-ll), in conjunction with the first-order equations 
(H-10b,c,d), gives one second-order differential equation: 

§ 5 (i-a)^ 

X« & X*« p*/tm W-*»V 

DW- V -*(&*' WX$r) • 


(tv-u) 

*0 

(fr04 

(»-«i 

(H-na 


The second-order even-parity differential equation for S=H^/r is then 

} 4 g. Oh*) 


Further details are provided by Edelstein & Vishveshwara (1970). It is 
safer to follow their development as it does not contain the errors made 
in the original Regge-Wheeler (1957) paper. 
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APPENDIX I: REGGE-WHELER GAUGE 

Recall from Appendix B that in Einstein's field equations for the 
first-order perturbations, h^ remain invariant under the gauge transforma- 
tion 



Odd- Parity Gauge 

Upon consideration of h yv in (3-57a), Regge & Wheeler let 

°= V = o 

S** 

whereby (B-5) transforms (3-57a) into (3-58a). A is arbitrary. 
Even-Parity Gauge 

Similarly, if for arbitrary M Q , M-j , and M one selects 

\ ~ MM&YjVo 

3 S « 

then (B-5) transforms (3-57b) into (3-58b). 
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APPENDIX J: EVEN-PARITY "SCHPDDINGER" EQUATION 


Zerilli's (1970b) derivation of the second-order "Schrodinger" equation 
for the even-parity Regge-Wheeler perturbation problem is as follows. He 
substitutes kR=H^ in the first-order even-pari.y equations (H-10b,c,d) 
and the algebraic relation (H-ll), and then uses (H-ll) to reduce (H-10b,c,d) 
to 


_ 



where a,B,Y> and 6 are functions of r,L, and M - but not k=w. 
following transformations are assumed u A 

K = 4fr)K * j(r)R T 

R= £(r)l? + J 

Jr s Mr) Jr* =. (i-W^Jr* . 


Next, the 

(») 

(X-3> 


n(r) in (J-3) is Eddington's (1920) refractive index. 

t ^ _ Lx^0r* ♦ 3X»r ♦ 4 m % \ 
+(f) * r 1 tXr;^) 


JW* » 

t(r>* i 


t-Xf x » SXmr t 3 w*T 

l?^ 55 (Sfo 3 r 

-ir* 

(r-2m) 


where 

X= 


Zeri 1 1 i 


lets 




noting that throughout the L and M subscripts on all radial functions have 
been suppressed. By virtue of (J-4), (J-5), and the integral of (J-3) 
which is (H-5b), then 



These may be re-written in second-order as 

* [If- VL»lK . - 0 


V L fr) 



(X+Ot^-V^XNnr 1 ♦ llWr + t6m J 

r»(Xr*3m)* 


q 


where angular momentum L and azimuth M have been re-instated. 


(J-7) 

(J-8) 
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APPENDS K: TENSOR SPHERICAL HARMONICS 

The following are the orthonormal set of tensor harmonics in 4-space as 
developed by Zerilli ( 1 970a ,c ) : 


00 

a s 

■*in 


Xi ° 0 ° 

/l\ • 

[oV° 1 


oo 0 o' 

o o o 0 


%,° 0 0 

a = 

OV ® 

0 o O 0 

> |F 

O O 0 0 

> ^Ul 

0 o o o 

. o o o o. 


.O 0 O 0 . 


.0 0 0 0 . 


.(•) . » 
by.-"* 


1L 






0 o o o 

• o o o 

*00 o 


b in -riW 


00 0 o 

0*0 0 


O 0 *«•«, 

0 o o o 

* o o o 

* o o o 




0 o o o 

o 

0*0 o' 

0*0 o 


* 4 

iH 


0 o 0 o 


ooo o' 

0 0 o o 


o o o o 

0 0 m^ih 

Q 

1 * 

00 1 O 

OO • tMlXg, 


,0 o o *»• 


■*in 


i* 


0 o o 
o o o 
o o 
0 o 


o 

o 




where 






a, b, and £ are orthogonal to the 2-sphere (o,4>). and £ are tangent. 

f-(e + h) and g«#(e - h) are, however, orthogonal. Thus an 

orthonormal set. Upon generalizing to a 4-dimensional pseudo-Euclidean space, 
we obtain 4 more tensor harmonics: a^ , Jb^, . Regge & Wheeler 

use e in place of f. Hence, theirs is not an orthonormal set. 
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APPENDIX L: COEFFICIENTS OF SOURCE DECOMPOSITION 

Any synmetric covariant tensor can be expanded In terms of the tensor harmon- 


Ics in Appendix K: 


ST = 

rr <0 0) A ril#V*> . o . 

H | 

M 

The coefficients A LM ,..., are determined by the inner products ^ 

► 

t 



(L-J) 

where 

T-.S = yV^.S* 

(L-3) 

and n is 

yv 

the background metric. Accordingly, 
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Now, the coefficient of a^ is , J, ). So from (L-2) through (L-4) 


Thus, from Appendix K 


\,«*-JjY>T, JO . <«> 

But from Appendix N, if the Schwa rzschild geodesic postulate Is assumed. 
Hence, letting y = dT/dS° (and similarly for the other coefficients): 


(L-8) 


' L rf * 
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APPENDIX M: DIVERGENCE CONDITION FOR SOURCES 

The divergence of the source stress-energy tensor 6T yv in Eq (3-61) must 


be zero. Since ( fi T )’ v Is a vector, expanding it In tensor harmonics gives 

ST*. £• 
/ 1 

■.IVl-WWt} 

Because of (3-61) 


(n-0 

then V L,H 

o 

II 

> 

> 

tz' 



*ui = Pm ~ " K* * 0 . 

(rvi) 

A simple calculation gives the following divergence conditions: 



(n-i) 


S m - i(^QSr o . 


DIVERGENCE CONDITIONS 

Upon Integration of these conditions (Appendix C of 2eri Ill's thesis) one 
obtains the equations of motion for a point particle falling along a geodesic 
of the Schwarzschlld geometry n - demonstrating consistency with the 
geodesic postulate. 
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APPENDIX N: TRAJECTORY OF PARTICLE & STRESS TENSOR 


The stress-energy tensor due to Peters (1966) in Eq ( 3- 62 ) is 


T /V = m. J 

where is the invariant delta function defined by f/f/S^ (x)/(-g)d^x - 1. 
s is the affine parameter along the world line z(s) of the particle. In 
the case of the Schwarzschild metric 



-V ▲ i_ 4 
’U r»wi0sr swiO 

S < 4 > (x) * S(x*)S(s) . 



Carrying out the integration (N-l), using (N-2): 



where ft : (0,4>) and 5^(ft) = 5 (cos0)6 (<J> ) . 

R(t ) , e(t), 4>(t), and T(t) are the position of the perturbing particle in 
the Schwarzschild coordinates r,e,$, t. 


This is the very critic tion in the P, „ers-Zeri! M analysis which 

adopts the Schwarzschild geodesics c.s t ho path of the radiating parti :le - 
in order to evaluate T as stated in (N-3). Radiation reaction is thus 

y v 

inaccounted for. 



131 


APPENDIX 0: RADIAL EQUATIONS WITH SOURCES 

The following radial equations are a consequence of the Zerilli tensor har- 
monic decompositions (3-64) and (3-65) of the Peters field equations (3-60): 

Magnetic-Parity Harmonics 

-X'lWoS 

v (k-0(t«4(r-ua Quj 

6-¥)§ - 0-9) 1 & * m K * xir’fitL^aul Ju . 

Electric-Parity Harmonics 

-0-Wi£ - (<- 

«,* fC 

*{&)$’ xAu, 

-&-??$ ♦*-*«&♦(>-?&$ -o-h^ >2^ -o-t)^ 

i(H.-H>-xr 1 lit^H'‘F'w, . 
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APPENDIX P: FOURIER TRANSFORMED FIELD EQUATIONS 


We may write the Fourier transform of the field equations in Appendix 0 as: 


Magnetic-Parity Equations 

«>2iw -XftMuof r'Qftv* 

* X ihUfOd-iXtn)! 1 r # # 


Electric-Parity Equations 

**»* AV* 

= -XiWLfc 4 *rB**r> 

»“H. * $ V 

+ ^ - r ('*f^ 4 S* -ir L 0 -» 0 (W* « -MRtGjp) 

t(H^- H, J * -i*Ja.wu- mVFJjv) t 
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APPENDIX Q: THE AUXILIARY RADIAL FUNCTIONS 


Having determined the auxiliary radial functions R^ as the solutions of the 
"Schrodinger" equations (3-66) and (3-67), one can derive the original Regge- 
Wheeler radial functions and determine h^ wit!, the following relations: 

Magnetic (Odd) Parity Pe^tui .Lions 



(o->) 

(Q-») 


Electric ( Even) Parity Perturbations 




(C) 

in 


in 


M 

(H 


U - T Xr (r-lm)- u A „ UMWrM .. £ 

°lm L (r- 2 m)(\r \ 3m) ]f\n [ior(Xr»3»)y ‘ua ut 

H n„ s H ** + . 


M 

(CHI 


Auxilia ry Source Functions 


C = r Xr* [j L(i.»ny * 1 pH) _ / tr \ ft 
l rn L iu) r-lm ] U Lr. \r-zm J *\n • 


(O-i) 

(0-9) 



APPENDIX R: THE BONDI 96 FUNCTION 


Bondi's representation of gravitational radiation begins with his inves- 
tigations into negative mass in General Relativity (1957a) - negative mass 
originally being considered by Foppl (1896) and Schuster (1898) although 
being discussed again by Einstein & Infeld (1949) - and his rediscovery of 
certain plane-fronted waves (1957b, 1959). He publishes in 1960 the first 
results of his multipole approximation method, and then in 196.2 in conjunc- 
tic- with van der Burg & Metzner. Bondi establishes a radiative solution for 
an isolated axisynmetric system whose mass decreases monotonically due to 
outgoing radiation. 


Bondi 's Metric &. The " Main Equations" 
Bondi's metric For axial symmetry is 


4s = (Vr ’ AuV**)!.* ♦ 2«*4»ir 
♦ 2Br** , *J»49 -•‘*(t**4®* ti > w944 < ) 


for coordinates (u=t-r, r, e,$). That is to say. 




UrY* o 

1*^ O O O 

rV* O -rV* O 

o o o 


«- 0 


(R-2) 


o o o 

q^ v -VjW* O (^-3) 

3 = O TIi* O Wi> 

O O O *i*^ t *i»6 # 

Bondi assumes throughout a stationary-radiative-stationary transition 
as depicted in Figure 6 where the radiative solution is represented by 
Einstein's field equations in vacuo, 


— 0 


(3-S^ 


just as is assumed in the Regge-Wheeler formalism. Using arguments based 
upon a contraction of the Bianchi identities 




- A 


(ft LA\ 
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with g ac i.i Eq (R-3), Bondi reduces the field equations R i y ■ 0 to the 
"main equations" * R^ = R^ * R 33 s 0 : 

Bondi 1 s " Main Equatio ns" 

0 = -tr'Ra* - 2r^-< tt »2*J[ t -2fcr H -2l 1 <*e| (RSI} 

0 = R a ^®-rR, 1 t'= 2V, ♦ r*Ua- 4»U,-i i ‘U 1 «-tO (**) 

Ox-lC/r - irfrWL +<y*m- K»WV-r(>-»J,Kt t V( 0 *e-| 1 )tt l (WU) 

+|-(ZrVlJ,»r2«t8-Wt«)U ♦ 

plus two supplementary conditions (Appendix 2 of Bondi, vd Burg, Metzner) 

^R«=0 (RO 

and the trivial equation R 01 " The various subscripts in Eqs (R-5) 
correspond to ordinary derivatives: ( 0 , 1 ,2,3)->-*-(u,r,0,<t>) . The Eqs (R-5a) 
and (R-5c) are referred to as the "hypersurface" equations because they in- 
volve only differentiation in the retarded hypersurface u=t-r=const, while 
Eq (R-5d) is known as the “standard" equation. Only the standard equation 
contains a derivative with respect to u. 

Numerical Integration Scheme 

Bondi's iterative solution of the "main equations" (R-5) consists of 

(R-7) 


1 - c.(u,81^ . 
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The functions of integration 

H(u,o) , uu,9) , n (u,e) , c(«,9) (R-® 

provide all the information or "news" as Bondi calls it, about the be- 
haviour of the source. The functions (R-8) can be reduced to one single 
"news function." 

Multi pole Expansion & Sonmerfeld Radiation Condition 

Based upon causality arguments, Bondi assumes an expansion of I in his 
metric (R-l) as 

V. t jfij) + . . . (*.*) 

■ r 

in terms of retarded time u = t - r, which is equivalent to an outgoing 
(Sommerfeld) radiation condition. 

Boundary conditions require that L(u,o) = 0, eliminating one of the func- 
tions in (R-8). Hence, the leading terms in the iterative scheme (R-7) are 
given by to. _ 

y 8 c(£» + e^... (ft-tta) 

P s H(u,6) •• • (R-< 

U * ZH 1 i? ,, /r ■*• ••• (R-lOcJ 

V = «'x*[l+2H J dt6+4H|+lH u ')+... (W04) 

Under coordinate or gauge transformations which do not change the Bondi metric 

(R-l), H(u,o) can be eliminated. Furthermore, by virtue of the Somnerfeld ra- 

2 

diation condition, Bondi argues that g(t-r) in Eq(R-9) or c(u,o) in (R-iOa) 
must vanish. Thus, he arrives at the following canonical form of the "main 
equations," noting that c = c(u,o) - and not the speed of light: 

Bondi ' s Canonical " Main Equations" 

J » c(*,e)/r ■*••••• (R-*W) 

TJ * - (c t * 2c uteW ♦ [iNCMHttc** 4 c 

•orcM) T 

\i =. r- iflcov- - cl-Acvote - tc (u***0]r 1 

4C,» 
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The News Function c(u,q) 

On the basis of the canonical form of the "main equations" (R-ll)» the 
supplementary conditions (R-6) reduce to 


Supplementary Conditions 

M 0 = -c\ + i (c„+ 3c t «*0-2c) o 
-3N„ = li 1 +3cc M +4cc o o40 + c # c t # (R-ab) 


If H and N are given for one value of u (that is, on some hypersurface 
u = const), and c = c(u,g) is given, the entire solution is determined. 

For this reason, Bondi calls c(u,q) the "news function." Another appro- 
priate name might be the "data function." 

Correspondence of Bondi's equations with a static, non-radiating case re- 
veals the physical nature and meaning of the functions of integration M, N, 
and C. An empty, axially symmetric metric such as Bondi's (R-l) can always 
be reduced to Weyl ' s form 


where 



*• * [l + ] 



(R-Oai 

(Sow) 


Letting a represent the transformation parameter between (R-l) and (R-l 3) , 
defined by , / . 

c'-T&t- <»«**) , (R-wl 

then one obtains (static case) 


Bondi's Monopole 
Bondi ' s Dipole 
Bondi's Quadrupole 


Mot,®* m / 

N (U.8) = DsinG - 

C (u,6)* iQunG - ^DsinG + ) 


where m, D, and Q are the mass monopole, dipole, and quadrupole respectively. 
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Depletion Of Mass Monopole &. Radiation 

M(u,e) in the "main equations" (R-ll ) , based upon the correspondence 
to the static case (R-15), is called the "mass aspect." The mass of the 
system is defined as the mean value of M(u,o) over the sphere 

= 2 ; J rl(u,9) sinOJS OMft 

Note that for the static case M Q = sM/^u * 0 and (R-16) gives m(u) = m, 
as it should. Differentiating (R-16) with respect to u and using (R-12a), 



That is ^ 

m ® = " * J C o S,n ® 

or re-instating the derivative 

The mass m of the source or radiating system, then, is constant if and only 
if there is no "news" : 



If c o ^0 on some hypersurface u=const (Figure 8), then m Q is negative and m 
decrees monotonical ly. 

Outgoing radiation depletes the mass of the source, then, according to this 
formalism. The news function c(u,o) of Bondi is generalized to include another 
polarization c(u,o,<t>) by Sachs (1962) while cylindrical gravitational news 
is addressed by Stachel (1966). 
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APPENDIX S: GRAVITATIONAL RADIATION & TIDAL FRICTION 

The inductive friction due to retardation and variations in the velocity 
of propagation in Figure G-l of Appendix G must be contrasted with tidal 
friction, which constitutes spin-orbit interac.ion: 



It occurs in the earth-moon system and is due to imperfect elasticity with a 
resultant energy dissipation which characterizes real matter. Neglecting the 
eccentricity of the orbit, the rate of change of expansion of the conic 
is approximately 


a T ~ [3A 1 J|R, 5 sm2|S|]a' }4 (s-») 


where the Love umber k 2 


is defined in terms of the inertias C and A as 



3G (C-fl 



Tidal friction (S-l) has the opposite effect as does gravitational radiation 
(that is, if gravitational radiation collapses the binary), which is given in 
Eg (5-5) as 



rus 1 
IT c r 


Mm(M«n)] 




The critical value of the apse a c where gravitational radiation and tidal 
friction remain in equilibrium is given by equating (S-l) and (S-3): 


_ cl[ is AwnllU lH R! 

C ■ G l U m . 


(S-4) 
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For a binary system of two solar masses, for example, this is 

<•-» 

or 

a c = 2. 25 X to ,4 [Ji sin R* . ($-6) 

An optimization of ksin2 1 6 1 with = 3/2 and 6 = ir/4 gives 

<l c * 128 R* (Sri) 

where R Q is the mean radius of the primary in Figure S-l. 

A binary system collapses under gravitational radiation until it reaches a 
size of the order of (S-7) and then it attains a stable orbit, radiating away 
the angular momentum of the primary component until the latter's angular 
velocity a, in Figure 5-1 is exhausted. 
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APPENDIX T: SCALAR, VECTOR, & TENSOR SYNCHROTRON MECHANISMS 

The work of Davis, Rufflnl, Tiomno, & Zerllli (1972) as well as Breuer, 
Rufflnl, Tiomno, & Vlshveshwara (1973) analyzes the power emitted from 
a particle In a circular orbit r = (3 + 6 )M about a Schwarzschlld "black 
hole." The treatment Is general, so the vector case can readily be modi- 
fied to treat vector theories of gravitation as discussed In Chapter 2. 
The even (-1)^ or electric and the odd(-l)^ +1 or magnetic components for 
the power In Eq (5-12) are respectively: 


Scalar 



Vector 



fm 

odd 


Tensor 

X* itt-OOM u** i 


The Green's function terms are 


where 


I 9 lunubiuii tci mo Ml c - - — m M ^ _ 

$gw>- Mk So: i**** ra ♦*«)»* 

€f«m ■ U4f ♦ MS > *r* S’W 
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APPENDIX U: THE PETROV-PIRANI CLASSIFICATION 


The Petrov types in the Penrose diagram break down as follows: 


f Cylindrical Waves (Einstein&Rosen, Rosen, Marder, 

I Bonnor, Weber & Wheeler) 

* | Axisymmetric field of Weyl (1918), Levi -Ci vl ta (1919) 

(Most General ) ^ * Empty space-times 

Isolated, bounded sources (oscillators) (Bonnor, 1963) 
All real fields (Piranl, 1962) 
k Field of an actual Isolated system 


I 

D 

nr 

N 

(Radiation) 

(Flat) 0 


J Robinson & Trautman (1960) 

[Peres (1959) 

I Schwarzschild Metric 
Robinson & Trautman (1960) 

( Semi-far fields 
Robinson & Trautman (1960) 

Plane Waves (Brinkmann,1925; Rosen, 1 937 ;Takeno, 

Kundt) 

4 "Plane-Fronted" Waves (Robinson, 1956) 

Any radiating system of sources 
Far radiation zone fields 
Robinson & Trautman (1960) 


Flat space-time 


It should be apparent from these results of the Petrov-Piranl classifica- 
tion that the method is much too general. 
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APPENDIX V: GRAVITATIONAL RADIATION & COMPOSITE FIELD THEORY 

One can consider the linearized graviton in the context ' composite field 
theory. The composite nature of electromagnetic quanta is suggested by 
Jordan (1928) using statistical arguments, de Broglie (1932, 1934, 1936) 
suggests that a photon is composed of a neutrino and an anti -neutrino, 
establishing the basis of the neutrino theory of light. This is one form 
of a more general Idea that elementary particles are composite particles 
built up out of a fundamental Spln-1/2 fermion, such as the neutrino. If a 
Spl n- 1/2 neutrino and a Spi n-1/2 anti-neutrino couple together and form a 
Spin-1 boson, It is a photon (emission or creation); and if they mutually 
couple and cancel out they annihilate (absorption of a photon). In terms 
of hole theory, a neutrino jumps the energy gap of the neutrino sea, crea- 
ting a neutrino-anti neutrino pair (emission of a photon). 

Jordan (1935) originally treats the de Broglie postulate in terms of a 
four- component neutrino theory, using a neutrino of momentum k and an anti- 
neutrino of momentum £ - k (producing a Raman effect for neutrinos) However, 
this four- component work Is brought to an abrupt end by Pryce (1938) when 
he demonstrates that it is not Invariant under spatial rotation o» arbitrary 
Lorentz transformation. Perkins (1965), however, overcomes this short- 
coming In the Jordan-Pryce treatment and establishes an interesting four- 
component neutrino theory of the photon, formulated from two two-component 
Weyl equations. Perkins' four-component model Is Invariant under spatial 
rotation, although his resultant photon operators are not strictly Bose 
commutation relations, due to additional terms. However, Planck's radia- 
tion law still follows as it does from Bose statistics. 

From the vector (Spl n-1 boson) electromagnetic theory of gravitational radia- 
tion, then, there Is a neutrino theory of the vector graviton - in direct an- 
alogy with the work of Jordan, de Broglie, Pryce, and Perkins. One can like- 
wise construct a neutrino theory of the graviton comprised of four coupled 
neutrinos or two coupled photons, for the linear Soln-2 tensor theory of 
gravitational radiation. This implies, of course, that linear gravitons can 
decay into four neutrinos or two photons - e.g. upon interaction with matter 
or absorption. Wave resonances due to interaction between gravitons and 
photons both travelling at tfr same velocity of propagation (Gertsenshteln, 
1962) are easily visualized from the aspect of composite field theory. Care 
must be taken, however, to account for the Interaction of photons with 
charge and gravitons with mass. 
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One can speculate, furthermore, in a candid way with the bold conjecture 
that given a fundamental Sp 4 n-l/2 quantum (such as a Spin-1/2 boson) he 
could construct a unified field theory using composite field theory. If 
this quantum were a nonlinear Spin-1/2 graviton, for example, he could ex- 
plain e-decay, y-meson and n-meson decay in terms of gravitational radiation 
'•• J ithout need for the neutrino. Weber's high fluxts might then be cue to 
u-m&son decay in the earth's atmosphere. 

There is another reason for pursuing this line of reasoning. The creation 
and annihilation of photons is completely compatible with the Dirac treat- 
ment of hi^ negative energy states, namely with positron-electron pairs. 
However, Dirac explicitly assumes that the infinite sea of negative energy 
exhibits no gravitational effect (if it does, we have Maxwell's postulate 
of an infinite sea of gravitational energy, in 1865). But if gravitation?! 
radiation exists and if we are to have a consistent quantum theory of radia- 
tion, then there must be a gravitational dual of Maxwell-Dirac spinor elec- 
trodynamics - for gravitons. Such a model may be afforded by neutrino 
theory, or by invoking (negative) mass conjugation rather than or as well as 
charge conjugation in the Dirac theory. 
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APPENDIX M: SPIN-2 & GRAVITATIONAL RADIATION 


The sixteen (16) independent variables h in Einstein's field equations 
(3-1) and (3-6) must be coupled to a terror source T . Generally speaking, 
a tensor T yv consists of Spin-2, three Spin-1, and two Spin-0 admixtures: 

FL,]-* (W-l) 

(D (®) (e) fo) H W * It 

where the numbers in parentheses under £q (W-l) represent the various spin 
components or helicities of each spin group. If T is symmetric, 

V <** 

then (W-l) loses six (6) degrees of freedom or two of the vector Spin-1 
admixtures and reduces to 



The spin parts are 

ITH 

[T/] -v I • 0 ftWtt 

[T./’l-*' © (***> 

The dot is placed in the zero components of (W-4b) and (W-4c) to distinguish 
them as the same tensor Spin-0 contribution, distinct and different from the 
scalar Spin-0 in (W-4a) which is coupled to the trace of T . Note that (W-l) 
through (W-4) are not equations but group- theoretical relations. It is (W-4s) 
which causes inconsistencies in massive radiation theory (Appendix X, Eq X-32), 
where T^jK) provides an excitation of the Spin-0 degree of freedom. 

If energy is assumed conserved, then 
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This means there is no source to which the Spin-1 and Spin-0 components in 
(W— 4b) can couple. Energy conservation, then, causes these helicities to 
vanish, reducing the degrees of freedom by four (4) - eliminating namely: 

+1, 0, and 0. (W-5) also dictates that 

0 (w-0 

is true. (See Fronsdal, 1958, on higher spins.) 

Symmetry ( h uv =h V( .) and energy conservation have thus reduced the 16 un- 
knowns to six (6). Invoking the Lorentz gaige condition 

*y/=o ( w - 7 ) 

on h ^ eliminates four (4) more v "vector" gravitons) and reduces the degrees 
of freedom to two (2). These are the pure massless Spin-2 helicities, with 
all other spin admixtures removed. This result constitutes the so-called 
"transverse-traceless" gauge (Appendix Z). 

In massive radiation theory for Spin-2 (Appendix X), on the other hand, 
the Lorentz (also called Hilbert-Lorencz) gauge condition does not produce 
four additional constraints, but rather it reduces the field equations in 
conjunction with energy conservation to provide a relation between the trace 
n = h p ^ and the trace T = T u ^ of the energy-momentum tensor. This gives 
only a fifth constraint, thereby reducing the six (6) independent variables 
to five (5), or 2S+1, where S is the spin. (See Eq X-28.) 

In addition to the above group- theoretic representation of the gravita- 
tional Spin-2 problem in terms of reducible groups, we can recall the tensor 

spherical harmonics of Appendix K. Spin-0 is comprised of a^, and a, 

(01 (01 9 W 9 

Spin-1 of and c,, and Spin-2 of d, <j, and This representa- 

tion, however, is often confusing because the respective harmonics consist 
of lower-spin admixtures and are reducible. 
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APPENDIX X: MSSIVE GRAVITATIONAL RADIATION 

Massive radiation theory is an embodiment of the screening potential e“ mr 
due to Laplace (1846). It concerns Itself with wave equations such as (2-7), 
(2-25), and (2-32) which include a finite rest-mass in the form of Kleln- 
Gordon equations (1\*c*l): 

(□-««?)$ »“V 

(□-m?)A* * -xJ* (X-0 

(o - O A'* = -xT** 

The subscripts on the masses m,. (i * 0,1,2) correspond to Spin-0, Spin-1, 
and Spin-2 for the respective scalar, vector, and tensor wave equations. 

Radiation theory with quanta of finite rest-mass has been well developed 
In electrodynamics, based upon the de Broglie postulate (1924) of a finite 
photon mass and its subsequent incorporation into Proca's vector meson 
theory (1936). The resultant massive electrodynamics is well-behaved since 
It reduces to Maxwell's theory as the photon mass goes to zero. 

However, such is not the case with conventional treatment of massive grav- 
itational radiation theory using Einstein's General Theory of Relativity. 

An assumption of a graviton rest-mass does not reduce to Einstein's theory 
In the massless limit. Indicating either that (i ) Einstein's theory is 
particularly unique, or that (ii) conventional representations of massive 
radiation theory are inadequate and bolder hypotheses may be justified. 
Einstein's theory is said to have no "neighbors," meaning it does not tolerate 
neighboring theories with a finite rest-mass. 

Massive radiation theory is the contemporary expression of the cut-off 
potential of Laplace, later investigated by Neumann (1874, 1896), Seeliger 
(1895), and Yukawa (1935). That is, an assumption of a finite rest-mass m Q 
for the scalar wave equation in (X-1) modifies the classical Newtonian po- 
tential with an exponential cut-off: 



In this fashion, the Newtonian or Cnulombic field is screened and does not 
have an infinite cross-section, an artifice also employed in plasma theory. 
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Massive El ectrodynamics 

Massive electrodynamics is important because it is manifestly covariant, a 
feature which is not true for the massless Maxwell theory. It has been 
developed by many authors, among them Wentzel (1943), Kallen (1972), Jauch 
& Rohrlich (1955), Sjorken & Drell (1965), and Bogoliubov & Shirkov (1959). 

de Broglie (1924) has demonstrated that if one assumes the photon to have 
a finite rest mass mj, then the energy relation 



necessarily implies that the vacuum has a refractive index 



and behaves as a dispersive medium. This idea is important because it is 
the basis of Proca's extension of Maxwell's equations for massive quanta. 


The massive Maxwell -Proca equations take the form 


where u = m 0 c/1i. These 
to a conserved current. 


- £ 4 > 

V*i 


VB =0 



equations represent a massive vector field coupled 




i£r 

C 






Their wave equation is 
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with a Poynting vector 

(x-8) 

and an energy density 

. l x ' 9) 

The addition of the massive quanta to the Maxwell theory increases the 
number of degrees of freedom in the field due to the fact that the Lorentz 
condition no longer reduces these degrees of freedom. Instead, 

r (O * ? r„ («-»> 

relates current conservation and the Lorentz condition. In addition to the 
two transverse polarizations of the photon, there are now a longitudinal 
photon and a scalar (time-like) photon due to the presence of a photon 
rest mass. (See Kaller..) 

From all of the above equations it is apparent that the massless limit 

7 7 . 

u * m-j * 0 reduces to the Maxwell theory. However, massive radiation 

theory is not as simple as this because the most significant relation above 

Is (X-9) for the energy density. The massive energy density must be 

positive semi -definite or everything will collapse into any negative energy 

★ 

states that happen to exist. 

For this reason, massive radiation theory is usually developed using an 
action principle. Provided the action I is expressed in a canonical 
"pq - H" form, 

I = * f(pi-»)A 

then the Hamiltorian energy density If can be readily identified in the Lag- 
ranglan density X and hence examined for negative energy states. 

*Th1s problematic feature of Lagrangian field theory can be alleviated by 
a massive exclusion principle, discussed in the final section of this 
appendix. 
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The first-order Proca action coupled to a prescribed conserved source J° 

* m V +J \] OmO 

has been decomposed by Deser (1972) into its transverse and longitudinal 
parts, demonstrating that the longitudinal mode constitutes a scalar field 
decoupled from the current. The energy density is 

t; • (*sY • "« H (it * iY-W} 

+$[»'(*)'+ &(7| -rf] . (*I2) 

The "T" and "L" superscripts represent "transverse" and "longitudinal" 
respectively. The second term in braces is non-NaxMelllan in that it re- 
mains coupled to gravitation in the m=0 limit, which may have implications 
about the frequency dependence of the gravitational deflection of light 
(See e.g. Woodward & Yourgrau, 1972). 

Massive Gravitational Radiation 

The current dilemna in the theory of massive gravitational radiation can 
be best characterized in terms of the source theory (Schwinger, 1970) 
associated with the Feynman amplitudes of quantum electrodynamics. One 
can represent the gravitational interaction as the exchange of a Feynman 
propagator between one source T pv (X) and another T otB (X* ) : 



FIGURE X-l 
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This propagator is given in a momentum space representation by 




p -nn -ie 


where Tl^ v ^ is known as the spin projection operator, 
by Ogievetsky & Polubarinov (1965). 



These are stated 


The resultant propagator for a massive gravitjon is (van Dam & Veltman, 
1970, 1972): 



which reduces in the massless limit (m=0) to 


if* 0 . 

>W/» ' p* 


In contrast, the graviton propagator for the massless case is 





(X-«0 

(X*) 


Comparison of the second term in the numerator of (X-15) and ( X-l 6) 
demonstrates t: 't the two theories are incompatible. If we take the case 
where the delta tun*. Jons are all unity, the numerators of (X-15) and (X-16) 
are respectively 



4 


and 





This means that 


T) n *° _ 3. T) mr0 



whereby the deflection of light 6 is changed as 


^ massive 



massWss 
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the massive radar time delays are likewise 3/4 of the massless theory, 
and the precession of the perihelion A6 is 


AS 




matsiw 


•t* 


nuasless 



The Jordan-Brans-Dicke (Brans & Dicke, 1962) scalar-tensor theory gives 
similar results when its dimensionless parameter w=0. 

KMtan . ( JH9 ^ 

w=«oin (X— 1 9) gives the Einstein theory. The Eddi ngton-Robertson metric 


comparatively speaking is 

a=B=Y=1 Einstein (Massless) 

a=e=l y s H Einstein (Massive, m=0) 

a=B=l y=(u+ 1 )/(u+2) Brans-Dicke 

for these three representations. Drawing further upon the u)*0 analogy be- 
tween the Brans-Dicke and the massive theory, the massive Eddi ngton-Robertson 
metric in the massless limit (m=0) is: 


lx = 0- W firf) (**»> 

t 


This is not the Schwarzschild solution. 


Massive Gravitational Action 

Many of the problems associated with massive gravitation theory have been 
addressed by van Dam & Veltman (1970, 1972), Zeldovich (1970), Boulware & 

Deser (1972a, 19 7 b), Deser (1973), van Nieuwenhuizen (1973), Vainshtein 
(1972), Sexl (1967), Freund, et al (1969), and Oglevetsky & Polubarinov (1965). 
Massive two-tensor theories have been investigated by Isham, Salam, & 

Strathdee (1971), Aichelburg et al (1971), Alchelburg & Mansourl (1972), 
Aragone & Chela-Flores (1972), Al chel burg(1973) , and Lawrence & Toton(1972). 



The definitive work at present has proceeded with the Hilbert-Palatini 
formulation of the action principle (Arnowltt, Deser, & Misner, 1962) 
which represents the full Einstein action as 
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](V> ■J(l'j ij -NR--H,Hyx ,M> 


where Latin indices vary over 3-space and where is the conjugate 
momentum: 


H «(>)** Hi -- 30; R* * 

-fC * *R ♦ 3* («*-*!»*¥«> 


(X-2la) 


The linearized massless action (including the coupling 
>sh T wv ) is 

MV ' 


•J 

to tru 


e source 


T l M = J A ft - (t% * W ♦ 


hyT' 


+HiT" + *»wr].(*»> 


In this action (X-22), and h Q() constitute Lagrangian multipliers - upon 
whose variation are recovered the auxiliary condition* (h»h^) 

V*k - = -T* 


>01 


(x-23a) 

(X-13W) 


IV -iT 

■j 

Alternatively, one can vary with respect to N and in (X-21) and recover 
the constraint 


R^=0 


M 


Hence, the field equations following from (X-21) or (X-22) consist of 
16 unknowns which are reduced to 10 upon the assumption of symmetry 
(hyv*h )* These 10 are reduced to 6 by energy conservation 
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T“> 0 

(X-2S) 

and further reduced to two (2) by the auxiliary conditions (X-23) or (X-24). 
Thence, the massless graviton has two (2) degrees of freedom (±2 heliclty). 

Now, the Fierz-Pauli (1939) massive action I m A^A 


*:] 

can be added to (X-22) to give a total massive action 


I- VI- ■ 

OWT) 

The immediate consequence is that variation of and h Q0 in (X-27) 
longer gives the massless auxiliary conditions (X-23) but rather 

no 

Vk-h,j +T“ . „*h 

(X*28a) 

2 TTV +n?M.--T" . 

(X-28b), as a constraint, simply recovers itself 

(X-28« 

Nj = « a (2T ryr) 

(X*28c) 


and does not eliminate three (3) other variables. Hence only (X-28a) 
constrains the six (6) degrees of freedom subsequent to energy conservation 
in (X-25). (X-28a), then, reduces the degrees of freedom by one, to 

five (5). The massive theory has five (5) or 2S+1 helicltles associated 
with a Spin-2 graviton of finite rest mass. 

Orthogonal Decomposition Of Massive Action 
Boulware & Deser (1972) have considered an orthogonal decomposition of h^j 
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which In turn decomposes I In (X-27) Into a ,: transverse-traceless" (Iyy). 
a "vector" (Iy), and a "scalar" (I s ) part: 

i - 1„* vi 5 . (**» 

This Is further decomposed by rescaling the action in order to eliminate 
problematic m~ c singularities and to get it in a canonical "pq - H" form: 

i ■<♦(«?♦ «♦«♦*) . (**> 

I c represents instantaneous Coulomb contributions, while the superscripts 
correspond to spin components or helicity. All contributions in (X-31) 
uncouple and vanish as the graviton mass nv+O except the first and last 

I(m=o)= I^+lf . (*32) 

Hence, the massive theory retains a scalar component coupled to the trace 
T° a of the energy-momentum tensor (See Appendix W). This massive theory, 
then, is a scalar-tensor theory. It is the Spin-0 or scalar contribution 
which creates the Inconsistencies. The right-hand term In (X-32) is that 
which derives from the trace in (W-4a). 

Speculations On The Problems Of Massive Limits 

The literature has investigated ways around the incompatibilities which 
arise In massless limits, using such techniques as indefinite metrics, 
indefinite probability, broken symmetries, Goldstone bosons, and the 
cosmological term (Appendix Y). As remarked earlier, bolder hypotheses 
may be warranted due to the inadequacies of the conventional massive theory. 
Therefore, some additional approaches are stated here. 

(A) One can suppress the Spin-0 helicity coupled to the trace T p ^0 by 
employing both tuvanced and retarded Green's functions in the determination 
of the graviton propagator In Figure (X-l). This Is tantamount to removing 
the Sommerfeld "adlatlon condition for the Spin-0 exchange, and constitutes 
a rejection of the Feynman propagator. 

(B) The Bel-Roblnson tensor has zero trace In empty space. It is possible 
that a massive gravitation theory (following Einstein) could be developed 
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around it, guaranteeing that no source can couple to the problematic Spln-0 
component. This tensor is fourth-order, but such a theory would be based 
upon its contracted, Spin-2 form. 

(C) Recalling the gauge conditions discussed Appendix W for the Spin-2 
problem, a new gauge prescription may provide possibilities. In 2q (W-5), 
the conservation of energy 

T'V 0 

forces the Spin-0 helicity of T uv to be zero: 

T'VO . W 

In other words, the Spln-0 contribution in the tensor source T is de- 

yV 

manded to be zero. Instead, this tensor Spin-0 component could be used 
to cancel out the scalar Spin-0 component coupled to the trace T^T* 1 ^. 

That is, let 

T-T^O . 

This gives a new gauge prescription 

"■(t; r; m 

instead of the Hilbert-Lorentz gauge studied by Cgievetsky & Polubarinov 
(1965). In their notation, the spin projection operator for Spin-0 

is broken up into two Spin-0 admixtures Pgcal ar anc * ^tensor IJS ^ n 9 
Hilbert gauge (q is arbitrary) 
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such that 



p% p“ - 1 


Owfc) 


The suggestion here is to use (X-34) and suppress P 


( 0 ) 


enti rely: 


whereby 


pC°) 

r s r " r r 


0«4J) 

(**) 


The price for this approach is that T uv 


is no longer conserved: 


* 





(D) Lastly, one can speculate on the existence of negative mass (Wilson, 
1972). In the Stuckel berg-Feynman representation of Quantum Electro- 
dynamics, the positron in a negative energy s.ate is equivalent to a 
positive energy electron travelling backward ir :ime. This view is oased 
upon the fact that changing .he sign of charge _ and time s in the equation 
of motion (X-36) does not. change the electromagnetic behaviour of the 
particle: 


m 


m/i-7* se #*F^ 


(X-34) 


However, the equation of motion (X-36) is also not changed under a mass-time 
reversal. By attributing the characteristic of negative mass to anti -matter, 
then, a particle of negative mass in a negative energy state behaves like a 
particle of positive mass in a positive energy state, electromagnetically 
speaking. 


Hence, by adopting an exclusion principle based upon intrinsic rest mass, 
only particles of negative rest mass can occupy the negative energy mass 
shell In Figure (X-2). Rest mass cannot change mass shells (If m Q =0). 

*Rastall (1972) treats the case T v . =*R . 

' Vju ,V> 
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E 



FIGURE X-2 


As a consequence, a Lagrangian such as (X-26) and (X-27) can have a 
Hamiltonian energy density containing negative energy states and not under- 
go collapse of positive matter into t l ose states. Instead of the Fierz- 
Paul i .ssive action (X-26), we can adopt the following one 

T„ - $J'x(Z) 

where 

t n H*) - i(h* »»H) (x-37) 

which gives the correct limit as m+0, namely Einstein's theory. Although 
the right-hand term in (X-37) is negative, no collapse of particles 
occupying positive energy states can occur into negative energy states, 
due to the exclusion principle invoked above. Problems associated with 
"ghosts" (as particles associated with these negative states are called) 
are eliminated. 

A more rigorous treatment of the above conjecture must treat gravitational 
behaviour under mass-time inversion. This invariance My be surmised from 
the radiation reaction equations of Lind, et al (1972, Eq 11). 
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APPENDIX Y: THE COSMIfEGtt. TERM £ GRAVITON REST MASS 


The field equations (3-1) of General Relativity have numerous cosmological 
implications among them gravitational collapse and an expansion of the 
Universe. Concerned about the latter, Einstein (1917) has considered a 
more general form of the field equations which still satisfy the contracted 
Bianchi identities and energy-momentum conservation: 



or in terms of the Einstein tensor G , 

yv 

+ ^ V y r * ^ ^ 


The most general spherically symmetric solution (by virtue of Birkhoff's 
theorem) of (Y-la) is the exterior Schwarzschi 1 d form with A^O, namely 




which for m=GM/c^=0 is 


a* A t* + jl Ar* + 

- !-ia- A.p* s * x 



is = AA* . (f*W) 

Because photons following geodesics in (Y-ld) do not travel at the speed of 
light c, A/0 implies a photon rest mass. Furthermore, the orbital equation 
is 

^ + a = -fa * 3«u - (*|^ 

and an additional precession of the pe>~!apsis « 

The cubic e v -0 now has two positive roots: r=2m and r*/(3/A). For a not to 
effect Mercury's orbit by more than one arc sec per century, 
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Teleologically speaking, the additional "cosmological term" governed by 
the "cosmological constant" a is meant to offset the expansion of the Uni- 
verse. However, Hubble's (1929) interpretation of the cosmological red- 

* 

shift as a Doppler effect has been invoked as a reason for rejecting the 
term, maintaining that General Relativity predicts a cosmological expansion. 


Nevertheless, a meaningful basis for retaining the cosmological term may 
be found in the problems associated with the theory of massive radiation 
(Appendix W). In such a context, the question before us is whether or 
not the cosmological constant A in (Y-l) is equivalent to or can be related 
to a graviton of non-zero rest mass. Some say yes (E.G. Tonnelet, 1965; 

Peak, 1972: Kurdgelaidze, 1965; and Freund et al, 1969) while others say no 
(Treder, 1968; Pol ievktov-Nikol adze, 1 967 ) . At the same time, some main- 
tain that A cannot alleviate the inconsistencies of the massless limit in 
Appendix W (Boulware S Deser, 1972). 

The wave equation for gravitational radiation on a non-flat background 

containing the cosmological term follows from the formalism of Peters (1966), 

Isaacson (1968), and Zerilli (1970) employed in Chapter 3. The variation 

** IQ) 

of Eq (Y-l) for a stable background n =g v ' is the following: 

(«-2A) - -wST,, . 

This equation can now be simplified by defining the function fT (this is 
the same es ^ in App. B, (B-7), and (3-10) except that the background 
is arbitrary) 

v * v*y M 




See, however, Eddington's (MTR, 1922, p. 154) Interesting argument for A^O. A 
provides a fundamental length. "An electron could never decide how large it 
ought to be unless there existed some length Independent of Itself for It 
to compare itself with." 

★* 

Stability of the background n, „ must be taken for granted in order that 
6T can be assumed small. 
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and Its divergence 



Re-grouping terms and substituting ( Y -3) and (Y-4) In (Y-2) v* get a 
relation similar to (3-60) 

. CM) 

Now we Impose the Hilbert (Elnstein-deDonder In Appendix B) gauge which 
sets (Y-4) to zero 

f sf >v = 0 tf-0 

r 7* 

and suppresses the vector gravitons. (f^O can be retained for further 
simplification in some cases of n uV t although problematic negative energy 
states may be associated with these degrees of freedom.) Wave Eq. (Y-5) 
hereby reduces to 

+^(R*2A)*-2xSTJ,. (Y-7) 

In an empty ( T uv = 0)» Ricci-flat (R^=0) space with no cosmological constant 
(R»4 a* 0), (Y-7) reduces to 



which Is the starting point of the Regge-Wheeler-Zerllll formalism In 
Chapter 3. 
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However, since Aj*0 we know that the field equations (Y-la) demand that 


4A - R - XT 

(M) 

whereby (Y-la) becomes 




For an empty Universe (T =0 and T=0), (Y-10) reduces to 

U*‘ 


and (Y-9) to 

(Y-to) 

R = 4A 



Substitution of (Y-lla,b) into (Y-7), using (Y-3), shows that the con 
tributions due to Aj*0 are now of second order in h . Neglecting these 
terms (particularly if A is very, very small ), (Y-7) simplifies to 


(HO 


Note that one can arrive at (Y-12) to first order in h^ v by using g yv as 
a raising and lowering operator rather than the background n wv - a result 
which leads Treder (1968) to the conclusion that a terms cancel out of the 
gravitational wave equation. Furthermore, note with caution that (Y-12) 
and (Y-8) are not the same wave equation. 

Overtly, the cosmological terms have vanished from (Y-12), like (Y-8) 

whoro A=n. but the character of the Riemanr • sor R a 6 Is significantly 

Mv 2 

different in the two relations. For a space - constant curvature K=l/R , 
the cosmological constant A^O is still manifest. That is, the Rlemann tensor 

*Tolman (1934) places a lower limj,t a£ 10 " 5 ^ cm“ 2 ^10“^ (light years)" 2 
upon A, which gives a ratio (Ar 2 /3)/(2m/r)*»10 * 9 at Neptune's orbit. 
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\ii r+K ^3 / a'3o3 / y) (t-o) 

reverts to 

for use In (Y-12). This substitution (raising and lowering with n yv ) into 
(Y-12) now gives a K and a A term contribution 

■ 2K [^' V** 

(\45) 

to second order in h^ v . Recalling that K is related to A by 

K = A/3 , (Y-lO 

(Y-12) is to first order 

iy,/ - JAiy M 


Noting that 

h * h (•* *»)) (*») 

then a traceless gauge h = 0 means that 

h =0 h*o «• > 1 ^ 2 . (f») 
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Whence, (Y-17) reduces (n a °V2) to 

-2X$T >«» 

r 

■ 

in a traceless-Hilbert gauge: 

(Y-6) 

0H9) 




At this point it is necessary to determine if the wave equation (Y-20) can 
be put in a Klein-Gordon form 



and determine if the A^O term in (Y-20) gets cancelled out. To simplify the 
computations, assume dn=0 and note the e" v -*-l and e v -*l, as r+0. We wish to 
determine if h yv . a ,a goes to a flat-space d 1 A1 emberti an In a locally flat 
region of (Y-22). We can write 



Define 
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where — 


Dh„« h 


L y* 




*/* % d (Y-Ma) 

-WW-SW'’- m 

V * -tyX - <&)% (vm.) 

B is the term of interest. A and C contain terms of second order, 

uv uv yv 

or terms that vanish in locally flat space (r«l). Furthermore, only the 

first-order second derivatives in B^ v remain as r+0. These terms are 


which we can define as 


®<^V = + tyu> 


where 


6 r**'*f fyrt/1y+ v/ lyel 

Next we note in this approximation ttat 7 


□*-?S?+V T ' 

D V^ 7 V“ _2A/3 ‘ 

Dn..-»-y a h..T3?+2A/3 . 


C«5c) 

(Y-25J) 

(*2Se) 

m 


D*J,| ^ ^711 7 {o + 2 A 

The results In (Y-26) are not ±2A because of the type of expansion in (Y-24). 
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We find that 


4 (M7a) 

whereby (all other terms do not contribute) 

Nex F »« = . Cf«c) 


, (f«d 

whereby (all other terms dc not contribute) 

€b, s *iAR,i 6 ,,-MHi. Gj-'f-Ah,, , Cf2 8W) 

And lastly, 


whereby ^ “0 

M \A.K 

Substituting (Y-29b), (Y-28b), (Y-27b,c) back 
G + H =0. That is, B* = F : 

Vjv yv yv yv 

8l * - HL 


WW, 

<**> 

into (Y-25b), we find that 

<~|Ah. . W-M> 


Substitution of (Y-30) and A =C =0 into (Y-24) determines that 

yv yv — 

(□- $h)Ko «-2HjT„ 

-*A)h„ =-2H5 T h . 

The d'Alembertian in (Y-31) is not arbitrary because curvilinear (e.g. spherical) 
coordinates cannot be used. Curvilinear transformations have been lost in the 
expansion (Y-24). What is important is that the sign of in (Y-30) and h 
in ( Y-20) is the same, with no cancellation. 




The significant result is that there is not a constant term hidden In a decom- 
position of h^.^’ 01 which removes the potentially massive term 2A/3 b In (Y-20). 
The cosmological term A^O does apoear to constitute a massive contribution to 
the linearized gravitational wave equation (Y-20) in a traceless-Hilbert gauge. 
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APPENDIX Z: TRANSVERSE-TRACElfSS GAUGE 

Upon consideration of the nature of the spin admixtures inherent in 
tensor radiation theory (Appendices W & X), we know that the linearized 
massless theory of Einstein (1916, 1918) and Fierz & Pauli (1939) is a pure, 
transverse Spin-2 with no longitudinal Spin-0 scalar contribution. This 
scalar Spin-0 couples to the trace of T , that is T«T W , as in (W-4a). 
Hence, the massless graviton h yv and its energy momentum tensor must 
be traceless in empty space. 

Recalling the Regge-Wheel er perturbations of the Schwarzschild metric 
(3-57) or (3-58), however, it is not readily apparent that h^ v is traceless. 
It Is, on the other hand, true for the plane-wave discussed in (3-33) and 
(3-37). Of course, as in electrodynamics, such a transversal ity condition 
holds only in a source-free region of space-time. 

The transverse- traceless (TT) gauge, then, is that coordinate condition 
(B-l) whereby the solution 



has zero trace. That is, 

h =h TT = 0. (2-1) 


Furthermore, it is transverse in the sense treated for the plane-wave of 
(3-24), (3-33), and (3-37). 

Note that in the transverse-traceless gauge, y yv of (3-10) and (B-7) and 
the Hilbert function (Y-3) are identical with h^ v , except for possible 
scale changes. The Hilbert-Einstein-deDonder coordinate condition (B-8) 
and the Hilbert gauge (Y-6) are likewise simplified. 
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Including a current J, (F-35) is 



Bateman's Equations, reminiscent of classical six-vectors and quaternions, 
provide a simple way of synthesizing vector electromagnetism and vector grav- 
itation, As such they are the teleological basis of a complex unified field 

th rv et * r-N- Vt 


6 2 h 


r m 

T iT 


(F-36) 
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^ ‘ «L ^ <SSft *«RP 38 18 B J» || 

Substitution of the complex fields (F-36) into the Bateman Eqs (F-35' ) gives 
the vector Maxwell Equations (2-11 ) for electromagnetism and (2-12) for 
gravitation. The convention 

* \*t _ rt . i*i (F-3 7) 


r e */ , e[ 0 +i « x <3 ] 


provides the complex Lorentz force relations, 
the real parts of F in (F-37 and F-38). 


EtiuM IZ 

One is interested only in 


The interesting feature of this approach is the implication that charge 
and mass are the real components of the same thing - a single complex 
quantity. That is, the Newtonian and Coulomb forces are unified into a 
single complex operation. The same is true for the Larmor-Lorentz term. 


The field intensities are determined by 

£= k Z - AV G* = (E 1 * H 4 ) (gS h a ) . (F*39) 

The above formulation is not derived from a Lagrangian, however. 


A Quaternion Formulation Of Maxwel 1 Vs Equations 

Recalling the rule of multiplication of two quaternions A=a+/^ and B=b+J3, 

AB s - A B *§* a§ * b ft + (F°40) 

then Bateman's representation of Maxwell's Eqs is simply 



M» E+iH 
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a*<L *6 * J m (F-41) 
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where 
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A Unified, Vector Field 


A derivation of the Maxwel 1 Eqs. (2-11 ) foil ows from the Euler-Lagrange Eqs 
for a variation of the Lagrangian density 

To get the gravitational Maxwell Eqs. (2-12) one must change the sign of either 
the first or last term in (F-42) . However, reversing the sign of F^F yv gives 
a negative energy density in the radiation field; hence, one must change the 
sign of the current coupling term But in order to derive the Lorentz 

force (2-15) and get a Newtonian rather than a Coulombic interaction, the 
signs of both of the last two terms in (F-42) must change. This argument, 
however, then results in the wrong sign for the middle, kinetic term. 

Nevertheless, if we adopt a complex vector potential 


where 


iy (y (f- 43) 


and where 


> r* 






results in the electromagnetic Maxwell Eqs (2-11) 

S" * 3 




and the gravitational Maxwell Eqs (2-12) 
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as well as the corresponding Lorentz foi 


Tf 

equations 
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(F- 44 ) 

(F-45) 

(F-46) 

(F-47) 

(F-48) 


The current for M 


-i(M* rf' 4 - (f* 49 > 


is conserved and does not transport charge or mass. The equations are linear 


as in the Maxwell theory. Interchanging M 


and M in (F-49) reverses the 

M 

tational interaction. To make 


direction of the current, basic for the gravitational interaction. To 
the theory nonlinear (transporting charge and mass) one simply adds the 
complex Klein-Gordon charge-mass terms (where p = P e + ip m ) to (F-45) 

VP I _*M*M>* * £ A 
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